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1.0 SUMMARY 


A new gas temperature sensor was developed that shows promise of suffi- 
cient ruggedness to be useful as a gas turbine temperature sensor. The sensor 
is in the form of a single crystal aluminum oxide ceramic, ground to a cone 
shape, and given an emission coating. A lens and an optical fiber conduct the 
thermally emitted light to a remote, near-infrared, photodetector assembly. 
Being optically coupled and passive, the sensor is highly immune to all types 
of electrical interference. Candidate sensors were analyzed for optical 
sensor performance, heat transfer characteristics and stress from gas loading. 
This led to the selection of the conical shape as the most promising for the 
gas turbine environment. One uncoated, and two coated, sensing elements were 
prepared for testing. Testing was conducted to an indicted 1750° C in a 
propane-air flame. Comparison with the referee optical pyrometer shows an 
accuracy of ± 25° C at 1700° C for this initial development. One hundred 
cycles from room temperature to 1700° C left the sapphire cone intact, but 
some loss of the platinum 6% rhrodium coating was observed. Several areas 
for improving the overall performance and durability are identified. 


2.0 INTRODUCTION 


Current gas temperature measurements in the hot flowpath of gas turbine 
engines has been limited by material properties and cooling techniques to 
approximately 1300° C. Recent advances in optical methods and refractory 
materials showed strong promise for extension of this temperature limit of 
the sensor itself. The use of fiber optic coupling was encouraged to provide 
the system with greatly improved immunity to electrical interference from 
nearby lightning strikes, radar pulses, and other noise sources. 

The program was planned to design, fabricate, and test an optical temper- 
ature sensor capable of measuring gas temperatures to an upper limit goal of 
1700° C. The program was carried out in four tasks: 

Task I - Preliminary Design 

Task II - Design and Fabrication 

Task III - Testing 

Task IV - Reporting 



3 . 0 PRELIMINARY DESIGN 


The various elements that made up this first task are shown graphically 
in Figure 1, "Preliminary Design Outline." 


3.1 BRAINSTORM SESSION 

A technical meeting was held to "brainstorm" various configurations and 
techniques to measure gas temperatures to 1700° C in a gas turbine engine. 
This resulted in the generation of additional sensor and configuration con- 
cepts beyond that described in the proposal. See Appendix A, "Brainstorm 
Session for New Sensor Ideas." 

A second result of this meeting was a list of key words and phrases that 
were used for an engineering literature search and for inquiries to vendors. 
See Appendix B, "Key Word List - Optical Gas Temperature Sensor." 


3.2 LITERATURE SEARCH 


From a review of approximately 700 abstracts, approximately 40 full 
reports and several patents were ordered and studied. The abstracts of the 
more important items are included, along with comments and criticism, as 
Appendix C, "Annotated Bibliography." Two subject areas are included: Gas 
Temperature Measurement and Refractory Materials. 


3.3 EVALUATION CRITERIA 

A criteria was developed for the purpose of placing both a value and an 
importance factor on each of the goals (parameters). The candidate sensor 
techniques were rated on a scale from zero to four on each parameter and 
value. 

Table 1 is the final version of the chart that was mutually worked out by 
the customer and contractor. Based on current and future gas turbine engine 
requirements, three temperature ranges were deemed of interest. Temperatures 
below 600° C are usually sub-idle and are of less importance than 600° to 
1700° C. Future engines may use combustion discharge temperatures up to 2230° 
C. Although not included in the original contract goals, this elevated 
temperature range was included because of its future significance. In each 
range, the two levels of accuracy reflect both the original contract goals and 
a second, more realistic goal based on state-of-the-art temperature standards. 

Two time responses are listed. Parameter 13 is a control system figure 
based on usual turbine discharge thermocouple characteristics. Parameter 14 
is a much faster (0.2 seconds) response goal that could be used for an engine 
stall evasion and recovery control loop. On other parameters on the chart. 
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Figure 1* Preliminary Design Outline 



Table 1. Evaluation Criteria - Parameter Value. 





Importance: 
Highest = 10 

No. 

Parameter 

Value 

Lowest = 0 

1 

Temperature Range, Subidle 

170® to 600® C 

5 

2 

Accuracy 

± 1® C 

5 

3 


+ 10® C 

5 

4 

Temperature Range, Idle to 

600® to 1700® C 

10 

5 

Full Power 

± 1® C 

7 

6 

Accuracy 

± 2® C 

10 

7 

Temperature Range, Future 

1700® to 2230® C 

5 

8 

Engines 

± 3® C 

6 

9 

Accuracy 

± 10® C 

5 

10 

Resolution 

± 1® C 

10 


Linearity* 

Repeatability* 

Hysteresis* 



11 

Environment : Oxidizing , 

Mach No. 0.5 

6 

12 

Oxidizing 25 Atmos. 

Mach No. 0.2 

10 

13 

Time Response 3 Sec. for 

63% of Step Changes 

10 

14 

0.2 Sec. for 63% Step Change (Stall 



evasion/ recovery) 


8 

15 

Endurance: 100 cycles RT 

to 1700° C 

10 

16 

Installed Life 

30,000 hours 

5 

17 


1,000 hours 

10 

18 

Source/Detector & Elec. 

Room Temp. 

10 

19 


-54 to 121° C 

5 

20 

Induced Engine Losses 

0.05% Maximum 

6 

21 

Low Risk to Engine (Foreign Object Damage) 

8 

22 

Vibration Endurance 


6 

23 

Low Installed Cost 


4 

24 

Low Operating Cost 


4 

Each 

of these items is included 

under the "umbrella" 

of Accuracy 


and will not be separately rated. 



two or more value listings represent either a different engine application or 
a contrast between the ideal and the practical application of the gas temper- 
ature sensor. 

Table 2 was applied to each of the sensor candidates and for each 
parameter and value. To be considered a candidate optical gas temperature 
sensor, the technique had to use a fiber optic cable(s), with connectors, for 
signal and power transmission. 


Table 2. Evaluation Criteria - Candidate Techniques Rating. 


How Well Are Parameter and Value 

Satisfied? 

Rating 

Full Satisfied 

(Excellent) 

4 

Satisfies Over Most of Range 

(Very Good) 

3 

Comes Fairly Close to Value 

(Good) 

2 

Barely Useful 

(Marginal) 

1 

Cannot Come Usefully Close 

(Poor) 

0 


The candidates’ final score was based on the sum of the products of the 
value and the rating . This is more fully explained in the section titled 
’’Evaluation of Candidate Techniques.” 


3.4 VENDOR EVALUATION CRITERIA- AND EXPERT CONTACTS 


The companies in the sensor market were reviewed. One company (Mikron 
Instrument Company) markets a fiber optic coupled radiation pyrometer, 
intended for surface temperature measurement in protected areas. Another 
company (Luxtron Company) sells a fiber sensor system that uses a phosphor 
coating on the tip of the optical fiber. The upper temperature limit is 
200° C. Finally, Dr. Ray Dils was just launching the Accufiber Company to 
manufacture and sell the optical fiber thermometer product that he developed 
while at NBS. 

Dr. Kenneth Kreider an Dr. Martin Reilly at the National Bureau of 
Standards were contacted. Dr. Kreider was continuing the work on development 
of the sapphire crystal thermometer to two areas: molten metal temperature 
determination (by immersion) and surface temperature measurement by direct 
contact to the surface. 

Dr. Reilly was investigating the sapphire crystal thermometer for use as 
a replacement for the platinum thermometer as used in precision melting point 
furnaces . 
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Helpful literature was received from Dr. Warshawsky of the NASA Lewis 
Research Center on the physics of gas temperature measurement. 

Several sources were helpful in the collection of single crystal aluminum 
oxide (sapphire) properties. Mr. Peter Warren of the Electronic Materials 
Division of the Union Carbide Corporation supplied a very extensive set of 
literature and referred the writer to Professor Arthur H. Heuer of the 
Department of Metallurgy and Materials Science at Case Western Reserve Uni 
versity. Professor Heuer supplied considerable data on the high-temperature 
strength properties of sapphire. His recent laboratory investigations confirm 
the highly directional flow stress phenomena that stem from the rhombohedral 
structure of the crystal. This information has allowed a "worst case" stress 
limit to be found for a particular orientation of the crystal. 


3 . 5 ANALYTICAL STUDIES 


Two types of gas-temperature measurement methods have been used in the 
industry: those that measure gas temperature directly from a thermodynamic 
property and those that measure indirectly from a secondary effect such as 
the heating of an immersed object (thermometer). The physical and chemical 
effects that underlie these two types of measurements are quite different and 
will be briefly summarized below. 


Direct Methods 


Sodium Line Reversal - An early direct method, sodium line reversal, is 
still used today as a standard. The method requires a source on one side of 
the gas and a spectrally sensitive detector on the other side. The method 
offers the feature of gas-path averaging, but it requires the presence of 
sodium (or other alkali metal salt) in the gas (Reference 2). The method has 
even been used in a probe to narrow the view to a localized region in a hyper- 
sonic tunnel (References 3 and 4) . 

Gas Molecular Band Radiation - This method is sometimes referred to as 
the emission absorption technique. It estimates temperature by measuring 
radiation in the optically active bands of the common gases, such as CO 2 or 
H 2 O, which are present in combustion gas. Long-path averaging is a feature 
of this method (References 2, 3, and 5). 

Carbon Particle Emission - This method is similar to the gas molecular 
band method, but is useful over a broader range of the visible/ inf rared spec- 
trum. A disadvantage is that modern fuel-efficient combustors do not generate 
enough carbon particles. Some unique applications of this method have been 
developed (References 6 and 7). 

Raman Scattering - Raman scattering has been demonstrated as a powerful 
research tool for flame diagnostics and combustor modeling. The method is 
suited to medium- or short-path averaging and can be tailored to microsecond 
transients or several seconds averaging. The equipment is expensive and deli- 
cate at the present time, and performance is limited in the presence of par- 
ticulates (References 8 and 9). 
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Velocity of Sound - Measuring the velocity of sound has been demonstrated 
as a potential method of long-path, gas-temperature measurement. The hardware 
is fairly rugged, but accuracy is highly dependent on knowledge of the gas 
constituents, which will depend on the local fuel-air ratio (Reference 10). 

To summarize, direct methods of gas-properties measurement involve some 
intrinsic property of the gas or some seeded tracer element or compound in 
the gas. The advantage of these methods is that no equipment must be main- 
tained at the high gas temperature; they are truly noncontact methods. The 
difficulty is usually in the complex equipment needed to generate and/or 
detect and process the signals. 


Indirect Methods 


Indirect methods involve the immersion of an object in the gas flow, 
attempting to cause the object to "reach" the gas-stream temperature and then 
measuring the temperature of the object by one of several means. 

Thermocouple - These devices are widely used (References 1 and 11). In 
the best probe designs they are placed in a holder to provide thermodynamic 
conditions that reduce the object-to-gas-stream temperature difference to an 
acceptable value. 

Other Sensors - Other devices have used properties such as thermal 
expansion of bimetal elements or (recently) of the Fabry-Perot cell (Reference 
12), the velocity of sound waves (Reference 13), and the velocity of light in 
an optical fiber (Reference 14). An optical absorption technique using rare 
earth doped glass optical fibers has been demonstrated (Reference 14) . The 
most accurate and reliable method, potentially, is the measurement of thermal 
radiation from a near-blackbody cavity. This is the technique used to inter- 
polate and extrapolate the international temperature scale (References 15 
and 16). 

The advantage of indirect methods is in the relative simplicity of the 
hardware to acquire and process the temperature information. The difficulty 
(challenge) lies in two parts: first, setting up and controlling the heat- 
transfer conditions from gas to sensor and, second, accurately measuring the 
sensor (object) temperature. 

Considering the first challenge, the temperature attained by an object 
immersed in a stream of hot gas generally is not identical with the tempera- 
ture of the gas; instead, it is characteristic of a steady-state condition at 
which the rate of heat transfer to the object is equal to the rate of heat 
transfer away from the object. The steady-state difference between object and 
gas temperature is commonly called an "error," but actually it represents the 
balancing of six well-defined phenomena: 
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1. Heat transfer from the gas stream to the object by convection 

2. Heat transfer from the object to the combustor walls or probe 
shield (and vice versa) by radiation 

3. Heat transfer along the object by conduction 

4. Conversion of kinetic energy to thermal energy at the boundary 
layer around the object described by the recovery factor, y) 

5. Heat transfer to the object by chemical reactions on the surface 

6* Heat transfer from the gas stream to the object by radiation from 
the gas stream. 

These processes are illustrated in Figure 2, 

Analytical studies of the heat transfer environment of various probe con- 
figurations were made. A computer program, normally used for gas temperature 
thermocouple data reduction and analysis, was found to be suitable for analy- 
sis of any cylindrical probe in one of two possible orientations to the gas 
stream; that is, either end-on flow or sideways (crossflow). 

The program calculates the difference in temperature between the gas 
stream and the probe tip, taking into account convective heat input and radi 
ation gains and losses, as- well as conduction loss. A complete copy of the 
program is included in Appendix D. 

The first case selected for analysis was that of an unshielded, simple 
cylinder shape placed in a flowing gas stream in crossflow. The gas temper- 
ature was given as 1682.5® C. The mass flow parameter was varied. The cylin- 
der diameter was also varied in binary steps, maintaining approximately 17 
diameters of length. Figure 3 shows the probe orientation and computer 
results. Note that the temperature difference is dependent on both mass flow 
and diameter. For any given mass flow, the smaller probe "reads closer" to 
the true gas temperature. This is considered beneficial since the smaller 
diameter probe also has lower aerodynamic blockage, and it would also do less 
damage if it broke off during engine use. Note also that the smaller diameter 
probe has a smaller change in the computed temperature difference for a given 
change in the mass flow parameter. 

Past experimenters have employed two different diameter probes, exposed 
to the same gas conditions, in an attempt to calculate the true gas tempera- 
ture, given only the probe-indicated temperature, when the mass flow is 
unknown. However, with known mass flow, the need for two probes is reduced. 

Figure 4 shows the probe in a cross section of a gas turbine combustor 
to better illustrate the radiation view factors that influence the computed 
temperature difference. The probe receives radiation from the very intense 
primary zone of the combustor, taken in this example to be at 1927® C. 
However, the view angle factor is only 0.05 of a sphere. Wall radiation 
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Figure 3. Calculated Temperature Difference Between Sensing Probe and Gas Stream. 


] ,1 111 I I .1 11 I J I ) 




represents a much larger view angle factor of 0.80, but the wall temperature 
is assumed to be at 949® C. The wall would consist of the combustor liner 
and the high pressure turbine nozzle assembly. The remaining view factor of 
0.15 represents the probe's view of itself; that is, the stem. This view fac- 
tor has a neutral (zero) contribution of net radiation exchange. The environ- 
mental and assumed parameters are listed in Table 3. 


Table 3. Environmental Parameters for Figures 1 and 2. 


Gas Flow Temperature 

1682.5® C 

Gas Properties Used 

Air Table 

Primary Zone Radiation Temperature 

1927° C 

Probe's View Angle Factor of Primary Zone 

0.05 

Wall Radiation Temperature 

949“ C 

Probe's View Angle Factor of Wall 

0.80 

Probe Emittance Factor 

0.30 

Probe Thermal Conductivity 

4 Btu/hr ft® R 

Calculated Temperature at 

Tip 


The above analysis shows clearly the effect of mass flow on the tempera- 
ture difference between sensor and gas. It was also important to look at this 
temperature difference as a function of throttle position on a hypothetical 
engine. This is plotted as Figure 5, where a unique mass flow was used at 
each gas temperature, simulating a typical gas turbine engine cycle at throt- 
tle positions from idle to full throttle. Table 4 shows the parameters used 
at the two throttle extremes. At intermediate throttle settings, the gas 
temperature and the mass flow were each linearly interpolated 


Table 4. Parameters Used at Two Throttle Extremes for Calculating 
Temperature Difference Between Probe and Gas Stream. 


Engine 

Average 

Mass 

Flow 

Wall 

Primary 

Zone 

View 

Factor, 

Throttle 

Gas 

Ib/sec 

Temp. 

Temp 

Primary 

Setting 

Temp ® C 

/Ft. 2 

° c 

® C 

Zone 

Full 

1700 

169 

871 

1927 

0.05 

Idle 

520 

25 

315 

1281 

0.02 


It can be seen from Figure 5 that the temperature difference is a 
strongly nonlinear function of gas temperature. 
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Figure 5. Calculated Temperature Difference Between Sensing Probe and Gas Stream 
Versus Gas Temperature. 





3.6 LAB TESTS 


No lab tests were done in Task I, Preliminary Design. 


3.7 EVALUATION OF CANDIDATE TECHNIQUES 

The Brainstorm Session, Literature Search, Vendor and Expert Contacts, 
and Analytical Studies were all brought together for the purpose of selecting 
a set of the most valuable and promising candidate temperature measurement 
techniques. Five candidate temperature sensors were chosen from the general 
field of sensors. All techniques appear to be usable with fiber optic cable 
for signal and power transmission. 

Two of the five candidates are of the direct type, in that a temperature- 
related property of the gas stream is directly measured. No physical probe 
is immersed in the gas stream. Instead, the sensor for the two direct methods 
receives infrared radiation from the combustor discharge region of the engine. 
The two direct sensor candidates are: 

♦ Gas Band Radiation Pyrometer 

♦ Particle Radiation Pyrometer 

Three of the five candidates are indirect types, which involves the 
immersion of an object in the gas flow, with the expectation that the object 
(probe) will "reach" the gas stream temperature. Some intrinsic property of 
the probe is then used to measure its temperature. For the three selected, 
this would be visible and/or infrared radiation. The three indirect sensor 
candidates are: 

♦ Hot Target Aspirated Ceramic 

♦ Hot Target Unshielded Metal 

♦ Sapphire Crystal Thermometer 

Each of the five candidate techniques will be discussed in detail. The first 
five characteristics listed under each discussion may be compared between each 
of the candidate techniques. 


Gas Band Radiation Pyrometer 

This gas temperature sensor uses infrared emissive radiation from the gas 
itself. Either the molecular CO 2 or H 2 O infrared bands are sensed. This is a 
direct-type sensor, in that no physical probe is immersed in the gas stream. 
See Figure 6. Infrared radiation is collected by the lens and focussed on a 
fiber optic cable. The cable conducts the energy to the Signal Processor and 
Display where a detector converts the energy to an electrical signal. Appro- 
priate spectral filtering must be used in order to minimize the interfering 
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Figure 6. Gas Band Radiation Pyrometer and Particle Radiation Pyrometer. 




radiation of metal engine parts. A purged sight tube is used to maintain lens 
cleanliness. Characteristics of the technique are: 

♦ Direct Measurements 

♦ Noncontact 

♦ Uses infrared emission from CO2 and H2O bands 

♦ May be applied as a monochromatic (single-band) or as a poly- 
chromatic (two or more infrared bands) system. 

♦ Straight path averaging characteristics, but limited by 
absorption. 

♦ Background metal surfaces present problem of biasing the 
measurement. 

♦ Cool boundary layer may also cause bias. 

♦ Species concentration and pressure level will affect spectral 
emittance and require calibration unique to engine model and 
sensor installation. 

♦ Longwave nature of infrared radiation may require a special 
fiber optic. 


Particle Radiation Pyrometer 

This gas temperature sensor uses infrared radiation from particles in the 
gas to measure gas temperature. There are two main sources of particles: 
externally and internally generated. The level of particles existing in the 
air inducted into the engine tends to vary with time and location. The in- 
ternally generated particles would be carbon from the combustion process. 
Figure 6 shows the basic components needed: a protected lens, a fiber conduit 

for the collected radiation, and a signal processor and display. A potential 
problem in the application of this technique is the uncertainty of having 
sufficient particle concentration. Engine designers are improving the com- 
bustion process with a goal of eliminating particulate emissions. Character- 
istics of the technique are: 

• Direct Measurement 

• Noncontact 

• Uses infrared emission of suspended particles 

• Particles are small and stay close to the gas stream 
temperature. 
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• May be applied as single color or two-color technique. 

• Straight path averaging characteristics, but limited at high 
particle concentration due to extinction. 

• Hot background metal surfaces present potential problem of 
biasing the measurement. 

• Near-infrared available spectrum allows use of current glass 
or silica fiber optics. 

• Particle concentration may not be sufficient at all power 
settings to obtain an accurate measurement. 


Hot Target Aspirated Ceramic Pyrometer 

This device measures the gas temperature indirectly, by aspirating a 
sample of the gas through a porous target of refractory material, and measur- 
ing the target radiametrically. Figure 7, Hot Target Aspirated Ceramic Pyrom- 
eter, shows the general arrangement. This has been reproduced from Reference 
17. The figure shows a photodiode could be located remotely and coupled by a 
fiber optic cable. The walls of the sensors could be cooled by compressor 
discharge air. However, some energy is used in both the aspiration and cool- 
ing processes. Characteristics of the sensor are: 

• Indirect Measurement 

• Immersed probe with some disturbance on the measured gas stream 

• Uses infrared radiation from a heated target. 

• Use of single color or two-color detection of the target 
radiance is possible. 

• Local sample of gas only 

• Can be ruggedly constructed, but depends on air cooling to 
prevent melting and loss. 

• Aspirated passage and cooling passage may need periodic 
cleaning maintenance. 


Hot Target Unshielded Metal Pyrometer 

This sensor measures the gas temperature indirectly by immersion of the 
sensor in the gas stream. Figure 8 shows the arrangement. The metal jacket, 
closed at the outboard end, is heated by the gas flow and the interior becomes 
a radiating cavity. A sapphire crystal is used as a high temperature light 
pipe to conduct the internal tip radiation to a fiber optic cable in a cooler 
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Figure 7. Hot Target 
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Figure 8. Hot Target Unshielded Metal Pyrometer. 
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region. The radiant energy is detected and processed at the Signal Processor 
and Display. 

The metal jacket could be made of a platinum-family alloy which has 
good impact energy absorption at elevated temperatures. However, these alloys 
have low creep strength and also must be protected against loss due to vola- 
tile oxide formation. Characteristics of the sensors are: 

• Indirect Measurement 

♦ Immersed probe with local (but small) disturbance on gas stream 

♦ Uses infrared radiation from protected cavity in the metal 

jacket. 

• Single-color or two-color detection is possible. 

♦ Local sample of gas only 

• Mechanically simple design, but materials selection and protec- 
tion are critical to success. 

♦ Sapphire bulk optical properties may affect high temperature 
accuracy. 

Sapphire Crystal Thermometer 

This candidate sensor is similar to the hot target unshielded metal sen- 
sor except that the radiating cavity is formed on the tip of the light pipe. 
Figure 9 shows the simplest arrangement possible. Here, no protection or 
aspiration is used. Figure 10 shows the cavity-forming coating and various 
other coatings to protect the cavity and to protect the light pipe proper- 
ties. No shield would be used. 

An alternative configuration would be to place the sensing element in a 
shielded aspirated housing (Figure 11). This housing also provides protection 
against particle impact. This configuration is mechanically complex; thus, 
materials selection presents a problem. For this reason, the first configu- 
ration discussed above (the unshielded crystal) is preferred and was selected 
as one of the five candidates. Its characteristics are: 

• Indirect Measurement 

♦ Immersed probe with very small local disturbance on gas stream 

• Uses infrared radiation from cavity on tip of light pipe. 

♦ Single-color or two-color detection is possible. 

• Local sample of gas only. 
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Figure 9. Sapphire Crystal Thermometer. 





Figure 11. A Shielded Gas Temperature Measuring Probe 
and System. 
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♦ Mechanically simple design, but materials selection and protec- 
tion is critical to success. 

• Sapphire bulk optical properties may affect high temperature 
accuracy. 

Table 5, the First Version Candidate Evaluation Chart, is explained as 
follows: The PARAMETER column (on the left) contains numbers that refer to 
Table 1, Evaluation Criteria-Parameter Value. The next column, IMPORTANCE, is 
a number between zero and 10, and appears on the right column of Table 1. 
Under each candidate are two columns. The left column of the pair is the 
rating number, between zero and four, that shows how well the parameter is 
satisfied by that particular candidate. The right column of the pair is the 
product of the importance number and the rating number. All products for a 
given candidate are summed to give a total score for that candidate. The 
candidate with the highest score, the Sapphire Crystal Thermometer, is 22 
points above the second highest. 

At the Review of April 1984, it was brought out that the Ratings for 
both the Gas Band and Particle Radiation Pyrometers may be too high. The 
wall radiation interference problem in the Gas Band Pyrometer could be more 
fairly rated in the temperature range Parameters 1, 4, and 7 by using ratings 
of 2, 2, and 2 in place of the 3, 4, and 4 ratings in Table 5. 

Likewise, the anticipated lack of sufficient particles in the gas stream 
should be reflected in a reduced rating in the temperature range Parameters 
1, 4, and 7 for the Particle Radiation Pyrometer- This could be more fairly 
rated by use of a 1, 1, 1 rating in place of the 3, 4, 4 rating in Table 5. 

Three new Candidate Evaluation Charts were prepared with the above 
rating changes made. The three charts differ in the importance given to the 
Subidle Temperature Range and the Future Engine Temperature Range. Table 6 
gives zero importance to both these ranges. Table 7 gives zero importance to 
the future engine temperature range, only. Table 8 uses the same importance 
for all temperature ranges as was used in Table 5. 


3.8 SELECTION OF CANDIDATE PYROMETER 

The contractor and customer mutually agreed upon completion of the Task 
I Review to pursue the Sapphire Crystal Thermometer as the selected candidate 
pyrometer. 
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Table 5. 1700^^ Optical Gas Temperature Sensor 

Candidate Evaluation Chart - First Version 


Parameter 

Importance 

Gas Band 

Radiation Pyrometer 

1 

5 

3 

15 

2 

5 

1 

5 

3 

5 

2 

10 

4 

10 

4 

40 

5 

7 

1 

7 

6 

10 

2 

20 

7 

5 

4 

20 

8 

6 

1 

6 

9 

5 

2 

10 

10 

10 

3 

30 

11 

6 

4 

24 

12 

10 

4 

40 

13 

10 

4 

40 

14 

8 

4 

32 

15 

10 

4 

40 

16 

5 

1 

5 

17 

10 

2 

20 

18 

10 

4 

40 

19 

5 

2 

10 

20 

6 

3 

18 

21 

8 

3 

24 

22 

6 

3 

18 

23 

4 

2 

8 

24 

4 

3 

12 

Total Score 

494 


Particle 

Radiation Pyrometer 

3 

15 

0 

0 

2 

10 

4 

40 

1 

7 

2 

20 

4 

20 

1 

6 

2 

10 

3 

30 

4 

24 


Hot Target 
Aspirated Ceramic 


Hot Target Sapphire 

Unshielded Metal Crystal Thermometer 


First Column: Rating (0 to 4) 

Second Column; Product (Importance x Rating) 
Total Score Line; Sum of Products for Candidate 





















Parameter 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Importance 

0 

0 

0 

10 

7 

10 

0 

0 

0 

10 

6 

10 

10 

8 

10 

5 

10 

10 

5 

6 
8 


Gas Band 

Radiation Pyrometer 


2 

1 

2 

2 

1 

2 

2 

1 

2 

3 

4 
4 

4 

4 

4 

1 

2 

4 

2 

3 

3 


0 

0 

0 

20 

7 

20 

0 

0 

0 

30 

24 

40 

40 

32 

40 

5 

20 

40 

10 

18 

24 


Particle 

Radiation Pyrometer 


1 

0 

2 

1 

1 

2 

1 

1 

2 

3 

4 
4 

4 

4 

4 

1 

2 

4 

2 

3 

3 


0 

0 

0 

10 

7 

20 

0 

0 

0 

30 

24 

40 

40 

32 

40 

5 

20 

40 

10 

18 

24 


Hot Target 
Aspirated Ceramic 


3 

3 

4 

4 

3 

4 

3 

3 

4 

4 

3 

3 

4 
3 

3 

1 

1 

4 

3 

2 

3 


0 

0 

0 

40 

21 

40 

0 

0 

0 

40 

18 

30 

40 

24 

30 

5 

10 

40 

15 

12 

24 


Hot Target 
Unshielded Metal 


3 

3 

4 

4 

3 

4 

2 

3 

3 

4 
2 
2 

4 

3 

3 

1 

1 

4 

3 

4 
3 


0 

0 

0 

40 

21 

40 

0 

0 

0 

40 

12 

20 

40 

24 

30 

5 

10 

40 

15 

24 

24 


Sapphire 

Crystal Thermometer 


3 

3 

4 

4 

3 

4 

2 

3 

3 

4 
4 
4 

4 

3 

3 

1 

1 

4 

3 

4 
3 


0 

0 

0 

40 

21 

40 

0 

0 

0 

40 

24 

40 

40 

24 

30 

5 

10 

40 

15 

24 

24 



3 18 
2 8 
3 12 

3 18 
2 8 
2 12 

2 12 
2 8 
2 8 

2 12 

3 12 

2 8 

2 12 
3 12 
2 8 

408 

398 

417 

417 

A49 
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4,0 DESIGN AND FABRICATION 


4,1 SENSOR STRUCTURAL INTEGRITY AND OPTICAL PROPERTIES 


Any indirect, immersed-type of sensing element will have to withstand the 
load of the flowing gas stream, engine vibration, and occasional particle 
impact. In addition, for the sapphire crystal thermometer the light pipe 
property of the crystal must be maintained. 

An in-depth analysis of structural loading stress and the strength of two 
designs was performed and is presented in Appendix E. The cylindrical sensor 
was analyzed first. Problems foreseen in the cylindrical sensor led to the 
invention of a conical shaped sensor, which was also analyzed both mechani- 
cally and optically. The thermal characteristics, found to be applicable to 
the cylindrical and conical sensors, was analyzed and presented in Preliminary 
Design, Analytical Studies. Thermal commonality is reviewed in Appendix E. 


4.1.1 Cylindrical Sensor 

The analysis in Appendix E showed that, for an immersion of 1 inch and a 
diameter of 0.040 inch, the stress at the base would approach 42,840 psi 
tensile. See Figure E-1. A shear stress of 570 psi was calculated. The 
cylindrical sensor depends on the light-pipe effect to provide the trans- 
mission of light from the tip cavity to the fiber optic cable. In turn, the 
light pipe effect works only with a clean or protected surface. 


4.1.2 Conical Sensor 

The structural analysis given in Appendix E was for 1 inch immersion, a 
cone angle of 6.25 degrees, and the same gas stream conditions of Mach 0.2 and 
400 psia. The maximum tensile stress was calculated to be 1016 psi. The 
shear stress was calculated to be 76 psi. See Figure E-2 of Appendix E. 
Optically, the conical sensor does not depend on the light-pipe effect but 
rather, the cone becomes a window for an image of the tip to be focused on the 
fiber optic cable tip. See Figure E-3 and the attendant text of Appendix E. 


4.1.3 Preferred Design 

The conical sensor was preferred over the cylindrical sensor for the 
following reasons: 

♦ Lower stress (40 times lower for the selected case) 

• Freedom from calibration changes due to deposits on the surface 
exposed to the gas stream. 
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4.2 SENSOR HOLDER 


The sensor holder supports the sensing element in the hot gas stream, 
supports the interfacing optic element (s), and the fiber optic cable. The 
assembly is supported by a cylindrical fitting. This is shown in Figure 12. 


4.2.1 Thermal 

In an engine application, cooling air would be available to limit the 
holder temperature. The actual temperature limit would be dependent on the 
engine model and cooling air source. For design purpose, the holder end 
nearest the combustion gas flow will be assumed to be hottest at up to 538® C 
(1,000° F) and the fiber optic termination will be limited to 80° C (176° F) 
so that plastic clad silica fiber may be used. 


4.2.2 Optical, Fiber Optic and Connector 

The theory behind the conical sensor was explained in Appendix E. The 
main elements, as shown in Figure E-3 "Conical Sensor Operation" may be com- 
pared with the cross-sectional illustration of Figure 12. In Figure 12, the 
fiber optic cable, 1, is threaded to attach to end piece, 3, with locking by 
nut, 2. The housing, 8, carries the solid single crystal aluminum oxide cone 
element 7, the aperture-spacer ring 6, the lens 5, and the image aperture - 
spacer ring 4. Housing 8 is attached to the simulated combustor wall 10 by 
compression fitting 9. A compliant sleeve of ceramic paper, 11, is positioned 
between the cone element 7, and the housing, 8. The tip of the cone element, 
12, is the actual sensing zone. This zone is 0.35 inches long and has a 
sputtered refractory metal coating with an aluminum oxide protective overcoat. 
The measured temperature will be the average temperature along the 0.35" 
length of this sensing zone. However, it should be pointed out that sapphire 
exhibits bulk optical absorption properties that may be significant at high 
temperature. This was reported in reference 20. This could have the effect 
of extending the sensing zone. 

Table 9 is a parts-list description and a specification for sensor and 
holder. Where a vendor part is listed, it is understood that an equivalent 
item from an alternative vendor could also be used. Figure 13 is a photo of 
the sensor assembled with uncoated cone next to uninstalled coated cone. A 
detailed mechanical drawing of the sensor assembly is shown in Appendix F. 


4.3 ELECTRONICS 


The electronics block diagram for the 1700° C gas temperature system is 
shown in Figure 14. The fiber optic cable from the sensor assembly terminates 
in an SMA style standard connector at the detector, a commercial two-color 
device, having two separate outputs. Each output is connected to an amplifier 
module. Module outputs are converted from analog to digital and input to a 
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Table 9- Sensor and Holder Parts, Description and Specification. 


Part No. 

Name 

Description/ Specification 

1 

Fiber Optic 
Cable 

200 micrometer core plastic clad silica, SMA 
style connectors, Newport Corp., part FCC-PC-2 
(2 meters), FCC-PC-10 (10 meters), or equivalent. 

2 

Lock Nut 

Part of inline adaptor, Newport part F-CA-001, 
or equivalent. 

3 

End Piece 

Stainless steel 300 series. 

4 

Image Aperture- 
Spacer Ring 

Stainless steel 300 series, 0.370 inch diameter, 
0.630 inch long, aperture 0.007 to 0.010 inch 
diameter to match lens used. 

5 

Lens 

Plano-convex, fused silica, diameter 0.30 inch, 
focal length 0.44 to 0.47 inch. 

6 

Aperture-Spacer 

Ring 

Stainless 300 series, 0.370 inch diameter, 0.125 
inch long, 0.260 inch aperture. 

7 

Cone Element 

Sapphire, low scatter optical grade, 0.300 inch 
base, 2.74 inch height, zero degree orientation 
base flat to one wave vis. yellow, 80-50 polish. 
Cone side to be left fine-ground. Tip truncated 
0.000 to 0.020 inch on height. 

8 

Housing 

Stainless 300 series, 2.72 inches length, 0.60 
inch maximum diameter. 

9 

Compression 

Fitting 

Make from standard commercial 3/8 inch fitting 
stainless 300 series. 

10 

Simulated 
Combustor Wall 

Stainless 300 series, 0.062 inch nominal thick- 
ness, size to fit bench experiment. 

11 

Sleeve 

Ceramic paper, 1/32-inch thickness free-standing, 
Cotronics Corp. No. 300-020 or equivalent cut to 
provide one layer against cone element. 

12 

Tip 

Length 0.35 inch nominal, defined by sputtered 
iridium coating with aluminum oxide overcoat. 
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Figure 13. Photo of Sensor Assembly With Uncoated Cone Next to Uninstalled 
Coated Cone. 
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Figure lA. 1700° C Gas Temperature System - Electronics Block Diagram. 
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conversion program in a computer. The computer calculates the sensor temper- 
ature based on the sensor calibration data and then calculates the gas tem- 
perature for the given heat transfer environment of the sensor. 

The components of this signal processing and computing system are listed 
in Table 10 with description and specifications. Where a specific manufac- 
turer's model number is used, it should be recognized that an equivalent item 
from an alternate vendor could also be used. An electronics diagram in 
detailed form is presented in Appendix G. 


4.3.1 Detector 

Several requirements became evident for the detector system. These are 
listed here: 

• Efficient use of infrared/visible spectrum for adequate signal 
strength at engine idle of 520° C 

• Spectrum within range of available fiber optics such as fused 
silica, 0.3 to 2.0 microns 

• Wide dynamic range of signal processing capability due to the 
temperature span of 520 to 1700° C 

• Capable of single or dual wavelength operation 

• Stable optically, electrically, and mechanically. 

The silicon detector for the above temperature range was chosen as having 
the best stability and other desirable features such as good linearity over a 
seven decade range. Stable operation from -55 to +150° C appears to be possi- 
ble. Two methods are available to achieve dual wavelength (two color or color 
ratio) operation: one uses external filters such as dielectric interference 

types. The second method uses the window property of the silicon slice. This 
second method appears to have the advantage of superior long-term stability. 
The near infrared filters needed for the first method are not stable in high 
temperature and/or high humidity environments. 

The selected detector uses the window filter property (second method) . 
One vendor, EG+G Electro-Optics, has been manufacturing a product called ”Two 
Color Detector** for several years. This unit has tandem silicon photodiodes 
in a hermetically sealed TO-5 package. Each photodiode has separate output 
wires. The top diode has a spectral range of 0.40 to 1,10 microns. This 
diode also acts as a window-filter for the second or bottom diode which has 
effectively a 0.95 to 1.10 micron spectral range. Overall, this detector is 
very efficient - it does not waste radiant energy in the 0.4 to 1.1 micron 
range . 

The long wavelength cutoff is environmentally temperature dependent, as 
is true for all detectors. This amounts to about a +1% signal output change 
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Table 10. Electronics for 1700® C Temperature Sensor. 


Name 

Fiber Optic 
Cable 

Two-Color 

Detector 

Temperature 

Controller 

Power Supply 
Amplifier 
Multiplexer 
A to D Converter 
Computer 

Disk Drive 
Monitor/Display 


Description/Specification 
See Table 9. 


Tandem silicon photodiodes, TO-5 package, 
parallel outputs EG+G, Electro-Optics Div. , 
or equivalent. 

Feedback loop based on AD590 Temperature 
Transducer Analog Devices Corp. or equivalent. 

±15 volt d.c. supply, Acopian DB15-20 or 
equivalent. 

Input range IpA to 1mA, Analog Devices AD515LH 
or equivalent. 

16-channel analog input, part of DIADAC 1, 
Microtech Co. or equivalent. 

12-bit, memory mapped, part of DIADAC 1, 
Microtech Co. or equivalent. 

BASIC in ROM, keyboard, accessories bus, 64K 
RAM, RS232-type output. Commodore C64 or 
equivalent. 

Single drive, 5-1/4-inch flexible disk 
Commodore 1541 or equivalent. 

CRT based, audio speaker. Commodore 1702 or 
equivalent. 



per degree C (diode case temperature) at 1.04 microns. Roughly, a 1° C change 
in diode case temperature would appear as a 1° C change in target (sensor tip) 
temperature when the target is at 1100° C. Thus, it is apparent that the 
diode case temperature must be either controlled or sensed and accounted for 
in the calculation of sensor tip temperature. Temperature control has been 
selected for this design. 


4.3.2 Amplifiers 

The dynamic signal range is expected to be approximately 250,000 to 1 for 
the target temperature range of 520 to 1700° C. Even a l6-bit A to D convert- 
er cannot handle this dynamic range, being limited to one part resolution in 
about 65,000 parts. Two options are available: range changing and range 
compression. Range changing makes use of either switched amplifiers of 
various gains or switched gain on one amplifier. Range compression makes use 
of a nonlinear, but stable, device such as a diode- or transistor-based loga- 
rithmic amplifier. Range compression has the advantage of no discontinuities 
in time or calibration data since no switching is used. But a logarithmic 
amplifier with an input current capability as low as l.OE-11 amperes was not 
available. Therefore, range changing was used by the technique of switching 
the gain on a state-of-the-art low input current amplifier. Supporting 
calculations of photodiode current versus temperature appears in Appendix H. 
The amplifiers were placed on the same temperature controlled-plate that holds 
the detector. Gain ranges are completely computer controlled or they may be 
controlled manually. Details appear in Appendix G. 


4.3.3 Analog to Digital Converter 

A 12-bit analog to digital converter was used to convert the log ampli- 
fier outputs to digital form. A resolution of one-third of a degree C was 
expected over the 520 to 1700° C range. The A to D converter card also 
included a 16-channel multiplexer for analog inputs, a 12-bit digital to 
analog output converter, a battery-backed clock calendar, and 12-digital 
input/output lines. 


4.3.4 Computer, Display, and Printer 

A personal computer, single-disk drive and CRT monitor, has been select- 
ed as the most efficient and versatile method to calibrate the gas temperature 
sensor and to provide the conversion from volts to degrees C. The computer 
also has the capacity to calculate the actual gas temperature from the sensor 
temperature when the sensors heat transfer environment is known. Note that 
the gas temperature will be equal to the sensor temperature only if the 
surrounding wall temperature is close to or equal to the gas temperature. 
This could be the case for an uncooled turbine, but it is not true for an 
engine area having cooled metal parts. Refer to the section on Analytical 
Studies and Figures 2 through 5 for a discussion and example calculations. A 
dot matrix printer was added in order to make permanent records of programs 
and data. 
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4. 3,5 Computer Programs 


Four programs were used in the development of the temperature sensor. 
These will be briefly described here. Complete listings of the programs and 
remarks corresponding to program line numbers are included in Appendix D. All 
programs will run as-is on the C64 computer. 

RESP-96D - This program is used to enter component calibration data from 
detectors, fiber optics, and spectrometer wavelength dial readings. The pro- 
gram also calculates spectral transmission of the individual fiber optics 
and detectors under test. The program was written by Paul M. Clark under this 
contract. 

WAT-CA96D - This program takes the calculated spectral transmission data 
generated in RESP-96D and combines this data to yield system calibration 
curves. Since the spectral data is in small increments (0.01 micrometers), 
Planck's equation can be used to calculate system output versus target temper- 
ature. This is explained in the CALIBRATION section of this report. This 
program was written under contract by Paul M. Clark, 

A/D READ MIO - This program operates the input multiplexer, the analog 
to digital converter, and looks up an indicated temperature using the system 
calibration curves. The program also looks up the time and date correspond- 
ing to the data acquisition cycle; it prints to either display monitor screen 
or to the dot matrix printer. This program is the only one needed to operate 
a probe that is already characterized. The program was written under contract 
by Paul M. Clark. 

PM-TEMP4A - This program takes the indicated probe temperature and calcu- 
lates the gas temperature based on input of the probe heat transfer environ- 
ment. The program has the capability to be used in reverse. If the gas 
temperature and heat transfer environment are put into the program, the probe- 
indicated temperature will be calculated. This program was adapted for the 
C64 computer by Paul M. Clark from the thermocouple heat transfer correction 
program written (not under this contract) by Robert C. Williamson of General 
Electric Company. See Reference 11, Williamson and Stanforth, for a complete 
description of this work. It should be noted that this program was written 
for corrections to cylinder-shaped sensors in either pure perpendicular or 
pure cross flow. The application to a conical-shaped sensor was done by 
approximating the conical tip to a cylinder. A better approach would be to 
calculate the heat transfer environment to a cone using finite element 
methods, but this approach was beyond the scope of this contract. 
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5.0 CALIBRATION 


Two approaches could be used to calibrate the sensor. Since the sensor 
uses radiant energy from a cavity to measure the temperature of the cavity, 
Planck's equation, which describes the radiant power at a given wavelength, 
could be used. This is the first approach. An alternate approach would be 
to calibrate the sensor against a reference sensor over the range. This 
alternate approach allows no extension of the calibration outside of the 
temperature range of the reference sensor. 

The first approach was chosen because it not only allows extension 
(extrapolation) outside of the reference sensor temperature range, but it pro- 
vides a self-check within the reference sensor temperature range. This will 
be described in detail. For the method to work correctly, Planck's equation 
must be applied over the entire spectral range of response - wavelength by 
wavelength. 

The calculated detector output is found from the sum of every spectral 
wavelength contribution to the photocurrent as follows: 


f ^2 

Itot ~ ^ J Total Photocurrent Equation 




(.» 


Planck's Equation 


where 


Itot detector output photocurrent, amperes 
K is the system overall efficiency constant (dimensionless) 
is the detector spectral response function (dimensionless) 
is the fiber optic spectral response (dimensionless) 
is Planck's equation for spectral radiation: Watts/CM^/CM 

A is the wavelength in centimeters (CM) 

T = temperature, Kelvin degrees 
Cl = 3.7413 X watt CM^ 

C 2 = 1.4388 CM degree 
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5 . 1 SPECTRAL RESPONSE CALIBRATIONS 


Although manufacturers usually give "typical" spectral response curves, 
these have been found to be not sufficiently accurate for use in calibration 
of a precision sensor. Therefore, the spectral response was calibrated using 
the system shown in Figure 15. A tungsten strip filament lamp, GE 18AT10, 
was operated on a regulated direct-current supply. A condensing lens illumi- 
nated the entrance slit of a Bausch & Lomb Model 33-86-25-03 Grating Mono- 
chromator. At the exit slit is a Kodak No. 89B filter to block overlapping 
orders from the grating. This equipment was calibrated in terms of radiant 
power versus wavelength, using an EG&G detector SGD-IOOA Serial No. 01-03-008 
which has an NBS traceable spectral calibration. Once calibrated, the spec- 
tral response of the holder, fiber optic, detector, and amplifier (as a 
system) was measured. The detector/ amplifier system was then separately 
characterized. Both a 5-meter and a 10-meter length fiber-optic cable were 
characterized. 

Computer program RESP-96D, described under the Design and Fabrication 
section and in the detailed listing in Appendix D, was written and used to 
enter the spectral data from the components in computer file form. 


5.2 PRECISION RADIOMETRIC CALIBRATIONS 


Next, the system overall efficiency constant, K, was found. This is 
illustrated in Figure 16. Figure 17 is a photo of the equipment set up to 
perform the precision radiometric calibration. The holder (without the crys- 
tal in place), fiber optic cable, and detector/ amplifier were set up to view a 
blackbody radiation standard. The standard used was the Electro-Optical 
Industries Model WS154, having a 1-inch-diaraeter port and an 11-degree field 
of view. A type "S" thermocouple system, Kaye Instruments Ice Point Refer- 
ence, and Mediator A-75A Precision Millivolt Potentiometer were used to set 
the radiation standard to the NBS traceable temperature levels. 

When the constant K from the above equation has been found at one known 
temperature, the radiation standard was reset to a different temperature and 
the 1700° C system output indicated the temperature was compared. This 
calibration physically does the total spectral integration called for in the 
Total Photocurrent Equation above. The blackbody radiation standard effec- 
tively substitutes for the crystal, thereby providing the most accurate source 
of radiant temperature outside of melting point standards. This testing has 
given information on the tracking accuracy of the computed temperature cali- 
bration curve versus the actual data. The two methods agree to about 5° C. 

The testing has also shown a repeatability of 2° to 3° C at 1000° C. 
Some temperature sensitivity on the optical fiber was discovered. This is 
correctable by proper design of optical stops and is discussed in the Results 
and Evaluation section. 
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Figure 15. Spectral Response Calibration. 




Ice Point Reference 



Figure 16. Precision Radiometric Calibration. 
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6.0 TESTING 


Tests were performed on components and on the system for the 1700® C gas 
temperature sensor. Several types of tests were considered necessary to fully 
characterize and evaluate the sensor. Two furnaces were constructed and a 
device to evaluate the coating integrity was built. These items will be 
described below. 


6.1 FLAME-POWERED CAVITY FURNACE 


An insulated furnace was built for the purpose of heating the sensor 
gradually and evenly. Figure 18 shows the setup used for testing. A propane- 
air burner was used for heating. A disappearing filament optical pyrometer 
was used to obtain a reference tip temperature. The Micro-Optical model by 
Pyrometer Instruments Company is designed to obtain radiant temperatures from 
very small targets, such as the tip of the cone sensor (0.005 to 0.015 inch 
diameter). This pyrometer has an NBS traceable calibration. The furnace 
operates by free convection with an internal velocity of 15 to 20 feet per 
second. The maximum temperature that the cone could be heated was 1350® C, 
even though the adiabatic flame temperature is 1920® C. This is because the 
low velocity results in a low coefficient of convective heat transfer. A 
second furnace was built for the purpose of heating the sensor in the range of 
1350® to +1700® C. This will be described next. 


6.2 VACUUM ASPIRATED FURNACE 


This furnace was built to achieve an improved convective heat transfer 
coefficient on the sensing element. Figure 19 shows the general arrangement. 
The sensor tip was placed in a ceramic nozzle attached to a steel chamber. 
The chamber was attached to a vacuum pump. A propane-air torch was located to 
fill the flow entering the nozzle. The reference optical pyrometer was locat- 
ed to sight on the tip. A cool air source kept the fiber optic cable at a 
safe temperature. A differential pressure gage was used to monitor and 
control the vacuum applied and, thus, the flow velocity in the nozzle. The 
furnace was operated at about 220 feet per second; testing with a sapphire 
cone element achieved 1750® C. It was observed that the sensor tip tempera- 
ture was very stable with time. Time response testing was done by chopping 
the torch flame with an aluminum plate "flag.” 


6.3 SENSOR AND COATING INTEGRITY TESTS 


The purpose of the coating on the sapphire crystal cone element is to 
provide the blackbody radiation cavity. When this coating is applied cor- 
rectly, it will have a very low light transmission. That is, it will be 
opaque. One can simply hold the cone against a uniform light source and check 
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Figure 18. Flame-Heated Cavity Furnace Test Arrangement. 
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Figure 19. Vacuum Aspirated Furnace. 





for transmission. However, it was desired to quantify this optical property, 
and a small fixture was designed. This is shown in Figure 20. A light source 
of high-intensity illumination projects light at an angle to the tip. A black 
paper mask prevents light from entering the cone in the uncoated region. A 
calibrated photodetector receives light from the base of the cone. An uncoat- 
ed cone of the same dimensions was used to provide a benchmark. It was 
necessairy to supply regulated power to the lamp to avoid false changes due to 
line voltage fluctuations. 

Other coating integrity tests or evaluations used successfully were 
examined by optical stereo microscope and by scanning electron microscope. 
Direct weighing was attempted but was unsuccessful for the following reason: 
the weight of the film was estimated to be 50 to 100 micrograms. For the 3.8 
gram cone, the electronic weight scale used exhibited a least count corre- 
sponding to 100 micrograms. Therefore, it became apparent that film weight 
loss or gain would be undetectable. 


6.4 SUMMARY OF TESTS 


Four sapphire crystal cone elements were procured. Each element was 
nearly identical in length, cone angle, and weight. All had an optically 
polished base and a fine-ground cone side. Two elements were given a nominal 
0.2 micron- thick platinum, 30% rhodium coating with a 1.5 micron aluminum 
oxide overcoat. The third element was given a 0.35 micron- thick platinum, 6% 
rhodium coating with no overcoat. The fourth element was left as-is for 
control puirposes and was not coated. The tests are summarized in Table 11 
and will be described here in chronological order. 

Referring to Table 11, the cone elements Number 1 and 2 were given coat- 
ings as specified in Figure 21. The subcontracted vendor used continuous 
rotation in the sputtering chamber. The coated cones had a very uniform 
appearance, but were not optically opaque. 

The as-received coated cone elements were tested for light transmission 
by comparing them to the light transmitted by an uncoated cone. One of these 
coated cones. Number 2 in Table 11, was set aside for use as a control. 

The other coated cone element. Number 1 in Table 11, was given a progres- 
sive elevated temperature heat treatment. Gradual heating was done first in 
an argon atmosphere to 1000® C, then in an air atmosphere, also to 1000° C. 
The cone was weighed before and after heat treatment. The change was not 
significant (see discussion above). However, a light transmission measurement 
showed increased transparency after the heating. Under the microscope, no 
cracks, chips, or spalling could be found. The coating appeared to be 
adherent but thinner optically. 

The system was tested with the heat treated, coated cone element No. 1 
installed in the probe holder. A propane-air burner was used to provide a 
variable temperature source. A muffle furnace was used to create a uniform 
cavity so that a reference optical pyrometer could be used. Testing was done 
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Figure 20. Coating Integrity Test Configuration. 
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Table 11. Test Summary Chart. 


Cone 

Element 

Coating/Overcoat 

Applied 

C-C 

Oven 

Argon at 1 ATM 

C-C 

Oven 

Air at 1 ATM 

C-C 

Propane Air Flame 

C-C 

1 

0.2 Microns 

G 

RT to 760° C and 

G 

RT to 871° C and 

F 

Heated in Steps to 

P 


PT-30% RH/ 


Hold 1 Hour 


Hold 1 Hour 


1350° 



1.5 Microns 


RT to 871° C and 


RT to 982° C and 





AI 2 O 3 


Hold 1 Hour 

F 

Hold 1 Hour 

G 






RT to 982° C and 







*Continuous Rotation 

Hold 1 Hour 







in Coating Chamber 








2 

Same as No. 1 

G 

None 


None 


None 





Set Aside for use 

as a Control Sample) 




3 

0-35 Microns 

E 

None 


RT to 260° C and 

E 

Heated in One Step 

G 


PT- 6 % RH/ 




Hold 1 Hour 


to 1750° C 



No Overcoat 




RT to 538° C and 


Heated RT to 1700° C 







Hold 1 Hour 

E 

and Return to RT 

F 


*Two Positions 




RT to 816 ° C and 


Total of 107 Cycles 



in Coating 




Hold 1 Hour 





Chamber, July 




RT to 1000° C and 

E 




1985 




Hold 1 Hour 

E 



4 

No Coatings 


None 


None 


Heated to 1700° C 

G 








and Return to RT 









Minimum of 25 Cycles 


C-C: Condition-Code 






E = Excellent 






G = Good 






F = Fair 






P = Poor 










Obj ective : Apply uniform and adherent coatings of Pt-30Rh (2000 A) 

^ 20 % thickness and overcoated with a thick coating of 
AI 2 O 3 ( 1.5 microns minimum thickness to provide protection 
of Pt-30Rh film from oxidation). 

Notes: 1. Masking must be provided to keep coatings confined to 

tip side of cones as shown. 

2. Vendor to be responsible for cleaning process - 

3. Samples should be rotated or, at a minimum, be 
manually turned at least once for each coating. 


Figure 21. Sputter Coating of Two Sapphire Cones. 
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in steps up to about 1350^ C. The platinum, 30% rhodium coating was examined 
visually after each temperature step (allowing the probe to cool first) . Some 
increased transparency was noted after each step. The most significant change 
was seen after the immersion at 1350°C. It is believed that the coating 
oxidized and evaporated. The original specification called for a thickness of 
0.2 microns. Based on the results so far, this may be too thin for the envi- 
ronment. There was some suspicion that the coating may have reversed sput- 
tered, resulting in a thickness of less than the specified 0.2 microns. This 
thickness could not be measured without destroying the part. 

One of the sapphire crystal cone elements, Number 3 in Table 11, was 
successfully coated in a new in-house facility. A platinum, 6% rhodium alloy 
was used. No alumina overcoat was used, so that the adhesion of the metallic 
coating could be evaluated. Some pinholes and scratches were evident, but 
are not in a location that would cause problems with the operation of the 
sensor. Coating was done with one manual rotation. Visually, a minor cir- 
cular nonuniformity could be seen, but only under high intensity illumination 
internally. 

This coated sensor, Number 3 described above, was given an elevated 
temperature soaking in air at four temperatures up to and including 1000° C. 
The holding period was one hour at each temperature. Preliminary inspection 
showed no distress from the heating. 

The Number 3 coated cone sensor was next tested in the new vacuum aspi- 
rated furnace. The sensor indicated that a temperature of about 1750° C was 
reached. The optical pyrometer used as a reference indicated the tip tempera- 
ture was within about 25° C of this value. After this, 107 cycles from room 
temperature to 1700° C were run. The platinum, 6% rhodium coating was 
examined under the stereo microscope. The coating appeared to have become 
thinner and darker. 

Next, the sensor was examined with the aid of a scanning electron micro- 
scope. The pictures confirm that the coating is adherent and free of chips, 
spalls, and cracks. It was desired to obtain thickness information from the 
SEM pictures, but this could not be accomplished without cutting and polishing 
a section of the sapphire cone. This was considered undesirable because the 
cone sensor would be unusable if cut. The SEM pictures are reproduced in 
Appendix J. 

Time-response tests were also run using the aluminum plate flag and the 
vacuum aspirated furnace. See the Results and Evaluation section. 



7.0 RESULTS AND EVALUATION 


Tests were designed to evaluate the parameters as listed in Table 12. 
Under "Accuracy" (Parameter 3) is included linearity, repeatability, and 
hysteresis. This table includes results in two columns: Evaluation Rating 
and Best Estimate and Comments. The evaluation rating is explained in the 
Preliminary Design section of this report. 


Calculation of Gas Temperature 

An example calculation of gas temperature was done by using sensor (indi- 
cated) temperature from a particular reading taken during the tests in the 
vacuum aspirated furnace. Table 13 is a compilation of sensor mechanical data 
and known heat transfer environment data. From this set, program PM-TEMP4A 
was used to calculate the gas temperature for various combinations of equiva- 
lent sensor diameter and sensor emittance. Table 14 shows the results. In 
this example calculation, it can be seen that the gas temperature as calculat- 
ed shows 80® C difference from extreme values taken for sensor external 
emittance and sensor equivalent diameter. In an actual engine, the mass flow 
parameter would be much higher than the 5 Ibm/sec sq ft used here, and the 
range of calculated gas temperature would likely narrow to less than 10® C. 
This is evident in Figure 5, appearing earlier in this report, where the mass 
flow parameter reached a value of over 160 inside the engine at full throttle 
calculation. 



Table 12. Test and Evaluation Results. 



Parameter 

Range or Goal 

Evaluation 

Rating 

(See Table 2) 

Best Estimate 
and Comments 


1 . Temperature 

170° to 1700° C 
Emphasizing the 
600° to 1700° C 
Range 

3 

510 to 1750° C 

— 

2. Resolution 

± 1° C 

4 

± 0.1° C 

— 

3. Accuracy 

4. Environment 

± 2° C 

2 

± 25° C 


a . Media 

Combustion Products 
of Hydrocarbon Fuel 
and Air 

3 

Some Loss of 
Coating Observed 

— 

b. Velocity, 
Pressure 

Mach 0.2, 25 atm 

3 

Gas Loading 
Calculations 

— . 

5 . Time Response 

3 Seconds for 63% 
of Step Change 

4 

Less than 0.25 
Seconds 

r 

6 • Endurance 

100 Cycles, Room 
Temperature to 
1700® C 

2 

Some Loss of 
Coating Observed 

— 

7- Installed 

1000 Hours or More 

1 

Additional Devel- 
opment of Coating 
Needed 

— 

8 . Induced 

Engine Losses 

0.05% Maximum 

4 

Very Low Aerody- 
namic Blockage, 
Area Ratio 0.0003 

. 

9, Risk to Engine 
(Foreign 
Object Damage) 

Estimate from 
Strength of 
Materials Used 

3 

Low Risk Predicted 
from the Stress 
Analysis 


10. Vibration 
Endurance 

Estimate Only 

3 

Thermal 


11. Fiber Optic 
Coupling 

6 Meters 

4 

5- and 10-Meter 
Lengths Used 
and Compared, a 

6- Meter Length 
Would Have No 
Problem With 
Attenuation 
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Table 13. Constant Input Data. 


Sensor Equivalent Diameter 

Variable (See Table 14) 

Sensor Equivalent Length 

0.25 inch 

Sensor Emittance (External) 

Variable (See Table 14) 

Flame Angle Factor 

0.05 

Flame Emittance 

1.0 

Mass Flow 

5.0 Ibm/sec sq ft 

Wall Temperature 

3281® R (1549® C) 

Sensor (Indicated) Temperature 

3581® R (1716® C) 


Table 14. Gas Temperature Calculations. 


Sensor 

Emittance, 

External 

Equivalent 
Diameter, inches 

Calculated Gas Temperature 

° R 

° C 

0.4 

0.016 

3721 

1794 

0.5 

0.016 

3754 

1812 

0.6 

0.016 

3786 

1830 

0.4 

0.032 

3780 

1827 

0.5 

0.032 

3824 

1851 

0.6 

0.032 

3865 

1874 




8.0 DISCUSSION 


The temperature range (see Table 12) found on this first sensor was 510° 
to 1750° C. The lower limit was determined by the decreasing signal strength; 
the upper limit is due to the melting point of the platinum, 6% rhodium coat- 
ing, approximately 1830° C. This range will be well suited for gas turbine 
engine combustor discharge conditions at idle throttle and above. 

The resolution result of 1° C is a good feature and stems from the low 
system noise and the use of the computer-controlled multiple-ranging ampli- 
fier and a 12-bit analog to digital converter. This is a real advantage over 
thermoelectric sensors where the very low output voltage makes the amplifier 
and system noise resolution-limiting. 

The accuracy estimate of ±25° C represents the difference between the 
reference pyrometer and the sensor taken at 90% confidence limit of error. 
This number is made up of several sources of error. This includes amplifier 
gain changes, noise, optical fiber temperature and stress effects, calibration 
coarseness, emissive coating physical changes, and sapphire crystal optical 
optical absorption effects. There are likely other effects present. The 
optical fiber temperature effect found was related to the particular fiber 
specification and the sensor optic design. It was discovered that the fiber 
transmission changed with fiber and end-fitting temperatures. It is theorized 
that the numerical aperture (NA) of the fiber was smaller than the NA of the 
sensor holder coupling lens and the two-color detector. Therefore, a change 
in the NA of the fiber due to temperature (or stress) would come right through 
the system. However, due to the characteristics of the two-color detector, 
most if not all of the transmission change is "cancelled" with regard to the 
indicated temperature. The residual effect seen was about 1° C of indicated 
target temperature change per 1° C of fiber optic temperature change. This 
problem has a solution which is described in the Recommendations section of 
this report. 

Another effect on accuracy could come from the bulk optical absorption 
properties of sapphire (mentioned in the Design and Fabrication section of 
this report). The effect could be to extend the sensing zone beyond the 
emissive coated region of the sensing element. This possibility was not 
tested or analyzed, but it may have a significance, depending on the tempera- 
ture gradient on the sensor. 

The environment parameter consists of the media that the sensor is 
immersed in (combustion products) and the velocity/pressure effects that lead 
to physical stress. Generally, some loss of the emissive coating was observed 
in proportion to the time the sensor was immersed in the flame. The surviv- 
ability of the sensor in a dense, high-velocity gas stream was estimated from 
gas loading calculations. 

The time response* was found by using a strip-chart recorder with a servo- 
pen. A 63% step change was seen as 0.25 seconds; this is also the nominal 
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response time of the recorder. It must be concluded that the sensor response 
time is less than 0.25 seconds. 

The endurance and installed life goals fell short primarily due to the 
loss of emissive coating. There are several possibilities to improve coating 
durability; these are discussed in the Recommendations section. 

The induced engine losses goal of 0.05% was met based on the calculated 
blockage area ratio of 0.0003. 

The risk to engine (foreign object damage) of this sensor (due to break- 
age) is very low because of the low stress of the cone design. For the same 
reason the vibration endurance would appear to be good. 

The fiber optic coupling requirement of 6 meters is quite acceptable, 
based on the low losses seen on both the 5- and the 10-meter fiber optic 
lengths . 
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9.0 RECOMMENDATIONS 


The gas temperature measurement technique pursued in this contract uses 
infrared pyrometry techniques to sense the temperature of a small probe. The 
probe must be in good thermal contact with the gas stream for accurate meas- 
urement. This measurement technique was selected at the end of the prelimi- 
nary design phase of the contract as the best choice. At this time of 
completion of the contract, there is no evidence that would change the evalua- 
tion criteria, or the rating that resulted, such that another technique would 
be favored. Nor has any novel technique emerged that would compete signifi- 
cantly with the chosen system. It is recommended, therefore, to not change 
from this measurement system approach. 

Within this measurement system approach, however, changes or improvements 
are recommended to two of the three components: the sensor and the fiber 
optic. An option of increasing the low-end temperature measurements range is 
suggested for the third component, the detector package. These ideas will be 
discussed. 


9.1 SENSOR 

A. For a solid sapphire sensing element, the effect of bulk optical 
absorption on the sensing zone length should be investigated. 

B. A more exact heat transfer analysis should be done for a conical- 
shaped sensor to improve the calculation of gas temperature. 

C. The durability of the sensor at elevated temperatures should be 
improved by one of the following approaches: 

♦ Use a thicker metallic emissive coating on the sapphire 
cone. Investigate protective coating(s) for the metallic 
coating. 

♦ Use an emissive material that is diffused into the sap- 
phire instead of a coating. 

♦ Use a hollow refractory material cone instead of the solid 
cone of sapphire. 

The third idea will be discussed in some detail. This concept has the 
potential of greatly extending the upper temperature limit of the sensor, 
while at the same time being more durable. The hollow refractory material 
cone for use as a sensing element is shown in Figure 22. The air or gas 
within the cone conducts the light rays (one is shown) that are emitted from 
inside the tip region of the cone. Just as in the solid sapphire cone, a lens 
collects and focuses the rays onto the optical fiber face. If the cone is not 
made of opaque material, then an emissive (and opaque) coating would be put on 
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Figure 22. Hollow Refractory Material Cone for Use as a 
Sensing Element. 



the inside of the tip region. This coating, if used, would be well protected 
since it is on the inside. 

A great extension in both the upper temperature limit and in durability 
would come about by making this hollow cone of a refractory ceramic such as 
yttria-stablized zirconia. With a melting point of 2590° C, the material 
could be usable to possibly as high as 2450° C. This would amount to an 
extension of 500° C over the usable limit of 1950° C for a sapphire solid 
cone. 


9.2 FIBER OPTIC 


The fiber optic energy coupling should be improved by design so as to 
reduce or eliminate temperature effects on the energy transmission. From 
analysis of the current design, this improvement could happen by redesign of 
the optical stops, particularly the final stop at the detector. The goal 
would be to make the final stop numerical aperture slightly less than the 
optical fiber's lowest value of numerical aperture. 


9.3 DETECTOR PACKAGE 


An option of extending the lower temperature limit is available. This 
limit could be extended downward by substituting a longer wavelength detector. 
The spectral response of the detector must still be compatible with the opti- 
cal fiber spectral transmission. 

As an example, an indium-gallium-arsenide detector was chosen for calcu- 
lation purposes. This detector is beginning to be used with certain fiber 
optics. It has a long-wavelength cutoff wavelength of 1.7 microns compared to 
the silicon longwave cutoff of 1.0 microns. Using this optical detector would 
extend the lower temperature limit downward to about 260° C. The current 
system with the silicon detector operates down to 510° C (see Table 12). 
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APPENDIX A - BRAINSTORM SESSION FOR NEW SENSOR IDEAS 


A brainstorm session was held on October 14, 1983 to look at all possible 
approaches to measuring the combustor discharge gas temperature in current 
aircraft gas turbine engines. Present at, and contributing to the list of 
ideas and concepts at this meeting, were Paul Clark, Paul Mossey, Ramayya 
Mulukutla, Wayne Shaffernocker , Chuck Stanforth, Bruce Stowe, and Bill 
Stowell. 

The list of ideas was sorted after the session (without dropping any 
ideas) into the following grouping: 

1. Basic Method: 

A. Direct (Thermodynamic property to gas or trace particle) 

B. Indirect (Immersed probe use) 

2. Modification, Addition, Improvement, or Materials 


Basic Method^ Direct 

• Direct radiation of gas stream 

• Carbon particle temperature 

• Use CO 2 at a narrow spectral region of maximum optical absorp- 
tion and shield pocket to obtain special resolution 

• Use the temperature dependence of the infrared emission 
spectral line width of CO 2 or H 2 O 

• Equilibria of NO 2 concentration (reference Jim Few or McGregor 
at AEDC) 

• IR line reversal (like Na "D” line) . 


Basic Method, Indirect 

• Sapphire fiber blackbody concept as published by R.R. Dils 

• Porous plug hot target - laser drilled 

• Temperature-dependent fluorescent effect in a crystal 

• Quarter wave stunt on end of sapphire 

• Solid state Fabry-Perrot (Rockwell Company) 

• Index of refraction versus temperature 
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• Optical transit in fiber or crystal 

• Panametrics, Inc., work on ultrasonic temperature sensor 

• Johnson Noise probe made with Pt Rh alloys or SiC or SiN^ 

• Bimetal effect, but use ceramics 

• Thermal expansion effect measured by optical interferometry 

(Fenwick) 

• Pulse probe with optical temperature sensor 

• Fiber doped with sodium, and use line reversal. 

Modification, Addition, Improvement, or Materials 

• Single crystal of Si, Ge, or for long IR operation 

• Time and temperature stability of doped target? 

• Need shield material for sapphire fiber thermometer 

- Thoriated dispersed platinum 

- Sapphire fiber strengthened (dispersed in) platinum 

- Sapphire tube 

- Silicone carbide tube 

Q 

- Silicomp or silicon nitride 

• Doped sapphire to make end opaque (as a cavity) 

• Contaminate the sputtering chamber to put opaque coating on 
fiber 

• Sapphire fiber C-axis oriented, then cut the end at an angle 

• Grown special end on sapphire fiber 

• Do fiber endurance test early 

• Flexible sapphire 

• A.D. Little’s work (bent crystal during pulling) 

• Polarization effects of infrared 

• Long wave infrared fibers that are flexible 
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Convert polycrystalline refractory material to single crystal 
for use as light pipe 

Jointed mirrors in tubes for use as (long wave) light pipe 
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APPENDIX B - KEY WORD LIST - OPTICAL GAS TEMPERATURE SENSOR 


Gas Thermometer 

Refractory Optical Material 

High Temperature Optical Sensor 

High Temperature Optics 

High Temperature Fiber Optics 

High Temperature Optical Coating 

Optical Temperature Measurement 

Optical Temperature Sensor 

Gas Temperature Sensor (Measurement) 

Optical Sensor (s) 

Spectroscopic Temperature Measurement 
Passive Optical Sensor 
Electromagnetic Interference (EMI) 
Fiber Optic Sensor(s) 

Infrared Radiation Thermometer 
Gas Radiation Temperature 
Sapphire Crystal Thermometer 
Hot Target Sensor 
Turbine Inlet Temperature 



APPENDIX C - ANNOTATED BIBLIOGRAPHY 


Gas Temperature Measurement 

1. Beretta, F., Cavliere, A., and Dalessio, A., "Visible and Ultraviolet 
Spectral Emission and Extenction Measurements in Oil Spray Flame," Com- 
bustion Science, April, 1980, p 1. Flame emittance and temperature pro- 
file were obtained with optical techniques. Spectral regions included 
the OH, CH, and continuum in absorption and emission from 0.28 to 0.8 
microns wavelength. A swirling burner with Number 2 fuel oil was used. 

2. Bundy, F.P., and Strong, H.M. , "Measurement of Flame Temperature Pressure 
and Velocity," Physical Measurement in Gas Dynamics and Combustion , 
Volume IX of the High Speed Aerodynamics and Jet Propulsion series, pp 
343-388, Princeton University Press, New Jersey, 1954. The use of the 
sodium "D" spectroscopic lines for temperature measurement, as well as 
pressure measurement is given. References important works relating to 
noncontract optical techniques. 

3. Coleman, H.W. , Hardesty, D.R. , and Cattolica, R. J. , "Diagnostics Assess- 
ment for Advanced Power Systems: Interim Report," Sandia Report SAND-77- 

8216, March 1977. Discussion and comparison of optical and nonoptical 
temperature measurements in hot gasses. Emphasis on Raman and CARS. 
Application towards pulverized coal and heavy residual fuel fired gas 
turbines. 

4. Daily, J.W., and Kruger, C.H., "Effects of Cold Boundary Layers on Spec- 
troscopic Temperature Measurements in Combustion Gas Flows," AIAA Paper 
76-134, January 1976. This is an-in-depth study of the use of spectral 
line reversal in gas temperature gradients. A valuable paper in this 
field. 

5. Davis, L.B. , Farrell, R.A. , and Hill, W.E., "Total Flame Radiation and 
Two Color Flame Temperature in Industrial Gas Turbine Combustion Sys- 
tems," ASME Paper 76-HT-57, August 1976. The method used assumes the 
flame to be a "gray body" in the radiation spectrum due to carbon par- 
ticles present in the flame. Commercial pyrometer used. 

6. Dils, R.R., "High-Temperature Optical Fiber Thermometer," Journal of 
Applied Physics, Volume 54 (3), March 1983, p 1198. This paper describes 
a thermometer with sensor made from a single crystal sapphire fiber, 
approximately one mm in diameter. A blackbody cavity is sputtered on the 
end of the fiber. A flexible glass fiber conducts the cavity radiation 
to a detector. A temperature range of 600° to 2000° C with precision of 
better than 1° C is claimed. 

7. Engel, M.E. , "An Investigation of the Thermal Characteristics of the Hot 
Target Sensor for High Temperature Gas Measurements," Master’s Thesis, 
Department of Aerospace Engineering, University of Cincinnati, 1972. A 
detailed analysis of the heat transfer to and from the sensing element. 
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8. Entine, G. , "Turbine Inlet Temperature Sensors," Final Report, Contract 
N00140-74-C-0583, Naval Air Propulsion Test Center, Trenton, New Jersey, 
08628, August, 1974, by Tyco Laboratories, Incorporated, Waltham, 
Massachusetts, 02154. A sapphire light pipe, immersed in the hot gas, 
becomes incandescent and emits light. The light is transmitted through 
the sapphire pipe itself to the photosensor. The claim is made that the 
emission comes not from the surface of the probe tip but from the crystal 
bulk, making surface contamination have "little effect." In the general 
case, this is not true since light pipe operation is impeded by surface 
contamination. The claimed "little effect" could only happen under a 
special case of uniform temperature along the length of crystal subjected 
to surface contamination. Even with this flaw, the report deserves a 
special mention because it was apparently a pioneering work that, unfor- 
tunately, did not receive attention in this industry. 

9. Fu, T., and Quan, M., "Gas-Temperature Measurement with Minimal Perturba- 
tion," NASA Tech Briefs, Spring 1983, p 305. This brief describes the 
method using two total heat flux calorimeters in contact with the gas 
and a third calorimeter that senses radiation only. The devices in con- 
tact with the gas are at differing temperatures. The measured gas tem- 
perature would be that at the calorimeter location, not out in the free 
stream. It would also seem important that the total calorimeters be 
physically located in such a way that each experiences a nearly equal 
heat transfer coefficient. 

10. Ghezzi, U., Zizak, G. , and Coghe, A., "Local Flame Temperature Measure- 
ments by Radiative Methods," AGARD Paper AD229, February 1978, p 35-1. 
Local flame temperature measurement at atomic fluorescence spectrometry. 
Line intensity ratios used. Gas stream seeded with gallium or indium or 
thallium. Typical errors of 2% to 4% of temperature were observed. 

11. Glawe, G.E., Johnson, R.C., and Krause, L.N. , "Intercomparison of Several 
Pyrometers in a High Temperature Gas Stream," Temperature, Its Measure - 
ment and Control in Science and Industry , Reinhold Publishing Corpora- 
tion, New York, 1962. (1) The potential and limitations of several types 
of pyrometers for gas temperature measurement are given. Discusses 
sources of errors for each of the techniques. 

12. Gottlieb, M., and Brandt, G.B., "Hot Spot Detection in Generators with 
Optical Fibers," ASME Paper 81-JPGC-PWR-18, October 1981. A joint EPRI- 
Westinghouse program led to the development of several methods of using 
optical fibers to detect the presence of thermal hot spots in remote 
locations in turbine generators. The fiber itself is used as the temper- 
ature sensor and also as the data transmission mechanism. One and two 
color pyrometry used. No test data. 

13. Gottlieb, M. , and Brandt, G.B., "Fiber Optic Temperature Sensor Based on 
Internally Generated Thermal Radiation," Applied Optics, October 1981, 
p 3408. This article is similar to the ASME paper, but it gives more 
detailed equations and actual test data using a germanium detector, 
starting at 140° C, and a silicon detector, starting at 240° C. A lead- 
sulfide detector was also proposed but not tested. 
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14. Goulard, R. , Mellor, A.M. , and Bilger, R.W. , "Combustion Measurements in 

Air-Breathing Propulsion Engines: Survey and Research Needs," Combustion 

Science, November 1976, p 195. General discussions of a broad range of 
combustion measurements. Extensive references given. Refers to Project 
SQUID workshop at Purdue University, May 1975, and published as "Combus- 
tion Measurements," R. Goulard, editor. Academic Press, New York, 1976. 

15. Hall, J.A. , "The Internal Temperature Scale," Temperature, Its Measure - 
ment and Control in Science and Industry , Proceedings of Symposia 2, 
1955, pp 115-139. 

16. Hill, W.E., and Dibelius, N.R. , "Measurement of Flame Temperature and 

Emittance in Gas Turbine Combustors," ASME Paper 70-GT-19, May 1970. An 
experimental method for measuring the temperature of luminous flames with 
a two-color pyrometer. Comparison with iridium versus iridium, 50% 

rhodium thermocouples. Emittance was calculated, various fuels were used 
from light through bunker "C." 

17. Kayukawa, N. , Aoki, Y. , and Ozawa, Y. , "Application of Light Polarization 
Techniques to the Generalized Line-Reversal Method for Gaseous Tempera- 
ture Measurements," Review of Scientific Instruments, November 1982, p 
1653. An improved line-reversal method, intended for fast transients, 
such as shock tubes. 

18. Kostkowski, Lee, Quinn, et al., "Intercomparison of the IPT568 above 
1064® C by Four National Laboratories," Paper 36, Temperature, Its Mea - 
surement and Control in Science and Industry , Volume 4, Part 1, Published 
by the Instrument Society of America, 1972. This paper is devoted to the 
high temperature end of the standards available. Correct and current 
values of all the melting point standards given. 

19. Lynnworth, L.C., and Carnevale, E.H., "Ultrasonic Thermometry Using Pulse 

Techniques," Temperature, Its Measurement and Control in Science and 

Industry , Volume 4, Part 1, published by Instrument Society of America, 
1972. This paper details the use of sonic velocity to measure gas tem- 
perature without immersed probes. Error sources mentioned. Transducer 
environmental sensitivity and delicateness a problem. 

20. Mossey, P.W., "Weighted Optical Temperature Measurement of Rotating Tur- 
bomachinery," U.S. Patent No. 3,696,678, assigned to General Electric 
Company. This patent has claims for a method of turbine blade surface 
temperature measurement that rejects interference from the hot turbine 
gasses that intermittently contain incandescent particulates. Another 
claim is made for a method to measure the gas temperature by holding and 
processing the peaks that occur when the same particulates are present. 

21. Rohy, D.A. , Duffy, T.E., and Compton, W.A. , "Temperature Measurements 
for Advanced Gas turbine Controls," AGARD Conference Proceedings, C57/ 
AD151, February 1975, p 27-1. Eleven concepts were analyzed. The three 
best concepts chosen were the beta-ray sensor (±2% of temperature at 
best), the optical immersed target radiation pyrometer, and the ultra- 
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sonic air-gap. The application was for a turbine inlet gas temperature 
flight measurement system. No details or references given on the optical 
immersed target pyrometer. 

22. Rohy, D.A., Seegall, M.I., and Compton, W.A. , "Turbine Inlet Gas Tempera- 
ture Sensor for Engine Control," Solar Division Research Report RDR-1748 
(RTD FDL TR 73-49), AD765193, February 1973. A beta-ray absorption tem- 
perature sensing technique was designed, built, and tested. See comments 
in above Rohy reference. 

23. (Sandia National Laboratories), "Ultrasonic Thermometer Measures High 
Temperatures," Machine Design, May 21, 1981, p 38. Thermometer measures 
high temperatures by sensing changes in the velocity of sound waves that 
travel along a thoriated tungsten wire. An accuracy of ±5° C to 600° C 
and ±1° C from 1000° to 2845° C is claimed. 

24. Stanforth, C.M., Engel, M.E., and Hall, R.L., "Hot Target Sensor for High 
Temperature Gas Measurements," General Electric Company TIS R73AEB483, 
December 1973. Describes the design, construction, test, and analysis of 
the Hot Target Sensor. The most promising target material was made from 
yttria-stabilized zirconia. Aspiration was used to maintain moderately 
high gas velocity through the porous target, bringing the target tempera- 
ture close to the gas temperature. 

25. Tichenor, D.A., and Dils, R.R., "A Fiber Optic Probe for Measuring High 
Frequency Temperature Fluctuations in Combustion Gasses," SAND83-8871, 
Sandia National Laboratories, Livermore, California, February 1984. This 
report compares an analytical solution to some test data in order to 
assess the frequency response of the sapphire fiber thermometer. 

26. Warshawsky, I., "Measurement of Rocket Exhaust-Gas Temperature," ISA 
Journal, Volume 5, No. 11, November 1958, p 91. A review of probe- and 
radiation-type temperature sensors, suitable for hot, high-velocity gass- 
es. Error sources and temperature ranges are given for each technique. 

27. Williamson, R.C., and Stanforth, C.M., "Measurement of Jet Engine Combus- 
tion Temperature by Use of Thermocouples and Gas Analysis," SAE National 
Air Transportation Meeting, Paper No. 690433, April 21-24, 1969. This is 
a comprehensive paper on high temperature measurement in gas turbine 
combustors. Error sources, system equations, and references are given. 


Refractory Materials 

1. Gryvnak, D.A., and Burch, D.E., "Optical and Infrared Properties of Alu- 
minum Oxide at Elevated Temperatures," Journal of the Optical Society of 
America, Volume 55, Number 6, June 1965, p 625. The absorption coeffi- 
cient of single crystal aluminum oxide (sapphire) is given versus temper- 
ature and versus wavelength. 

2. Larsen, D.C., and Adams, J.W. , "Property Screening and Evaluation of 
Ceramic Turbine Materials," Semiannual Interim Technical Report No. 12, 
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Contract No. F33615-79-C-5100, U.S. Air Force, by IIT Research Institute, 
Chicago, Illinois, May 1982. Strength Characteristics given for a wide 
variety of ceramics. 

3. Materials Selector, published by Materials Engineering (Journal), Penton/ 
IPC, Cleveland, Ohio, December 1984. Has tabular listings of significant 
materials properties of refractory metals and ceramics. 

4. Musikant, S., "Optical Materials ," Volume 6 in the Optical Engineering 
Series, Marcel Dekker, Inc., New York, 1985. The selection and applica- 
tion of materials for optical applications. 

5. (Union Carbide), "Optical Properties and Application of Cz Sapphire," 
Electronic Materials Division, San Diego, California, 1980. This bro- 
chure compiles the most useful optical, physical, chemical, and thermo- 
dynamic properties of single crystal aluminum oxide, grown by the 
Czochralski (Cz) process. 
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APPENDIX D - COMPUTER PROGRAMS 


D1 - RESP-96D 
D2 - WAT-CA96D 
D3 - A/D READ-MIO 


D4 - PM-TEMP4A 



APPENDIX D1 


COMPUTER PROGRAM RESP-96D 



2 REM PROGRAM 'RESP-56B ' 

4 REM date: 4-26-85 

5 REM CALC OF DETECTOR RESPONSE CURVES 

6 REM 96 ITEM FILES 

10 DIM Rfa8,8),R<18,100)iZTa3il00),2S<lS,100) 
16 GOTO100 

18 PRINT" INITIALIZING TABLES PLEASE WRIT" 

19 GOTO 100 

20 F0RU=1T018 

21 PRINTU 

22 F0RV=1T096 
24 A<U/V)=1 

26 NEXT 

27 NEXT 


100 PRINTCHRi:< 147) "»#*** RESP-96P MENU *mm>' 

102 print: PRINT 

104 PRINT" 1 - INPUT DATA 
106 PRINT" 2 - EDIT DATA FILE" 

103 PRINT" 3 - NORMALIZE DATA FILE " 

110 PRINT" 4 - DIVISION OF FILES " 

112 PRINT" 5 - SAVE ON DISK" 

113 PRINT" 6 - LOAD FROM DISK" 

114 PRINT" 7 - PRINT FILE 

113 PRINT" 8 - CALC. LIGHT PIPE TRANS" 

116 PRINT" 9 - CALC. DET. RESPONSE - CUT #1" 

117 PRINT" 10 - EDIT FILE HERDER " 

118 PRINT" 11 - READ DISK ERROR CODE" 

130 K»0 : INPUT"COMMAND" ; K 

132 IFK<1THEN100 
134 IFfOllTHEN100 

136 ONKOOTO200.260^400^500^ 600^ 700.- 800^900.- 1000.1100.- 1230 


200 PRINTCHR*< 147) "♦**#***!**;♦!*■ INPUT DATA **»♦#*#*#" 

202 print: PRINT 
206 PRINT 

208 GOSUB1208 ' 

209 PRINT 

210 Lf=""! INPUT" FILE NUMBER = ";L$ 

211 IFVAL<L$)>18THENGOTO210 

212 HH»a"": INPUT"ENTER LETTER DESIGNATION OF FILE NUMBER ",:HH$ 

213 L=VALaS):QQM"": INPUT" INIT TRBLE IF V ";QQ* 

214 GOSUB2100 

216 A$<L, 1 )=L» : A»<L. 3)»"96" : R$<L. 2)="RESP-"+L»+HHf : U=. 3 

217 IFQQ»0"V"THEN219 

213 F0RQ-1T096 : 2T<L/ Q)=0 : ZSCL . Q ) =0 : NEXT 

219 P=1 : INPUT"STRRT AT 1= ";P 

220 PRINT SPRINT" INPUTS ARE:" 

221 PRINT" VALUE » INCHES SCALE = VOLTS/INCH" 

222 PRINT" OR" 

223 PRINT" VALUE » VOLTS SCALE = MULTIPLIER" 

224 PRINT :PRINT"ENTER E+LINE HO. TO MAKE CORRECTION" 

225 PRINT :PRINT"LINE WAVELEN VALUE SCALE ENTRY 

226 FORI=PT096 

223 TJ="":iNPUT"VALUE a"JTJ: 

230 IFLEFT#<T$.1)="E"THEN234 

231 PRINT"n"TABa5); :1NPUT"SCALE =";S$ 

232 IFLEFT$<S*.l)a"E"THEN254 

233 IFTSa""THEN238 

234 W2=W+<.01H«I-1)) :T=VAL<T*) :S=VAL<S:f) 

235 A<L.I)=T»S 

- 237 ZT<LiI)aT:ZSa.I)»S 
233 PRINT"n 

239 PRINT"n"TAB<2); SPRINTI 

240 PRINT"n"TAB<6); :PRINTW2 

241 PRINT"n"TAB<14)J :PRINTZT<L. I) 

242 PRINT"n"TAB<22); sPRINTZSCL. I) 

243 PRINT"n"TAB<30); :PRINTA<L.I) 

248 NEXT 

249 FORK* 1 TO4000 : NEXT 

250 GOTO100 

254 IFLEN<TfX2THEN258 

256 P*VAL<MID*<T* , 2) ) : IFP>0THEN22S 

233 LLS»"":iNPUT"MENU IF V ";LL4:: IFLLf*"V"THEN100 

259 IHPUT"RESTART AT LINE ";p:G0T0226 

260 PRINTCHR$<147) 



261 PR I NT "***#»*•♦!# EDIT DRTH FILE : PRINT: PRINT 

262 PRINT'TILE NO » "L: INPUT'THflNGE TO ";L 

263 IFL>10THENPRINT:GOTO262 

264 U=1 : INPUT"STflRT AT "JU 

265 PRINT" VALUE SCALE ENTRV" 

266 FORP=IJTOU+10 

267 PRINTP" "2T<L,P)" "2S<L^P)" "ACL/P) 

268 IFP=96THEN270 
263 NEXT 

270 I=0!INPUT"CHANOE LINE NO "JI 

271 IFK1THEN231 

272 IFKUTHEN270 

273 IFI>U+10THEN270 

274 K2=0!INPUT"CHRNGE VALUE ";K2 

275 IFK2>0THENZT<Li n*K2 

276 K3=0!INPUT"CHANGE SCALE ";K3 

277 IFK3>0THEN2S<L,I)=K3 

278 R<L/I5=K2*K3 
273 IFK2>0THEN265 

280 IFK3>0THEN265 

281 LL#=""! INPUT’TIENU IF V "ILLf 

282 IFLL5="V"THEN100 

283 U=U+1 I ! IFU>36THEN264 

284 G0T0266 
283 GOTO 100 


400 PRINTCHR$<147)")K***m^* NORMALIZE DATA FILE mmmm” 
402 PRINT! PRINT 

404 INPUT"ENTER SOURCE FILE NO.’SL 

405 INPUT"ENTER DESTINATION FILE NO.";LD 

406 PRINT: PRINT 

407 PRINT"SOURCE FILE HEADER: PRINT 
410 F0RI=1T03:PRIHTI" '•R*<L/ 1) : NEXT 
420 print: PRINT 

422 LL»="":iNPUT"MENU IF V";LLf 
424 IFLL*="V"THEN100 
430 A<L.0)=0 

432 F0RI=1T036 

434 IFACL. I)>A<Li0>THENACL/0>=A(L.I) 

436 NEXT 

433 PRINT :PRINT"MAXIMUM VALUE » "A<Li0) 

440 IFACL/0)=0THENPRINT"NO DATA FOUND IN FILE" : G0T0472 
430 F0RI=1T036 

452 n<LD.I)=RCL,n/'R<L,0):PRINTI" "RCLD, I) 

454 NEXT 

460 Rf(LD,l)=STRi:<LD) 

462 RS:<LDiZ)=RJ<L>2)+"-N" 

464 R«<LDi3)=Rf<Li3) 

466 R*CLD,4)=H*<L,4) 

470 PRINT’TILE "LD" IS NOW THE NORMALIZED VERSION OF CURVE "L 
472 FOR I = 1 TO4000 : NEXT 
474 GOTO 100 


500 

502 
504 
506 

503 
503 

510 

511 

512 

513 

514 
516 
520 
522 
524 
526 
528 
530 
532 
540 


PRiNTCHRfa47)“***Hi»m DIVISION OF FILES 
PRINT : PRINT !L=0:INPUT"FILE NO. OF DENOMINRTOP=" ; L 
PRINT!LLi:="":iNPUT"ENTER V FOR RECIPRICAL";LL:f 
IFLL*»"V"THEN550 

print: INPUT"FILE H0. OF NUMERATORa";LN 
PRINT :PRINT"FILE "LN" IS TO BE DIVIDED BV FILE "L 
print: INPUT"DESTINRTION FILE N0.*";LD 
PRINT"FILE "Rf<LN/ 2) "DIVIDED BV "Af<L.2>; 

PRINT"WILL BE LOCATED IN FILE "LD: PRINT 

H*<LD,4)=Rf<L,2)+" / "+Ai:<LfL2) : ACLD. 1 >=STR$<LD) :RJ:<:LD. 3)=96 
LLf=""!INPUT"MENU IF V";LL$ 

IFLLJ="V"THEN100 

F0RI=1T036 

G2=R<L>I) 

IFRCL, I)«0THENO2=lE-20 
R(LD/I)=R<LN,nXG2 
PRINTI" "RCLD^n 
NEXT 

FORI=1TO4000:NEXT 
GOTO 100 



550 print: IHPUT'TESTIN flTION FILE NO.="JLD 

552 PRINT"RECIPRICnL OF FILE ''L"IS TO BE STORED IN FILE "LD 
554 PRINT'TILE "L" IS "Flf<L>2) 

556 LL3:="":iNPUT"riENU IF V'JLL* 

553 IFLL«a"V"THEH100 

566 F0RI=IT096 

567 02=fl<L.I) 

563 IFfl<L/I)=0THENG2*lE-20 

569 fl<LD>n=l/G2 

572 PRINTI" ''FKLD/I) 

574 NEXT 

578 FORI=1TO4000:NEXT 

579 flf <LD, 1 )=STRf <LD) :R$<LD. 2)®fl* <L^ 2)+"-R" : flT<LD/ 3>=95 
■IRR RnTni00 


600 PR I NTCHR$< 147 SAVE TO DISK mmm%” 

601 print: PRINT :PRINTK$” "B*:PRINT 

602 GOSU61604 

620 Lf='"' : INPUT"ENTER FILE NUMBER 

622 HH$="":iNPUT"ENTER LETTER DESIGNATING VERSION ";HHT 

631 L=VflL<L$) 

632 PRINTCHRf<147) 

633 PRIHT"FILE NAME IS RESP-'SLTJHHf 

634 print: PRINT 

636 PRINT"HEADER FOR FILE % "L*" * IS:" : PRINT 

637 F0RI=1T03 

633 PRINT"Af(''r*) "Af<L,I) 

639 NEXT 

640 LLa0:iNPUT"CHANGE LINE ";LL 

641 IFLL<1THEN643 

642 INPUT“ENTER CHANGE" JAJC Li LL) :G0T0632 

643 IF V ";LLT: IFLLT="V"THEN10O 

646 LL#="":iNPUT"DELETE OLD FLIE IF V ";LLf 

647 IFLLf="V"THEH630 

648 K$»"RESP-"+L$+HH$ 

649 L=VAL<Lf) 

650 0PEN15,8il5 

651 0PEH2i8i2i "0: "+Kl:+"iSi W" 

652 INPUT#15iAfiBfiC$iDf 

653 IFA4:»"00"THEH657 

654 PRINTCHR»a47> : PRINT : PRINT : PRINT : PRINT = PRINTTRB< 12>BS 

655 CL0SE2 : CL0SE15 : FORI=1TO3000 : NEXT : GOTO600 
657 F0RU=1T08 

653 PRINT#2i"R"+A$<LiU) 

659 NEXT 
665 F0RU=1T096 

667 PRINT#2iA<LiU) 

668 PRINTU" "R(LiU) 

669 NEXT 

670 PRINT#2i "REND" 

672 CL0SE2: CLOSE 15 
674 GOTO100 

680 0PEN15i8i 15i "S:RESP-"+Lf+HHf 

632 CLOSEIS 

684 G0T0648 

690 Bf="":GOTO100 


700 PRINTCHR*<147>"r4(**#:tfi*! LORD FROM DISK ***4:****'' 

702 PRINT! PRINT :PRINTKf" "BI!:PRINT 

703 GOSUB1604 

704 Lf="":iNPUT"ENTER FILE NUMBER "ILt 

705 HHf="":iNPUT"ENTER LETTER DESIGNATION OF VERSION ";HHf 

706 LL$="":iNPUT"MENU IF V ";LLf 

707 IFLL*="V"THEN100 
710 KlP="RESP-"+LT+HHf 
713 0PEN15i8il5 

720 0PEN2i 8i 2i "0 : "+K^+" / Si R" 

721 INPUT#15iR^.BTiC«iD$ 

722 IFRJ:="00"THEN725 

723 PRINTCHR$<147)"M!WW)!I!!li!n’RB<12)BT 

724 CL0SE2 : CLOSEIS : FORI= 1TO3000 : NEXT : GOTO 1 00 



725 F0RU=1T08 

726 INPUT#2iC$<U) 

727 PRINTC$<U) 

723 NEXT 

73G CL0SE2 

732 L=VflL<MID»<Cfa)/2>):PRINT"FILE N0.= "L 

734 PP=VRL<riID5<Cf<3).2)):PRINT"FILE CONTFIINS "PP" ITEMS" 

740 LL»=""!INPUT"HEHU IF V ”;LL$ 

741 lFLLf=»V"THEN100 

750 0PEN2. 8, 2, "0: "+K$+" , S, R" 

752 F0RU=1T08 

753 INPUT#2,Cf 

754 R$a/U>=MI13fCW,2) 

755 NEXT 

760 FORIJ=lTOPP 
762 INPUT#2,R(LiU) 

764 PRINTU" "fl<L>U) 

766 NEXT 

770 CL0SE2:CL0SE15 

771 FORUs1TO500:NEXT 
730 B$»""!GOTO100 


800 PRiNTCHRfa47)":mf**i» PRINT FILE mm'mmmmw." 

802 print: PRINT 
804 GOSUB1600 

806 B$<n="FILE NUMBER 

807 B$<2)="FILE NAME 

808 B#<3)="BflTR ITEMS IN FILE 

810 LL=0:INPUT"ENTER FILE NUMBER ";LL 
012 P=1 : INPUT"STRRT LIST RT N= ";P 
820 OPEN4,4,0:CMI)4 

330 PRINT"PRINTOUT OF DRTR FILE PROGRAM * RESP-69R * " 

332 PRINT : PRINT"BESCRIPTIOH = " 

833 print: print: PRINT 

834 F0RI=1T08:PRINTI" "BJ<I);fl$(LL. D" " :NEXT 
836 print: PRINT" ITEM IVL VALUE 

333 F0RI=PT096 

340 PRINTI" “.3+<.01H«I-l>)" "R<LL, I) 

845 NEXT 

848 PRINT#4:CL0SE4 
850 GOTO 100 


300 PRINTCHRi:<147)"»#*»!»f* CALC LIGHT PIPE TRANSMISSION »**»:♦;**" 
902 print: PRINT 

904 PRINT-TOP * FILE 17 / FILE 11 DESTINATION FILE 4 

905 PRIHT"BOT- * FILE 13 / FILE 12 DESTINATION FILE 4" 

910 INPUT-LIGHT PIPE FILE";N 

911 INPUT "DETECTOR FILE";D 
920 RJ<4,l)a"4" 

922 R$<4i3)=*96" 

924 R3!<4, 2)«A*<N, 2>+'V"+Ai:<D, 2) 

926 Rf<4.4)*"FIRST CUT TRANSMISSION CALC. " 

940 F0RI=1T096 

942 IFRCD, I)=0THENA<4, I)=0:GOTO946 
944 R<4,I)bR<N/I)/R<D,I) 

946 NEXT 
962 GOTO 100 


1000 PR1NTCHR$< 147) "***:+:*** CALC DET. RESPONSE - CUT#1 
1002 print: PRINT 

1004 C1=0: INPUT-DETECTOR FILE NO = ";C1 

1005 PRINT : PRINT : PRINT"RB=< IE/IA)*RB" 

1006 PRINT" lA IS EG4G OUTPUT FILE NO. 13" 

1008 IFC1=12THENH1=9 

1009 IFC1=11THENH1=10 

1020 FORI=1T096 

1021 IFAa3,I)=0THEN1030 

1022 R<Hl,l)=<R<Cl/I)/fla3/I))*Ra4,I) 

1028 PRINTI" -A<H1.I) 

1030 NEXT 



r 


1032 FORU=1TO1000=NEXT 

1034 flf<Hl/l)=MID$<STR$<Hl)>2) 

1035 fl$<Hl;3)=“96" 

1036 flJ<Hl , 2)=''RESP-"+riID3:<STRi: (HI ) , 2) 

1037 flf(H1.4)="CUT ONE CRLC OF RESPONSE OF DETECTOR" 

1040 B$a)=" s FILE NUMBER" 

1041 B$<2)*" a FILE NOME" 

1042 B$(3)a" a DATA ITEMS IN FILE" 

1060 PRINT'TILE HEADER : " 

1061 FORIalTOS 

1062 PRINTI" "Rf<Hl/I)" "B^a) 

1064 NEXT 

1068 Sa0.‘INPUT"CHANGE LINE'SS 

1069 IFS>0THEN1071 

1070 GOTO 100 

1071 K$a««:iHPUT"MAKE CHANGE" JK$ 

1072 IFK#>""THENAf<Hl.S)aKf 

1073 GOTO 1060 


1100 PRINTCHRf< 147) "********* EDIT FILE HEADER mmmm" 

1102 print: PRINT 

1104 C5(l)a"HCFILE NUMBER]!)" : CS<2)a"«CFILE NAME]!)" 

1105 C$<3)a*'||CFILE LENGTH]!)" :Cf<4)a"ilCTVPE OF DATA]!)" 

1106 C*<7)a«'':CJ<8)="" 

1107 C*(3)="HCRESERVED FOR COMPUTER ENTRV]!)" 

1103 C$<6)a«r»C0PERAT0R NOTES]!)" 

1110 L=0:iNPUT"ENTER FILE NUMBER" ;L 

1120 F0RI=1T08 

1121 PRINTCJCI) 

1122 PRINTI" "A$<LiI) 

1124 NEXT 

1130 L8=0:INPUT"CHANGE LINE";L8 
1132 IFL8>0THEN1140 
1134 GOTO100 

1140 LL$a"":iNPUT"MAKE CHANGE "JLL* 

1142 A»<L,L8)aLL# 

1146 GOTO1120 


1206 REM **»»»*# FILE NUMBER LIST *mmm 

1208 PRINT"ENTER 11 FOR TOP DET. DATA" 

1209 PRINT"ENTER 12 FOR BOT DET. DATA" 

1210 PRINT"ENTER 13 FOR EG.£G DET. DATA" 

1211 PRINT"ENTER 14 FOR EGfiG DET. MFG.RESP." 

1212 PRINT"ENTER 15 FOR TOP DET. MFG.RESP." 

1213 PRINT"ENTER 16 FOR BOT DET. MFG.RESP." 

1214 PRINT"ENTER 17 FOR TOP DET. + LIGHT PIPE DATA" 

1215 PRINT"ENTER 18 FOR BOT DET. + LIGHT PIPE DATA" 
1219 PRINT 

1230 RETURN 


1280 PRINTCHRi:<147)"**»!^**** DISK 
1231 PRINT:pRINT:pRINT:pRINT 

1282 0PEN15,3,1S 

1283 INPUT#15iAi:iB*iC*,Df 
1234 CL0SE15 

1236 PRINTA$,Bf/C«,Df 

1288 INPUT"MENU ON C.R. ";X 

1290B*a“«!GOTO100 


1600 

1604 

1605 

1606 

1607 

1608 

1609 

1610 
1611 
1612 


REM *»**»*** LIST OF FILES ********* 

PRINT :PRINT"FILE 1 - TOP DET RESP 
PRINT'TILE 2 - BOT DET RESP" 

PRINT-FILE 3 - LIGHT PIPE RESP" 

PRINT'TILE 4 - MISC. A 

PRINT"FILE 5 - MISC. B 

PRTNT"FILE ^ - flISC- c 

PRINT-FILE 7 - BOT DET INTERMEDIATE CALC #2" 

PRIhI-FILE 3 - TOP DET INTERMEDIATE CALC #2" 

PRINT-FILE 9 - BOT DET INTERMEDIATE CALC #1" 
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1613 PRIHT"FILE 

1614 PRIHT'TILE 

1615 PRINT'TILE 

1616 PRINT"FILE 

1617 PRINT'TILE 

1618 PRINT-'FILE 
1618 PRIHT"FILE 

1620 PRINT’TILE 

1621 PRINT'TILE 
1629 RETURN 


10 - TOP DET INTERMEDIRTE CRLC #1" 

11 - TOP DET DRTR" 

12 - EOT DET DRTR" 

13 - EG40 DET DRTR" 

14 - EO&O NFG. RESP." 

15 - TOP DET MFC. RESP. DRTR" 

16 - EOT DET riFG. RESP. DRTR" 

17 - TOP DET - LIGHT PIPE DRTR" 

13 - EOT DET - LIGHT PIPE DRTR" 


2099 REM RSSIGN TITLE TO FILE HERDER 

2100 IFLfa"H"THB'IR*<Li4)="T0P DET DRTR" 

2102 IFL#="l2"THENfl»<L/4)="B0T DET DRTR" 

2103 IFL#="13"THENR$<L.4)="EG«G DET DRTR" 

2104 IFLJ:="14"THENR*<L.4)s"E0*G MFO. RESP CUR’/E" 

2105 IFL*="15''THENflf<Li4)="TU0 COLOR DET # TOP * MFG. RESP. DRTR 

2106 lFL*a"16"THENRf<L>4)="TW0 COLOR DET * EOT * MFG. RESP. DRTR 

2107 1FL#="17"THENR$<L/4>="T0P DET. m LIGHT PIPE DRTR " 

E108 IFL#="18"THENR*<L/4)="B0T DET. m LIGHT PIPE DRTR " 

2120 RETURN 



PROGRRM ^RESP-9611 ^ : - 12/^03^85 - REMflRKS_ 

THE nniN FUNCTIONS OF THE PROCRRM RRE: 

INPUT DRTR. 

SAVE DRTR ON DISK. 

PRINT DRTR ON PAPER. 

SOME OF THE tIRTH OPERRRTIONS RRE BEST DONE IN THE IMEDIRTE NODE. 
RN EXAMPLE IS THE PRODUCT OF TWO TABLES: 

F0RIelT09S:fl<5/ I>sRa3» I)*Ra4. 1? iNEXT 


THE SUPPORT PROGRAMS RRE USED TO OPERATE THE I/O PORT. 

MADE FOR fl SET OF 18 FILES. SOME RRE LISTED STARTING RT LINE 
1600 IN THE PROORRM. 


ERCH FILE MUST HAVE H DIFFERENT NAME. WHEN R FILE IS 

WHEN RUN THESE PROGRAMS DISPLAY R LIST OF REM STRTEMENTo FOR MAIN 

PROGRAMS #1 #2 a «3. 

THE REM PROGRAMS CAN BE LISTED OR RUN. 

WHEN RUN THE PROGRRM GIVES THE OPTION OF DISPLAY OR HRRD COPY. 


ALL FILES ARE SEQUENTIRL STARTING WITH R HERDER 8 LINES LONG 

FOLLOWED BY R 96 LINE DRTR SET 

LINES 3 THROUGH B OF THE HERDER MRY BE USED FOR UoER NOTES. 
LINE 4 IS RESERVED FOR PROGRRM USE. 

LINE 3 INDICATES THE HUMBER OF DATA ITEMS IN THE FILE. 

LINE 1 ll ThI ^lI NUnIeR - (STORAGE LOCATION IN THE PROGRRM), 


WAVELENGTH - IS RELATED TO DHTA 
W2=0.3+(.01*CI-1) 


LINE NUMBER RS FOLLOWS •’ 
WHERE I IS THE LINE NUMBER. 


WHEN ENTERING DATA TWO ENTRIES ARE RSKED FORE 

1 - SCALE IN VOLTS PER INCH 

2 - READING IN INCHES 

THIS SCHEME IS INTENDED TO FRCILITRTE ENTRY OF GRRPHED DRTR. 


IF RN INCORRECT ENTRY HRS BEEN MADE - ENTER THE CHRRRCTER 'E' IN 
THE"pROGRRM^WILl'''r^ for the NUMBER OF THE ENTRY TO BE CORRECTED. 
IT MAY BE EASIER TO MAKE CORRECTIONS WHILE IN THE EDIT FUNCTION. 
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P RDGR RM ^RE5P-9 6D^ j- 12/09/85 - REMARKS 

mmmw. calculation of tetector response m*mm* 
assumptions: 

FILE 13 IS THE OUTPUT OF A SPECTROMETER AS READ BV A REFERENCE 
DETECTOR. 

FILES 11 & 12 IS THE OUTPUT OF THE SPECTROMETER AS REAHD BV 
TOP AND THE BOTTOM DETECTOR. 

FILE 14 IS THE SPECTRAL RESPONSE CURVE OF THE REFERENCE DET. 

THE CALCULATION IS AS FOLLOWS •* 

(CTOPCOR BOT. JDETVCREF DET.)?*(REF DET RESPONSE) 

((FILE IIVCFILE 13))»FILE 14) 

SEE LINE 1022 


FILE numbers: 

IF THE SOURCE FILE IS NO. 11 THEN THE DESTINATION FILE IS 
NO. 10 

IF THE SOURCE FILE IS NO. 12 THEN THE DESTINATION FILE IS 
NO. 9 

TABLE OF FILE NUMBER ASSIGNMENTS wmmmm 

1 TOP DET. - NORMALIZED SPECTRAL RESPONSE. 

2 EOT DET. - NORMALIZED RESPONSE. 

3 NOT ASSIGNED. 

4 LIGHT PIPE - NORMALIZED RESPONSE. 

5 NOT ASSIGNED. 

S NOT ASSIGNED. 

7 NOT ASSIGNED. 

3 LIGHT PIPE TRANSMISSION. 

9 EOT DET. - CALCULATED SPECTRAL RESPONSE. 

10 TOP DET. - CALCULATED SPECTRAL RESPONSE. 

11 TOP DET. - OUTPUT OF SPECTROMETER. 

12 EOT DET. - OUTPUT OF SPECTROMETER. 

13 REFERENCE DET. - OUTPUT OF SPECTROMETER. 

14 REFERENCE DET. - SPECTRAL RESPONSE - MANUFACTURERS CAL. 

15 TOP DET. - SPECTRAL RESPONSE - MANUFACTURERS CAL. 

16 EOT DET. - SPECTRAL RESPONSE - MANUFACTURERS CAL. 

17 DATA - TOP DET. READING SPECT. OUTPUT THRU LIGHT PIPE. 

18 DATA - BOT DET. READING SPECT. OUTPUT THRU LIGHT PIPE. 

note: suffix 'N' ON A FILE INDICATES NORMALIZED DATA. 



m*m REMARKS FOR PROGRAM RESP-96D mm:* 


' GOTO8000 10/27/85 

3 PROGRAM RESP-96D-R3 

4 . THE FOLLOWING REMARKS ARE BV LINE NUMBERS WHICH 

5 . CORRESPOND TO THOSE IN THE PROGRAM. 

6 . 

10 DIMENSION FOR 18 FILES OF 100 ITEMS EACH 
IS GOTO MAIN MENU, BVPASS LINES 13 THRU 27. 

18 

20 . . SETS ALL VALUES OF ALL . 

21 . .18 TABLES TO 1.0 . 

22 • 

24 ! ! EACH TABLE IS A FILE, i 

26 . . 

27 


100 . main menu 

102 . 

104 . 

106 . 

108 . 

110 . 

112 . 

114 . 

115 . 

116 . 

117 . 

118 . tttttmt END OF MENU tttttTttTtt 

130 . . ■ 

132 CHECK THAT COMMAND IS WITHIN LIMITS. 
134 . t T t t t t 

136 STARTING ADDRESSES OF MENU ITEMS. 


200 . mmM MANUAL INPUT OF DATA ***** 

202 . 

206 . 

203 SUBROUTINE DISPLAYS FILE NUMBERS AND NAMES. 


209 

210 
211 
212 

213 

214 
216 

217 

218 

219 

220 
221 
222 

223 

224 
-225- 

226 

229 

230 

231 

232 

233 

234 

235 

237 

238 

239 

240 

242 

243 
243 
249 


ENTER FILE NUMBER. 

CHECK THAT FILE NO. IS WITHIN LIMITS. 

ENTER UNIQUE NAME OR CODE DESCRIPTION. 

INIT TABLE ONLY TO ERASE PREVIOUS FDATA IN FILE. 
PUTS DATA TYPE IN HEADER LINE #4 FOR FILES 11 THRU 
LORDS FILE HEADER LINES 1/ 2, & 3. 

SKIP NEXT LINE IF FILE INIT NOT WANTED. 

SET ALL VALUES TO ZERO. 

SELECT SUBSCRIPT OF STARTING ENTRY. 


13. 


P IS THE STARTING SUBSCRIPT. 

VALUE IS A RAW READOUT ENTRY. 

CHECK TO SEE IF ENTRY IS A REQUEST FOR A CORRECTION. 

SCALE IS A MULTIPLIER OR CONVERSION NUMBER. 

CHECK TO SEE IF ENTRY IS A REQUEST FOR A CORRECTION. 

IF VALUE NOT ENTERED - SKIP TABLE ENTRY AND PRINT TABLE CONTENTS. 

CALC. WAVELENGTH 

CALC ••ENTRY- 

ENTER DATA INTO TABLE. 

DISPLAY SUBSCRIPT 
DISPLAY VALUE <0N ABOVE LINE). 

DISPLAY VALUE <0N ABO^v’E LINE). 

DISPLAY SCALE 
DISPLAY 


"TIS CRSR UP 
"TIS CRSR UP 


DELAY TO PERMIT READING OF LAST ENTRY. 


250 GOTO MAIN MENU. 

254 SKIP CORRECTION ROUTINE IF SUBSCRIPT NOT SPECIFIED. 
256 SET P = TO LINE NUMBER TO BE CORRECTED. 

253 ESCAPE TO MENU IF DESIRED. 

259 . 

260 CLEAR SCREEN. 



261 #******.4!#* EDIT DfiTfl FILE -Umm-mM 

262 SELECT RHOTHER FILE IF WANTED. 

263 CHECK INPUT. 

264 INPUT STARTING SUBSCRIPT <LINE NUMBER). 

265 . 

266 DISPLRV 10 LINES OF FILE. 

267 . 

263 CHECK FOR END OF FILE. 

269 

270 ENTER NUMBER OF LINE TO BE CHANGED. 

271 . 

272 BOTTOM OF DISPLRV. CHECK FOR VALID LINE HUMBER. 

273 TOP OF DISPLRV. 

274 . 

275 CHANGE VALUE ONLV IF ENTRV HRS BEEN MADE. 

276 . 

277 CHANGE VALUE ONLV IF ENTRV HRS BEEN MADE. 

273 RECALCULATE "ENTRV" 

279 DISPLRV NEW VALUES IF ENTRV HRS BEEN MADE. 

280 DISPLRV NEW VALUES IF ENTRV HAS BEEN MADE. 

281 . 

282 . 

283 ^T UP TO DISPLRV NEXT SET OF LIST. 

284 . 

283 GOTO MAIN MENU. 


400 ******** NORMALIZE DATA FILE 
402 . 

404 NUMBER OF FILE TO BE NORMALIZED. 

405 NUM. OF FILE THAT WILL BE THE NORMALIZED VERSION. 

406 . 

407 . 

410 . 

420 . 

422 ESCAPE IF DESIRED. 

424 . 

430 INIT "MAX. VRL." 

432 . . • 

434 SEARCH FOR MAX. VRL. 

436 . 

433 . 


440 . 

450 NORMALIZE TABLE. 

452 . 

454 . 

460 PUT FILE NUM. IN HEADER LINE #1. 

462 HDD "-N" TO FILENAME IN HERDER LINE #2. 

464 ENTER NUM. OF ENTRIES IN NEW FILE. 

466 SET HERDER LINE #4 OF NEW FILE = TO OLD FILE. 
470 . 

472 TIME DELAV. 


500 MmmmMm division of files mmmmm* 

502 . 

504 . 

506 . 

508 . 

509 . 

510 . 

511 . 

512 . 

513 LINE #4 OF NEW FILE HEADER TELLS WHERE DATA CAME FROM. 

514 ESCAPE IF DESIRED. 

516 . 

520 START DIVISION ROUTINE. 

522 . 

524 PREVENT DIVISION BV ZERO ERROR. 

526 DIVISION. 

528 DISPLRV RESULT. 

530 . 

532 TIME DELAV. 

540 GOTO MAIN MENU. 


554 • 

556 ESCAPE IF DESIRED. 


558 . 

566 START RECIPRICAL ROUTINE. 

567 . 

563 PREVENT DIVISION BV ZERO. 


56S . 

572 . 

574 . 

578 . 

579 INPUT DESTINATION FILE HERDER. 

580 GOTO MAIN MENU. 


600 SAVE FILE ROUTINE 

601 PRINTPRINTPRINTPRINTPRINTPRINT 

602 DISPLRV LIST OF FILES. 

620 FILE HUM. TO BE SAVED. 

622 . 

631 . 

632 . 

633 DISPLRV FILE NAME. 

634 . 

636 . 

637 DISPLRV FILE HERDER. 

633 . 

639 . 

640 EDIT HERDER LINE OPTION. 

641 . 

642 EDIT LINE. THEN REDISPLRV HERDER. 

643 ESCAPE IF DESIRED. 

646 IF FILE OF SAME NAME IS ON DISK - DELETE IT. 

647 . 

643 ASSEMBLE FILE NAME. 

649 GET FILE NUMBER AS INT. 

650 OPENCOMMRND CHANNEL ON DISK. 

651 OPEN FILE. 

652 READ DISK STATUS CODES. 

653 IF NO PROBLEM - CONTINUE 

654 IF PROBLEM --DISPLRV IT AND CLOSE FILES. 

655 START OVER. 

657 SAVE FILE HERDER. 

653 HDD CHARACTER TO EACH STRING IN CASE OF NULL STRING. 
659 . 

665 START DATA SAVE. 

667 . 

663 DISPLRV DATA. 

669 . 

670 END OF FILE MARKER. 

672 CLOSE FILES. 

674 GOTO MAIN MENU. 

680 DELETE FILE ROUTINE. 

632 DELETE FILE ROUTINE. 

684 GOTO SAVE FILE ROUTINE. 

690 CLEAR ERROR MESSAGE - GOTO MAIN MENU. 


700 »*»***» GET FILE FROM DISK ROUTINE Ht****** 

702 . 

703 PRINT LIST OF FILES. 

704 ENTER FILE NUM. TO BE LOADED (FROM DISK) 

705 . 

706 ESCAPE IF DESIRED. 

707 . 

710 ASSEMBLE FILE NAME. 

718 OPEN ERROR CHANNEL. 

720 OPEN FILE. 

721 READ DISK STATUS. 

722 IF NO ERROR/ CONTINUE LORD. 

723 DISSPLHV ERROR. 

724 CLOSE FILES - GOTO MAIN MENU. 



723 INPUT FILE HERDER FOR INSPECTION 

72S FIRST CHAR. IN STRING IS NOT REMOVED SO THAT NULL STRINGS ARE ODVIOUS. 
727 . 

723 ■ 

730 CLOSE FILE. 

732 CHECK FILE NUMBER. 

734 CHECK HUMBER OF ITEMS IN FILE. 

740 ESCAPE IF DESIRED. 

741 . 

750 INPUT HEADER. 

753 . 

734 REMOVE FIRST CHAR. OF STRING. 

755 . 

760 INPUT DATA. 

762 . 

764 . 

766 . 

770 CLOSE FILE AND COMMAND CHAN. 

771 DELAV. 

790 INIT ERROR MESSAGE - GOTO MAIN MENU. 


800 ******* HARD COPV OF FILE ********802 . 
804 DISPLAY LIST OF FILES. 

806 USED BELON. 

807 USED BELOW. 

803 USED BELOW. 

810 FILE TO BE PRINTED. 

312 CAN START AT SUBSCRIPT OTHER THAN ONE. 

820 OPEN PRINTER. 

830 .• 

832 . 

833 LINES. 

834 PRINT DESCRIPTION OF FILE. 

836 . 

333 PRINT DATA. 

840 LINE ~ WAVELENGTH — VALUE. 

843 • 

348 CLOSE PRINTER. 

850 GOTO MAIN MENU. 


900 ******** CALC LIGHT PIPE TRANS. ********* 

902 . 

904 . 

903 . 

910 INPUT ONE OF ABOVE TWO SOURCE FILES - DETECTOR OUTPUT THROUGH FIBER. 

911 INPUT FILE NO. FO DETECTOR USED FOR ABOVE MEAS. 

920 INPUT FILE HEADER INFO. 

922 INPUT FILE HEADER INFO. 

923 INPUT FILE HEADER INFO. 

924 INPUT FILE HEADER INFO. 

926 INPUT FILE HEADER INFO. 

940 CALCULATE TRANSMISSION. 

942 . 

944 . 

946 . 

962 GOTO MAIN MENU. 


1000 *** FIRST STEP CALC. OF DET. RESPONSE *** 
1002 . 

1004 . 

1005 CALCULATION TO BE MADE. 

1006 . 

1008 ESTABLISH DESTINATION FILE. 

1009 ESTABLISH DESTINATION FILE. 

1020 CALC RESPONSE. 

1021 . 

1022 . 

1028 DISPLAY RESULT. 

1030 . 



'01 


1032 TI,ME DELFiV TO VIEW LAST HUM. 

1034 INPUT TO DESTIRHTIOH HEADER. 

1035 INPUT TO DESTIRHTIOH HERDER. 

1036 INPUT TO DESTIRHTIOH HERDER. 

1037 INPUT TO DESTIRHTIOH HERDER. 

1040 DISPLAY INFO. 

1041 DISPLAY INFO. 

1042 DISPLAY INFO. 

1060 DISPPLRY DESTINATION FILE HERDER. 

1061 . 

1062 . 

1064 . 

1068 EDIT HERDER OPTION. 

1069 . 

1070 GOTO MAIN MENU 

1071 CHANGE CONTENTS OF LINE. 

1072 IF LINE HAS BEEN CHANGED - ENTER NEW CONTENTS. 

1073 REDISPLAY FILE HERDER. 


1100 #»**»■* EDIT ANY FILE HERDER mnm 
1102 CREATE HEADER DISPLAY TITLES. 

^ 1104 . "ir IS COLOR CHANGE - CTRL 6 

1105 . PRINTPRINT 

1106 . 

1107 . 

1108 . 

1110 . 

1120 DISPLAY HERDER. 

— 1121 . LINE TITLE. 

1122 . HERDER LINE CONTENT. 

1124 . 

1130 . ' 

_ 1132. 

1134 IF NO LINE INDICATED ABOVE/ GOTO MAIN MENU. 

1140 CHANGE LINE. 

1142 . 

1146 REDISPLAY HERDER. 


1206 *m SUBROUTINE -TO DISPLAY LIST OF FILES- 
1203 . 

1209 . 

1210 . 

1211 . 

1212 . 

1213 . 

1214 . 

1215 . 

1219 . 

1250 . 


1231 *m ROUTINE TO CHECK DISK ERROR CODE 

1232 OPEN COMMAND CHANNEL. 

1233 READ CODES. 

1234 CLOSE CHANNEL. 

1236 DISPLAY STATUS INFO. 

1293 . HOLDS DISPLAY UNTIL READY TO GO OH. 
1290 CLEAR STATUS MESSAGE. -• GOTO MAIN MENU. 


1600 . 
1684 . 
1605 . 

1607 . 

1608 . 

1609 . 

1610 . 
1611 . 
1612 . 
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1613 . 

1614 . 

1615 . 
1617 . 

1616 . 
1619 . 
1626 . 
1621 . 
1629 . 


2099 . 

2100 . 
2102 . 

2103 . 

2104 . 

2105 . 

2106 . 
2107 . 
2103 . 
2120 . 


5001 PRIHTPRINTPRINT 

8000 PR I NTCHR* < 1 47) ! PR INT : PR I NT : PR I NT : PR I NT 
3001 PRINT"REMflRKS FOR PROGRAM RESP-96D" 

8002 H*='"':PRINT:PRINT:PRINT:iHPUT"HRRD COPV IF V 

8003 IF Hf=''V"THEN3050 
3010 LIST 

8015 END 

3050 print: PRINT 

8052 PRIHT'TLERSE CLOSE PRINTER WHEN LIST IS FINISHED" 

8035 OPEN4>4,0!CMD4 

8060 PRINT REMARKS FOR PROGRRM RESP-96D 

8062 print: PRINT" 10/'27/35" 

8070 PRINT :PRINT"THE FOLLOWING REMRRKS ARE BV LINE NUMBERS WHICH" 
8072 PR I NT "CORRESPOND TO THOSE IN THE PROGRRM." 

3075 print: print: PRINT 
8080 LIST 
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APPENDIX D2 


COMPUTER PROGRAM WAT-CA96D 
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2 REM ** WflT-CH96F m REV. 12/09/35 

3 PRINTCHR»<147)'' INITIALIZING FILE AREA PLEASE WAIT" 

5 REM CHANCED TO PERMIT DATA TABLES OF 99 ITEMS (WAVELENGTHS) 

6 REM CHANGED TO PERMIT QUICK CHANGE OF CONSTANTS 

7 REM BETTER HEADER EDITING ADDED 

10 DIM DT(8<10O)iAfC8.3)iEi:<2^3)/CCi:2;35) 

50 ZZfo"" 

100 PRINTCHR$<147):PRINT 

101 PRINT" PROBLEM 77 - TO PREVENT LOSS OF DATA RESTART WITH GOTO100" 

102 PRINT!PRINT:PRINTZZ|::PRINT 


103 PRINTTAB<8) "•*****» MAIN MENU : PRINT: PRINT 

104 PRINTTAB<6)"1 - CALC. SVSTEM OUTPUT" 

105 PRINTTAB<6)"2 - CALC. COMPOSIT RESP. " 

106 PRINTTAB<6)“3 - EDIT" 

107 PRINTTAB<6)"4 - SAVE TO DISK" 

108 PRINTTAB<6)"5 - INPUT FROM DISK" 

109 PRINTTAB<6)"6 - KEVBOARD INPUT COMPONENT DATA" 

110 PRINTTAB<6)"7 - PRINT TABLE DATA" 

111 PRINTTAB<6)"8 - DISK DIRECTORV" 

113 PRINTTAB<6)"11 - SVSTEM OUTPUT TIMES CONSTANT" 

114 PRINT:PRINT!PRINTTABa5)J 

115 INPUT"COMMAND "iK! IFK<1THEN100 

116 IFK>11THEN100 

118 ONKGOTO4O0i200n0OOi800il2OO<60Oi 1400; 130. 100. 1500.2400 

130 OOSUB2O06 

131 INPUT"CONT. ON C.R. "IX 

132 GOTO 100 

200 PRINTCHRf<147) 


201 PRINT :PRINT"*»» CALCULATE COMPOSIT RESPONSE m*." 

202 PRINT 

204 print: PRINT 

209 D 1 B0 : D2«0 : D3»0 : D4»0 : DS^O : D6-0 

210 Kfa""iINPUT"USE DETECTOR RESPONSE IF V PRINT 

211 IFK5»"V"THEND1=1:D2=1 

214 K#a""!INPUT"USE LIGHT PIPE RESP IF V ".’KJfjPRINT 

215 IFK»*"V"THEND3»1 

218 UPPER DET. TRIM IF V ";Kf:PRINT 

219 IFKJa"V"THEND5»l 

222 KJa""!INPUT"USE LOWER DET. TRIM IF V "jK*: PRINT 

223 IFK*="V"THEND6al 

226 SVSTEM RESPONSE IF V "iK* 

227 IFK$si"V"THEND4“l 

228 PRINT"Dls"Dl" D2=»"D2" D3="D3" D4 b«D4" D5="D5" D6="D6 

240 PRINT:K$«""!INPUT"MAKE CHANGES IF V "iK* 

241 IFK*a"V"THEN200 

242 PRINT! print: PRINT SPRINT 

252 PRINT"* CALCULATING COMPOSIT RESPONSE CURVES *" 

255 A*C7.4)="D1«»"+STR*CD1)+" D3n"+STR*<D3)+" D4a"+STR*<D4)+" D5="+STR$(D5) 

256 R*<8.4)a"D2="+STR*CD2)+" D3="+STR*<D3)+" D4="+STR*(D4>+" D6="+STR*(D6) 

257 G0T0312 

260 PRIHTCHR*<147) "PRESENT FILE NAME IS "R*<K.2) 

264 print: PRINT 

266 F0RI=1T03:PRINTI" "A*<K. I ) : NEXT 
263 L=0:INPUT"CHAHGE LINE "iL 

269 IFL=0THENRETURN 

270 PRINTA*<K.L) 

271 Kfa"": INPUT" CHANGE TO "iK* 

272 IFK*>""THEHA*<K.L)=KS :GOTO260 
274 RETURN 

280 PRINTCHR*<147) "PRESENT FILE NAME IS "B*<0.2> 

284 PRINT (PRINT 

286 F0RlalT03:PRINTI" "B*<0. 1) .'NEXT 
233 L=0:INPUT"CHANGE LINE "iL 

289 IFL=0THENRETURN 

290 PRINTB*<0.L) 

291 INPUT" CHANGE TO "iK* 

292 IFK*>""THENB*C0.L)=K* :G0T0280 
294 RETURN 



t-z 


312 F0RIalT096 

313 DT<7»I)b1 

314 DT<8,I)=1 

315 NEXT 

320 F0RI-1T096 

323 IFDl»lTHENDT<7iI)=DT<7^n*DTa/n 
325 IFD3*lTHENDT<7i I)=DT(:7. 1)*DT<3^ 1) 

327 IFD4=1THENDTC7, 1)»DT<7^ D*DT«4/ 1) 

329 IFD5=lTHENDT<7i I)»DT<7/ 1)#DT<5 j I) 

330 PRINTI" "DT<7,D 

331 FORL=1TO100!NEXT 

332 NEXT 

333 PRINT"UPPER DET. CflLC. COMPLETE" 

340 F0RI»1T096 

343 IFD2=1THENDT<8/ D=DT<8^ I)*DTC2. 1) 

345 IFD3=1THENDT<8, D»DT<8j I)*DT<3/ I? 

347 IFD4»1THENDT<8; I)=DTC8i I)*DT(4. D 
349 IFD6»1THENDT<8/ I)»DTC8< I)#DT<6i I) 
'350 PRINTI" "DT<8/I) 

351 FORL»»1TO100:HEXT 

352 NEXT 

353 PRINT"LOWER DET. CRLC. COMPLETE" 

355 FORLelTO1000:NEXT 

360 22*=>"C0MP0SIT RESP. CURVES COMPLETE" 
370 GOTO I 00 


400 REM 

401 PRINTCHR$C147) 

402 PRINT! PRINT"** CflLCULRTE SYSTEM OUTPUT CURVE **" 

403 PRINT !PRINT"DEFflULT VALUES:" 

410 PRINT"STRRTING TEMP. = 340 DEG. C" 

411 PRINT" INTERVAL » 20 DEG. C (MINIMUM VHL. )" 

412 PRINT :PRIHT"ENDING TEMP = 2000 DEG. C" 

413 PRINTiPRINT 

416 TS=340!INPUT"ENTER STARTING TEMP. ";TS 

417 PRINT iTD=»20!INPUT"ENTER INCREMENT" JTD 

418 IFTD<:20THENPRINT"USE VALUE GREATER THAN 20 ":G0T0417 

419 PRINT !TE»2000!INPUT"ENTER ENDING TEMP. ";TE 

420 PRINTiX7=7.114E-07:PRINT"TOP DET CONSTANT = "X7 

421 INPUT"CHANGE TO "iX7 

422 PRINT!X8»1.4423E-07:PRINT"BOT DET CONSTANT = "X8 

423 INPUT"CHANOE TO "iX8 

424 PRINT :PRINT"W93. 30 " :W9=.30: INPUT"CHANGE WAVELENGTH TO ";W9 

425 PRINT :PRINT"K9=1 " :K3»1 : INPUT"CHANGE CONSTANT MULT FACTOR TO ";K3 

428 L=0:M»1 

429 B$<0i3)«’"84" ' B»<0i4 J=STRS (TS) : 5)»STRJ <TD) : Bf (0/ 6)=STRT <TE) 

430 R=1 

431 PRINT!L$e"":iNPUT"ENTER V FOR PRINTOUT 'SLf 

432 X9=0!IFL#="V"THENX9a4 

434 IFX9*4THENOPEN4,4,0 

435 PRINT"CREATING TABLE - PLEASE BE PATIENT" 

436 Cla37418 

437 C2-14337.9 

433 GOSUB560 

439 IFX934THENGOSUB502 

440 T«TS 

441 K»T+273.15 

442 REM WR-TOP DET WB=BOT DET. H9=STRRTING WAVELENGTH 

443 REN ORIGINAL CONSTANT IN LINE 452 WAS 7.114E-07 

444 REN CRLC AMPS FOR EACH 0. I MICRION WAVELENGTH INCREMENT 

445 WT=0!WU=0 

446 S7«.01*X7*K9!S3=.01*XS*K9 

447 REM SPEED CRLC COMBINE CONSTANTS 

448 F0RD=1T096 

449 V=W9+<D-1)*.01 

450 WXa<ClXVt5>/<<EXP<C2/'(V*K) ) )-l > 

452 WR=WX#DT<7,D)*S7 

453 WBeWX*DT(8<D)*SS 89 

454 WTaWT+WR 

455 WU«WU+WB 

456 NEXT 



453 CC<0,R)=WT 
460 CCaiR)=WIJ 
467 IFXS=4THEHGOSUB550 
463 GOSUB570 

474 IFT=TETHBM82 

475 IFTMETHEN432 
473 T*T+TD 

479 R=R+1 

480 G0T0441 

482 Z2*='' OUTPUT CURVES COMPLETE" 

484 IFX9=4THENPRINT#4!CL0SE4 

485 INPUT"CONTINUE ON C.R. ";WIJ:GOTO100!REM ALT TO LINE 484 
490 IFI=0THENJ=7 

492 IFI*1THENJ*8 
494 RETURN 


500 REM 

502 PRINT#4< "SYSTEM OUTPUT CURVES" :PRINT#4 
504 PRINT :PRINT"PREVIOUS FILE NAME® "BJCL.2) 

506 K*="":iNPUT"ENTER NEW FILE NAME ";Kf 

507 INPUT" IF NOT OK ENTER N ";LI 

503 IFLf»"N"THEN506 

509 IFKJ=""THENKf=Bf<:0,2) 

510 B#<0i2)=Kf :Bf<0,l)="0CC" 

511 PRir'IT#4,TRB<5) "SOURCE FILES= "A$<:7,2)" "A$<3/2) :PRIHT 

512 PRINT#4,"2 FILE NAME: "Bf C0. 2) : PRIHTM 

513 PRIHT#4<" COLUMN "VAL<Bf<0/ 1)>" "ETCG/S)" LINES START AT "ET(0.4 

514 PRINT#4i" DELTA T "Bf<0i5)" STOP AT "Ef(0.6> 

515 PRINT#4i"7 "BJ<0/7) :PRINT#4/ "8 "E$(0i8) 

516 PRINT#4 

517 PRINT#4,"T0P BET CONSTANT <X7)= "X7:PRINT#4> "EOT DET CONSTANT <X3>= "X3 

518 PRINT#4i "STARTING WAVELENGTH® "W9" SYSTEM CONSTANT K9="K9 

519 PRINT#4 

520 PRINTtt4iCHR*<27)CHR$a6>CHRJ<0)CHRJ<000) "TEMPERATURE 

521 PRINT#4i CHR$<27)CHR$a6)CHRJ <0)CHRTC 100) "PHOTOCURRENT" ; 

522 PRINT#4,CHRf<27)CHRf<16)CHR$<0)CHR3:':230) "PHOTOCURRENT" 

523 PRINT#4>CHR*<27)CHR*C16)CHRfC0)CHRTC000) "DEGREES C "; 

524 PRINT#4,CHR|:<27)CHR$<16>CHR#C0)CHRf<100)"RMPERES * TOP"; 

525 PRINT»4,CHR$C27)CHR*a6)CHR$<0)CHRT<230)"RMPERES * EOT"; 

526 PRINT#4> CHRJ:<27)CHRf (16)CHRf ( 1 JCHRSC 100) "RATIO" 

527 PRINT#4 
529 RETURN 

550 PRINT#4iCHR*<27)CHRJ:<16)CHRf<0)CHRJ<000)T; 

551 PRINT#4iCHR*<27)CHRfC16)CHR*<0)CHRf<100)CC<0iR); 

552 PRINT#4, CHR»<27)CHRT< 16)CHR$<0)CHRT<230)CC< 1 , R) ; 

553 PRINT#4, CHR*<27)CHRJC16)CHRT< 1 )CHRf < 100)CC(0i R)/CC< 1 . R) 


554 

RETURN 




560 

561 

562 
564 

PRINT" 

RETURN 

PRINT 

RETURN 

T TOP I 

BOT I 

RATIO" 

570 

572 

574 

PRINTT" 

PRINT 

RETURN 

"CC<0.R>" 

"CC<1,R)" 

"CC':0.R)7CC<1.R) 


600 PRINTCHR$<147) 

602 PRINTTAB<3) "*»•***** DATA TABLE MENUE #»*»**#*" 

604 print: PRINT 

606 PRINTTAB<5)"1 - UPPER DETECTOR RESPONSE" 

607 PRINTTAB<5)"2 - LOWER DETECTOR RESPONSE" 

605 PRINTTAB<5)"3 - LIGHT PIPE RESPONSE" 

609 PRINTTAB<5)"4 - SYSTEM RESPONSE" 

610 PRINTTABC5)"5 - TRIM CURVE - UPPER DET." 

611 PRINTTAB(5)"6 - TRIM CURVE - LOWER DET. " 

612 PRINTTAB<5)"7 - INPUT FROM DISK": PRINT 

613 PRINTTAB(5)"8 - EDIT DATA" 

615 PRINT 

616 PRINTTRB<5)"9 - MAIN MENU" 

624 PRINT 

626 PRINTTAB<5)"10 - SAVE DATA ON DISK" 

628 print: PRINT :K=0 

630 PRINTTAB<12); :iNPUT"COMMAND ";K 



631 IFK=0THEN600 

632 IFKM0THEN600 

633 PRINTCHR* < 147> : PRINT : PRINT 

634 0NKG0T0651^633/6S5/657/653;66l. 1200, 1000, 100,800 

651 PRINT”ENTER RESPONSE OF UPPER DETECTOR" :fl$<K,4)="T0P DET RESP" 

652 fl$<l,l>="l":COTO700 

653 PRINT"ENTER RESPONSE OF LOWER DETECTOR" •flf<:K,4)*"iiOT DET RESP" 

654 flf <2, 1 )="2" ; GOTO700 

655 PRINT"ENTER TRANSfllSSION OF LIGHT PIPE" :fl|:<K,4)="FIEER TRANS" 

656 Rt<3,l)="3":GOTO700 

657 'PRlHr’ENTER SVSTei RESPONSE" : 4 5* "SYSTEM R'ESP" ' 

650 fi$<4, 1)*"4" !GOTO700 

653 PR1NT"ENTER UPPER DET. TRIM CURVE" :flf<K,4)="T0P DET TRIM" 

660 n$<5,l)="5":GOTO700 

661 PRINT"ENTER LOWER DET. TRIM CURVE" •■flf<K,4)="E0T DET TRIM" 

662 flf <6, 1 )="6" : GOTO700 
700 print: PRINT 

710 1*1 

712 PRINT" I*" I" "DT<K/D" : V*0: INPUT"NEW VfiLUE = ";V 

714 IP/>0THENDT<K, D*V 

716 IF V ";Lf : IFL**"V"THEN730 

718 PRINT 

720 1*1+1 

721 IFI>36THEN100 

722 G0T0712 

730 Q*0:INPUT"RESTRRT RT I* ",Q 

731 IFQ=0THEN600 

732 IFQ>0THENI=Q 
734 G0T0712 

733 GOTO600 


800 PRINTCHR$< 147) 

801 PRINT"*#***#** SRVE DATA ROUTINE ******#*« 

802 PRINT : PRINT : PRINT : PRINT : GOSUB2000 

803 PRINTCHR$a47) 

804 PRINTCHR*< 147) "*********** SAVE DRTR CURVE #**#**#****" 

808 PRINT 

312 PRINTTRBC3)"! - UPPER DET. RESP. " 

313 PRINTTRB<3)"2 - LOWER DET. RESP." 

314 PRINTTRB<3)"3 - LIGHT PIPE RESP." 

315 PRINTTRB<3)"4 - SYSTEM RESPONSE " 

816 PRINTTRB<3)"5 - UPPER TRIM CURVE" 

317 PRINTTRB<3)"6 - LOWER TRIM CURVE" 

313 PRINT 

324 PRINT:PRINTTRB<4)"3 - RETURN TO MAIN MENU" 

826 PRINT 

827 print: PRINTTAB<4)" 11 - RETURN TO DATA TABLE MENU" 

823 PRINT :PRINTTAB<4>" 12 - SAVE SYSTEM OUTPUT CURVES" 

823 print: PRINT 

330 K*fl:iNPUT"SAVE CURVE NO. ",K 

832 IFK*0THEN800 

833 IFKM2THEN300 

834 ONKGOTO836,836,836,33S,836,836,300,300, 100, 1510,600,500 

836 GOSUB260 

837 PRINT :PRINT"CURRENT FILE WITH ABOVE NAME WILL BE ERASED FROM DISK" 

842 PRINT:L*="":INPIJT"MENU IF Y "iLf 

843 IFL*="Y"THEN600 

844 Af<K,l)*STRf(K):RJ:<K,3)="3S" 

345 0PEN15, 8, 15, "S : "+A$<K, 2) = INPUT# 15, A$, B# , Cf , D$ : PRINTED 

346 CL0SE15 

843 OPEN2, 3, 2, "0 : "+Af <K, 2)+" , S, W" 

850 F0RI=1T03 

351 Af<K, I)*"A"+A$<K, 1) :PRINT#2,Ri:<K, I) 

852 NEXT 

360 F0RI*1T036 

862 PRINT#2,DT<K,1) 

864 NEXT 

866 PRINT#2, "END" 91 

363 CLOSES 
870 GOTO 100 



S00 PRINTCHRf<147) 

904 print: PRINT! print: PRINT 
906 PRINT"*#*#* SAVE - SYSTEM OUT 


CURVES ***.*#*" 


908 print: PRINT 
910 PRINTTflB(4)"l 

912 PRINTTflB<4)"3 

913 PRINTTRB<4)"4 
'914 PRINTTRB<4)"7 

915 PRINTTRB<4)"8 

916 PRINTTRB<4)"9 
929 print: PRINT 


*# SYSTEM OUTPUT **" 

M 

COMPONENT RESPONSE CURVES" 
COMPOSIT RESP. - UPPER BET 
COMPOSIT RESP. - LOWER BET 
MAIN MENU" 


932 K= 0 ! INPUT"COM M RNB" ; K 


II 


tl 


: PRINT 


933 IFK=0THEH900 

934 IFK>10THEN90O 

935 GOSUB2000 

936 ONKOOTO973,900/900..800/900/900^933;933^ 100. 1520 
933 PRINT"CURRENT FILE NRME IS-' "flf<K.2) 

939 G0SUB260 

940 PRINT :PRIHT"EXISTING CURVE OF THE SAME NAME WILL BE ERASEB" 


544 C*="":iNPUT"MENU IF Y "JCf 


945 IFC#="Y"THEN900 

946 R$<K, 1)=STR«:<K) :fl|:<K.3)*"96" 

947 0PEN15/3.15."S:"+Rl:<K.2) 

943 INPUT#15.R».Bf.Cf :PRINT" "B$:CL0SE15:L$="" : INPUT"MENU IF Y".:LT 

949 IFL$="Y"THEN900 

950 0PEN2. 8.2. "0: "+R^<K. 2)+" . S, W" 

954 F0RI=1T03 

955 Cfa"R"+R3:<K. I > : PRINT#2. C* : PRINTC* 

956 NEXT 

960 F0RI=1T096 

962 PRINT#2. DT<K. I ) : PRINTBT<K. I ) 

966 NEXT 

970 PRINT#2. "ENB" 

972 CL0SE2 

973 GOTO900 

974 GOSUB230 

975 PRINT :PRINT"EXISTING CURVE OF THE SAME NAME WILL EE ERASEB" 

978 PRINT"FILE NRME NOW IS "Bf<0.2> 

980 C«*>"":INPUT"MENU IF V ";C$ 

982 IFC*n"V"THEN900 

983 Bf<0.1)="0":B$<0.3)a"34" 

934 OPEN15.8.15."S:"+Bf<0.2) 

985 INPUT#15.fl:f.B$.Ci:.II*:pRINTBf 

986 CLOSE 15 

983 0PEN2. 3. 2. "0 : "+Bf <0. 2)+" . S. W" 

990 F0RI=1T03 

991 CJ:="R"+B*<0. 1) :PRIHT#2.Cf :PRINTC* 

992 NEXT 

993 F0RI=1T034 

994 PRINT#2.CC<0. n :PRINTCC<0. D 

995 PRINT#2. CC< 1.1): PRINTCCC 1 . I ) 

996 NEXT 

997 PRINT#2. "END" 

993 CL0SE2 

999 GOTO 100 


1000 REM 

1010 PRINTCHRf<147) 

1012 PR I NT "***********»** EDIT MENU *****************" 
1014 print: print: PRINT 

1016 PRINT"! - UPPER DET. RESP." 

1017 PRINT"2 - LOWER DET. RESP." 

1018 PRINT"3 - LIGHT PIPE" 

1019 PR I NT "4 - SYSTEM RESP." 

1020 FRINT"5 - UPPER TRIM CURVE" 

1021 PRINT"6 - LOWER TRIM CURVE" 

1022 PRINT"7 - *# SYSTEM OUTPUT **" 

1023 PRINT"8 - *** EXAMINE FILE ***" 

1024 PRINT"9 - MAIN MENU":PRINT 

1025 PRINT" 10 - COMPOSIT RESP. * UPPER" 

1028 PRINT"11 - COMPOSIT RESP. * LOWER" : PRINT: PRINT 

1029 K=0 : INPUT"COMMRND" J K 



1030 IFK-8THEN31O0 

1032 IFK-0THEN1000 

1033 1FK>11THEN1000 

1034 IFK=9THEN100 

1035 IFK<7THEND*96 

1036 IFK*7THEMD*84 = GOTOl 067 

1037 IFI09THENIFK<12THEHD»96 

1038 1FK*10THENK=7 

1039 IFK=11THENK*8 

1042 IFDS84THEN1067 

1043 PRINTCHR*<147) 

1044 PRINT 

1045 PRINT'TILE: "R*<K,2) 

1046 F0RI=lT08:PRINTr' I) :NEXT 

1047 P»0!lNPUT“CHflNCE LINE NO, "JP 

1048 PRINTfl*<K,P) 

1049 IFP»0THEN1052 

1050 lNPUT"CHnNGE TO ";fl*<K,P> 

1051 GOTO1044 

1052 E»0:INPUT"STnRT AT ITEM ";E 

1053 FORIsETOE+lQ 

1054 IFI>96THEN1060 

1055 PRINTI" "DT<K,I> 

1058 NEXT 

I860 H»8:iHPUT"CHRNGE LINE NO. "JH: IFH*0THEN1062 

1061 IHPUT-CHRNGE TO “;DT(K.H) 

1062 IFH<1THEN1064 

1063 GOTO1053 

1064 E=E+1S 

1065 IFI>96THEN1000 

1066 GOTOl 053 


1067 PRINTCHR»<147) 

1068 IFK*7THENK*0:PRIHT«K="K 

1069 F0RI»IT08: PRINTI" I) :HEXT 

1071 C=0:INPUT"CHRNGE LINE NO. "JC! IFC=0THEN1075 

1072 PRINT B*<0.C) 

1073 INPUT"ENTER CHANGE 

1074 GOTO1069 

1075 C*»"":1HPUT"HENU IF V ";C#: IFCf="V"THEN10lO 

1076 PRINT : PRINT : E=1 : INPUT"STflRT AT ITEM NO. ";E 

1078 PRINT 

1079 PRINT"ITEM UPPER DET LOWER DET":PRIHT" RATIO" 

1080 F0RI=ET0E+9 
1032 IFI>84THEN1096 

1034 PRINTI ; 

1035 PRINTTRB<4)CCC0,I>; 

1036 PRINTTRB<4)CC(lj D^IFCCCIj I>=0THEHPRINTTflE<15)"#****:+:^^T*" 

1037 PRINT THB<15)CC<K,I>/CC<R+1,I5 
1088 IF1>03THEH1O96 

1039 NEXT 

1096 D=0:1NPUT"CHBNGE LINE NO. ;D: IFD=0THENU 10 

1097 IHPUT"CHRNGE UPPER VALUE TO " J CC<0.D>RDV. 

1099 INPUT"CHRNGE LOWER VRLUE TO "JCCCI/D) 

1105 PRIHTCHR$a47):G0T01G79 

1110 IFI>84THEN1000 
nil E=E+10 

1112 C*="": INPUT"NEHU IF V ";Cf 

1113 IFCf="Y"THEN1000 

1121 FRlNTCHR*a47):PRIHT:PRIHT 

1122 GOTO 1070 


1200 PRlHTCHRJt l47)"**Sf****** INPUT CURVE FROM DISK " 

1203 PRINT :PPINT"TO EXAMINE DISK FILE CO TO EDIT FUNCTION OR * G0TQ3 1 00 PRINT 

1204 print: PRINT 

1206 PRINT''DIRECTORV: " = PRINT 

1208 COSUB200O 

1209 PRINT 

1210 K»="":iNFUT"ENTER FILE NAME ";K$ 

1220 E*="": INPUT"CORRECTIClN IF V "JET 93 

1221 IFET="V"THEH1224 

1223 GOTO 1230 

1224 IMPUT"EHTER CORRECTION "JKT 

1230 ET="": INPUT-MENU IF V "JET 

1231 IFET*"V"THEN100 



1232 0PEH15,8,1S 

1234 0PEH2. 8, 2, "0 : ••+KJ:+" , S. R” 

1235 F0RI=1T08 

1236 IHPUT#2.C»<I):PRINTr' ■■MID|:<Cf<I).2> 

1237 NEXT 

1238 Dl=VflL<MID$<C$<l),2>) 

1239 D3»VRL<MII3»<C$<3),2>) 

1240 GOTO1310 

1250 0PEN2i 8, 2, "0 : ’•+K*+" / S; R" 

1251 IFD1=0THEN1271 
1254 IFD3=34THEN1271 

1258 FOR I = 1 T08 : 1 NPUT#2 , C$ : fl$ < D 1 , D =M I DJ C C$ . 2 : NEXT 
1260 F0RI=1T096 

1262 INPUT#2iDT<Dlin 

1263 PRINTI" ''DT<D1,D 

1267 NEXT 

1268 CL0SE2 

1270 GOTO 100 

1271 FORI=1TO8:IHPUT#2,C$:BI:<0, n=MHJf<W.2) :NEXT 

1272 F0RI=1T084 

1274 INPUT#2,CC<0i I) : INPUT#2.CC<1 / I) 

1279 NEXT 
1281 CL0SE2 
1284 GOTO 100 

1310 lHPUT#15,flf,Bf/Cf.Df 
1312 CL0SE2:CL0SE15 
1314 IFflf="00“THEN1320 

1316 PRINT : PRINTB$ :FORI=1TO20O0 : NEXT : GOTO1200 

1320 REN 

1321 Ef*"":iNPUT"MENU IF V ";E3:: IFEf="V"THENGOTO1200 

1329 I)l=VflL<riIB*<W<l),2)) 

1330 PRINT !PRIHT"DESTINfiTION IS FILE NO. "D1 

1332 PRINT :PRINT"FILES RRE flS FOLLOWS:" 

1333 PRINTTflB<8)"0 - SYSTEM OUTPUT m" 

1334 PRINTTRB<3)"1 - TOP BET RESP" 

1335 PRINTTRB<3)"2 - BOT BET RESP" 

1336 PRINTTRB<3)"3 - LIGHT PIPE RESP" 

1337 PRINTTRB<3)"4 - SYSTEM RESPONSE" 

1338 PRINT 

1339 PRINTTRB<8)"7 - COMPOSIT RESPONSE - TOP" 

1340 PRINTTRB<8)"8 - COMPOSIT RESPONSE - BOT" 

1344 PRINT :PRINT"FILES GENERRTED WITHIN THIS PROGRAM HAVE THE PROPER NUMBER. 
1346 PRINT : PRINT :1NPUT"CHRHCE DESTINRTION TO ";D1 
1348 L$="":INPUT"MENU IF V ";Lf : IFLf="V"THEN1200 
1350 GOTO I 250 


1400 PRINTCHR*<147) 

1402 PRINT"**********'*** PRINT DRTR *************" 
1404 PRIHT:PRINT:PRINT 

1406 PRINTTRB<6)"1 - UPPER BET. RESP." 

1407 PRINTTRB<6)"2 - LOWER BET. RESP." 

1408 PRINTTRB<6>"3 - LIGHT PIPE RESP." 

1409 PRINTTRB<6)"4 - SYSTEM P£:SP. " 

1410 PRINTTRB<6)"5 - UPPER TRIM" 

1411 PRINTTRB<6)"6 - LOWER TRIM" 

1412 PRINTTRB<6)"7 - COMPOSIT RESP. - UPPER" ' 

1413 PRINTTRB<6)"8 - COMPOSIT RESP. - LOWER" 

1415 PRINTTRB<6)“9 - MENU" 

1416 PRINT 

1417 PRIHTTRB<6)"10 - ** SYSTEM OUTPUT **" 

1420 INPUT"COMMRNB ";K 

1422 IFK=0THEN1400 

1424 '1FK>10THEN1400 

1426 IFK*9THEN100 

1428 IFK<9THEN1434 

1430 IFK=10THENK=0: GOTO 1460 

1434 OPEN4,4.0:CMB4 

1436 PRINT"FILE: "R$<K.2) :PRINT 

1440 PRINT"COL "K; :PRINT" "flf<K.3V BRTR ITEMS" 

1441 FORI=4TO0: PRINTI" I) :NEXT 

1442 GOTO1700 

1443 PRINT :PRINT"WfiVELENGTH "RD* 



1444 F0RI=1T096 
1446 V=.3+<I+1)*.02 

1443 PRINTVCHRf<16)CHR$(43)CHR$(:53?DT(K, I) 

1450 NEXT 

1452 PRIHT#4:CL0SE4 
1456 GOTO1400 

1460 OPEH4,4>0:CMD4:PRIHT''SYSTEH OUTPUT CURVES" : PR I NT 

1461 PRINT" SOURCE CURVES* "fif<7^2)“ S: "fl$<8..2> = PRINT 

1462 P0RI=lT08:PRINTr' "BJ(0. DNIEXT 

1463 PRINT!TS=VRL<B$<0.4))!TD=VfiL<Bf<O,5))!TE*VflL<Bf<0..6>> 

1464 PRINT"STflRT RT "B^:<K.4>’S DELTR T "B*<K.5)".. END RT "EJCK.6) 


1465 PRINT 




1466 PRINT" 

UPPER DET. 

LONER DET. " 


1468 PRINT"TEMP C. 
1470 PRINT 
1472 F0RI=1T084 
1474 T*TS+<TD*<I-1)) 
1476 PRIHTTJ 

CURREN 

CURRENT 

RATIO" 


1478 PR I NTCHRf <16) CHR$ <49)CHR3:<49)CC<K.I); 

1480 PRINTCHR$<16)CHR|:<51)CHRT<49)CCCK+1, D; 

1431 IFCCa>I)*0THENCC<l.l)=l 

1482 PRINTCHRf<16)CHR$<52)CHRf<45)CC<K/ I)/CC<K+1.. I) 
1433 IFT=TETHEN1483 

1484 IFT>TETHEN1483 

1485 NEXT 

1438 PRINT:PRINT 
1490 PRINT#4:CL0SE4 
1492 GOTO1400 


1700 0NKG0T0171121.. 1712. 1713.. 1715. 1716. 1717. 1713. 1719 

1711 RDJ="UPPER DET RESP. RETURN 

1712 RDf="LOWER DET RESP. " .RETURN 

1713 RD$="LIGHT PIPE TRRHS" :RETURN 

1714 RDf="SVSTEM TRRNS" = RETURN 

1715 RD$="UPPER TRIM TRRHS" = RETURN 

1716 RDf="LOWER TRIM TRRNS" : RETURN' 

1718 RDf="COMP. RESP. UPPER" = RETURN 

1719 RDf="COMP. RESP. LONER" : RETURN 
1730 RETURN 


2000 H$="":INPUT"LIST DIRECTOPV IF V ";Hf 
2002 IFHf*"V"THEN2006 
2004 RETURN 

2006 LS="": INPUT" IF PRINTOUT NRNTED PRESS V ";Lf 

2007 KK=l:REt1 PRINT COUNT 

2008 OPEN1.8.0."«0" 

2009 IFLf="V"THEH2015 

2010 COTO2020 
2015 OPEN4.4.0 
2020 GET#l.R*.Bf 
2030 GET#1.RJ.B$ 

2040 GETttl.RT.Bf 
2050 C=0 

2060 IFRfO""THENC=flSC<Rf) 

2070 IFBJO" "THENC=C+RSC < E$ ) *256 
2030 PRINTriID#<STRf <C) . 2) ; TRB<3) ; 

2081 lFLf*"V"THENPRINT#4.riIDi:<STR3:<C).2); 

.._2S230..QEX»liBf lIESK>ilTJ:lEhl2200 1. — 1 

2100 IFB$OCHR$<34)THEN2090 

2110 GETttl.B*: IFBfOCHR$<.34)THEHPRINTBf.: 

2111 IFL*="V"THEN2114 

2112 IFBfOCHR#<34)THEN2110 

2113 GOTO2120 

21 14 1FBK>CHR*<34>THENPRINT#4. CHRf < 16>CHRf <000)CHRT<70>Bf ; ' G0T021 10 
2120 GET#l.Bf : lFB*=CHRf<32)THEN2120 

2130 PR1NTTRB<18); -CT*"" 

2140 Cf =C*+B$ : GET# 1 . Bf : I FBf O " " THEN2 140 

2150 PRINTLEFT5<Cf.3):KK=KK+l 

2151 IFLJ="V"THENPRINT#4.CHRf<16)CHRi:(:050)CHRT<:000)LEFTf<:CT..3:. 95 

2160 GETTf ; IFTfO""THENGOSUB2300 

2 1 64 I FKK*20THENKK= I : G0T02 1 63 
2166 GOTO2170 

2168 INPUT'TONTINUE ON C.R. ",:K 
2170 IFST=0THEN2030 \ 



2200 PRINT" BLOCKS FREE" 

2201 IFLJ="V"THENPRINT#4/ " BLOCKS FREE" 

2205 CLOSE 1 

2206 IFLf="Y"THENPRIHT#4 : CL0SE4 
2210 RETURN 

2300 1FTJ="Q"THENCL0SE1 =00102205 
2310 OETT$=IFT#s"“THEN2300 
2320 RETURN 


2400 PRINTCHRfa47)"f***#!^** SVOUTifCONSTflNT 
2402 PRINT=PRINT 

2404 PRINT"DEFflULT CONSTANT IS 1.0 
2406 KT=1=KB=1 

2408 INPUT "ENTER TOP CONST. ";KT 
2410 1NPUT"ENTER BOT CONST. "JKB 
2412 flf<0.7)="KT*"+STRf<KT)+" KB="+STR$<KB> 
2420 F0RI=1T084 

2422 CC<0, 1)=CC<0/Di*!KT:CC<l, I>=CC<1, I).T:KB 

2423 PRINTCC<0iI>" "CCCl,!) 

2424 NEXT 

2425 PRINT"FINISHED" 

2426 FOR I =» 1 TO800 = NEXT 
2460 GOTO 100 

3000 0PEN15,8, 15 

3010 INPUT#15,fl*,B»iC*;D* 

3020 CLOSE 15 

3030 PRINTfl*,B».Cf,DJ= 

3040 RETURN 


3100 PRINTCHR$<147?"#*** EXOrilNE FILE ****4i*»:m:**r***" = PRINTS 
3110 INPUT"FILE NflriE";F$ 

3120 INPUT"FILE TVPE";TS= 

3130 T*=LEFT$<Tf,l> 

3140 IFT5O"S"THENIFT*O"P"THENIFT$O"U"THEN3120 
3145 0PEN15,3,15 

3150 OPENS. 8, 5, "0= "+Ff+". "+T$+"R" 

3160 GOSUB3200 

3170 GET#5,nJ 

3180 IFST=0THEN3190 

3185 IFST=64THENCL0SE5. 15=GOTO3290 

3187 PRINTST=GOTO3290 

3190 PRINTR*; 

3191 GOTO3170 

3200 INPUT#15iR*,BS,CJ.DJ 

3210 IFVnL<R$»0THEN PRINTHf.Bf.Cf .Dl= 

3220 RETURN 

3290 1NPUT"RETURN ON C.R.";X:f 

3291 GOTO 1000 
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1 GOTO 8000 REM fiUTO LIST INSTRUCTION 

2 * WflT-Cfi96F-RD m 

3 INITIRLI2ING DONE BV RUN COMMAND. 

5 DIM STATEMENTS: 

6 . DT # 8 MAIN DATA TABLES. 

7 . AT * HERDER LINES FOR ABOVE DATA TABLES. 

10 . CC * SVSTEM OUTPUT CURVES. 

12 . * HERDER LINES FOR OUTPUT CURVES. 

50 . 22 IS A FLAG # SEE LINES 102 AND 201. 

100 CHR«:a47) IS CLEAR SCREEN CMD. 

101 itiIMPORTRNT!* RESTART AT MENU TO PRESERVE DATA. 

102 PRINTS COMPOSIT CURVE COMPLETE * IF TRUE. 


104 . 

105 . 
108 . 
107 . 


103 . 

109 . 

110 . 

113 . 

114 POSITION COMMAND REQUEST. 

115 . 

113 ADDRESSES OF MENU SELECTIONS. 


130 LINE 2006 IS DIRECTORV PRINTOUT SUBROUTINE. 

131 END OF DIRECTORV LISTING RETAINS SCREEN DATA. 

132 . 


200 . 


201 . 

202 . 

204 . 

209 IHITIRLI2E CURVE SELECTION FLAGS. 

210 2ERO IS "DONT USE" ONE IS "USE". 

211 . 

212 . 

214 . 

215 . 

218 . 

219 . 

222 . 

223 . 

226 . 

227 . 

228 VERIFICATION OF SELECTIONS. 

240 PERMITS START OVER IN CASE OF ERROR. 

241 . 

242 . 

252 . 

255 LOADS ITEMS 7 « 8 OF HERDER WITH LIST OF CURVES USED. 

256 . 

257 LINE 312 IS START OF CALCULATIONS. 

260 START OF A HERDER EDIT SUBROUTINE. 

264 EDITS DATA FILE HERDER. 

266 . 

268 Rf<R/B) EDIT 

269 . mwsmmmmmmmmmmmm 

270 . * THIS SUBROUTINE CALLED * 

271 . * AT LINE 836 IN - SAVE * 

272 . DATA ROUTINE : ALSO # 

273 . * CALLED AT LINE 939 * 

274 . 

280 START OF HERDER EDIT SUBROUTINE 
284 EDITS SVSTEM OUTPUT CURVE HEADER. 

286 B$<A,B) EDIT 

283 . *>:!^****itf»»r*'*****it!!f:»»*if**** 

289 . * LINE 280 CALLED FROM * 

290 . * LINE 973 * 

291 . mmmmmmmmmmmmm* 

292 . 

294 . 



312 IN1TIRL12E CURVE #7 TO 1.0 

313 INITIALIZE CURVE #8 TO 1.0 

314 . 

315 . 

320 CALCULATE VALUES FOR CURVE #7 
323 UPPER DETECTOR. 

325 . 

327 . 

323 . 

330 . DISPLAV DATA AS CALCULATED. 

331 TIME DELAV FOR READABILITV. 

332 . 

333 . 

340 CALCULATE VALUES FOR CURVE #8 
343 LOWER DETECTOR. 

345 . 

347 . 

343 . 

350 DISPLAV DATA 

351 TIME DELAV 

352 . 

353 . 

355 TIME DELAV FOR VIEWING SCREEN. 

360 PUT MESSAGE OH LINE 102 OF MAIN MENU. 
370 . 


400 . 

401 . 

402 . 

403 . 

410 DEFAULT VALUE. 

411 DEFAULT VALUE. 

412 DEFAULT VALUE. 

413 . 

416 CHANGE VALUE IF DESIRED. 

417 CHANGE VALUE IF DESIRED. 

413 REJECTS INPUT OF LESS THAN 20. 

413 CHANGE VALUE. 

420 OUTPUT CONSTANT * TOP DETECTOR * AMPS OUT PER WATT IN. 

421 CHANGE IF DESIRED. 

422 OUTPUT CONSTANT * DOT DETECTOR if: AMPS OUT PER WATT IN. 

. 423 CHANGE IF DESIRED. 

424 W3 IS STARTING WAVELENGTH OF TABLE. 

425 K3 ALTERS CURVES OF BOTH DETECTORS BV SAME MULTIPLIER. 
428 . 

423 PUT INFO IN HERDER LINES 3/ 4, 6. 

^ 430 . 

431 ENTER V FOR HARD COPV. 

432 SET HARD COPV FLAG. 

434 TURN ON PRINTER * USE PRINT#4 COMMAND TO PRINT DATA. 

435 • 

436 Cl FOR EQUATION ' LINE 450. 

437 C2 FOR EQUATION LINE 450. 

438 DISPLAV PART OF DISPLAV HERDER. 

433 ROUTINE TO PRINT HERDER INFO ON HARD COPV. 

440 T « STARTING TEMPERATURE. 

441 K = STARTING TEMP. IN DEG. K . 

442 . 

443 . 

444 

445 WT s TOP DET ACCUMULATOR * WU ■= EOT DET ACCUMULATOR. 

446 . 

447 . 

443 CALCULATE DET. OUTPUTS FOR EACH TEMPERATURE. 

443 CALC WAVELENGTH = V. 

450 PLANKS LAW CALC FOR BLACK BODV OUTPUT. 

452 MULT BV REL RESP. AND BV RESP AT PEAK. 

453 MULT BV REL RESP. AND BV RESP AT PEAK. 

454 ADD TO SUBTOT. 

455 ADD TO SUBTOT. 

456 . 
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438 ENTER TfiBLE VRLUE * OUTPUT CURRENT OF TOP DET. 
460 ENTER TABLE VALUE * OUTPUT CURRENT OF BOT DET. 

467 HARD COPV PRINTOUT. 

468 DISPLAV DATA. 

474 CHECK FOR END OF TABLE. 

475 CHECK FOR END OF TABLE. 

478 INCREMENT TEMPERATURE. 

479 INCREMENT TABLE LINE POINTER. 

480 CALCULATE NEXT VALUE IN TABLE. 

482 CHANGE DISPLAV MESSAGE FOR LINE 102. 

484 TURN PRINTER OFF IF IN USE. 

485 RETURN TO MAIN MENU WHEN READV. 

490 REM ?7?7 

492 REM 7777 
494 . 


■ 

> 


o 


500 . 

502 START OF SUBROUTINE FOR HARD COPV OF SVS. OUTPUT CURVES. 
504 . 

506 . 

507 CHANCE TO CORRECT ERROR. 

508 . 

509 . 

510 . 

511 . 

512 . 

513 . 

514 . 

515 . 

516 . 

517 . 

518 . 

519 . 

520 . PRINTER TAB COMMANDS * SEE BOOK. 

522 . ASSUMES USE OF COMMODORE 801 PRINTER. 

523 . 

324 . 

523 . 

526 . 

527 

523 END OF HEADER SUBROUTINE. 

550 . 

551 . 

552 . 

353 

554 END OF DATA PRINT SUBROUTINE. 

560 DISPLAV HEADER. 

561 . 

562 . 

564 . 

570 DISPLAV DATA. 

572 . 

574 . 


600 CLEAR SCREEN. 

602 mm*m data table menu ***iK**# 

604 . 

606 INPUT DATA'ROUTINE 0 LINE 700 VIA LINE 651. 

607 INPUT DATA ROUTINE 0 LINE 700 VIA LINE 653. 

608 INPUT DATA ROUTINE 0 LINE 700 VIA LINE 655. 

609 INPUT DATA ROUTINE 0 LINE 700 VIA LINE 657. 

610 INPUT DATA ROUTINE 0 LINE 700 VIA LINE 659. 

611 INPUT DATA ROUTINE 0 LINE 700 VIA LINE 661. 

612 INPUT FROM DISK 0 LINE 1200 

613 TO EDIT MENU 0 LINE 1000. 

615 . 

616 tttttTft END OF MENU tttttttttTttt 
624 . 

626 . 

623 

630 MAKE MENU SELECTION. 
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633 

634 MENU SELECTION RDDRESSES. 

651 ASSIGN HERDER TITLE. 

652 ASSIGN FILE NUMBER. 

653 ASSIGN HEADER TITLE. 

654 ASSIGN FILE HUMBER. 

655 ASSIGN HERDER TITLE. 

656 ASSIGN FILE NUMBER. 

657 ASSIGN HEADER TITLE. 

653 ASSIGN FILE NUMBER. 

659 ASSIGN HEADER TITLE. 

660 ASSIGN FILE NUMBER. 

661 ASSIGN HEADER TITLE. 

662 ASSIGN FILE NUMBER. 

700 START DATA EHTRV ROUTINE 

710 INIT ITEM SUBSCRIPT NUMBER. 

712 INPUT VALUE. 

714 KEEPS OLD VALUE IF NO ENTRV MADE. 
716 ESCAPE OR MAKE CORRECTION 
718 . 

720 INCREMENT SUBSCRIPT. 

721 TO MAIN MENU IF TABLE COMPLETE. 

722 TO - ENTER NEXT LINE. 

790 INPUT LINE TO BE CORRECTED. 

791 IF NO ENTRV ABOVE - GOTO MENU. 

792 RESET SUBSCRIPT NUMBER. 

794 TO - ENTER NEXT LINE. 

793 NOT NEEDED. 


300 . 

801 . 

802 SUBROUTINE AT LINE 2000 PERMITS LISTING. 

803 OR PRINTING OF THE DISK DIRECTORV. 

804 *** 

803 . 

812 . 

313 . 

314 . 

815 . 

816 . 

317 . 

819 . 

824 . 

825 . 

826 . 

327 . 

823 . 

829 tttTtttt END OF MENU tttttttt 

830 INPUT NUMBER OF CURVE TO BE SAVED. 

332 ASK FOR COMMAND IF PREVIOUS INPUT 

333 WAS OUT OF LIMITS. 

834 COMMAND ADDRESSES. 

836 SUBROUTINE AT 260 PERMITS CHANGE OF FILE NAME. 

837 . 

842 EXIT TO DATA TABLE MENU IF DESIRED. 

843 . 

844 ENTER VALUES IN HEADER STRINGS 1 AND 3. 

845 SCRATCH FILE FROM DISK AND CHECK RESULT. 

846 CLOSE COMMAND CHANNEL. 

848 SET UP TO SAVE FILE. 

850 SAVE FILE HEADER. 

851 . 

852 . 

860 SAVE 96 ITEM DATA FILE. 

862 . 

864 . 

866 END OF FILE MARKER. 

863 CLOSE CHANNEL 2. 

870 . 


SAVE DATA SET MENU 


************************* 

* UPDATE FILE RECORD - * 

* CALLS LINE 1510 WHICH * 

*• CALLS SUBROUTINE AT * 
■« LINE 1600 . LINE M 
» MENU * 

mmmmmmmMmmmMm* 





900 . 

904 , • 

905 MENU m SRVE SVS OUT CURVES 'mm* 
908 . 

910 . 

912 . 

914 . 

915 . 

916 . 

929 tttTmt END OF MENU ftttttttttt 

932 • 

933 PROTECTS flORINST OUT OF LIMIT COMMAND. 

934 . 

935 CALLS DIRECTORV LIST SUBROUTINE. 

936 COMMAND ADDRESSES. 

938 i 

939 CALLS NAME CHANGE SUBROUTINE. 

940 . 

944 RETURN TO ABOVE MENU OPTION. 

945 . 

946 PUT VALUES IN HEADER LINES 1 AND 3. 

947 SCRATCH DISK FILE. 

948 VERIFV AND CLOSE COMMAND CHANNEL. 

949 . 

950 OPEN FILE «2 TO SAVE FILE. 

954 SAVE HEADER. 

955 INSERT EXTRA CHARACTER IN LINE TO. 

956 AVOID NULL STRING. 

960 SAVE 96 ITEM DATA SET. 

962 . 

966 . 

970 END OF FILE MARKER. 

972 CLOSE FILE - RETURN TO MENU. 

973 CALL HEADER EDIT ROUTINE. 

974 . 

975 . 

978 . 

980 ESCAPE TO MENU. 

982 CONDITIONAL RETURN TO MENU 

983 ENTER HERDER INFORMATION 

984 SCRATCH EXISTING FILE 

985 READ DISK DRIVE STATUS 

986 CLOSE COMMAND CHAN. 

988 OPEN FILE TO WRITE 

990 WRITE HERDER TO DISK 

991 . 

992 . 

993 WRITE DATA TO DISK 

994 . 

995 . 

33S 

997 END OF RECORD MARKER 

998 CLOSE FILE 

999 . 


1000 . 

1010 . 

1012 WUllllitDtlWltiiKltliKitt EDIT MENU 
1014 . 

1016 . 

1017 . » LIST OF ALL OF THE DATA * 

1018 . n USED IN THIS PROGRAM. * 

1019 . 

1020 . 

1021 . 

1022 . 

1023 EXAMINE FILE READS FILE FROM DISK AND 

1024 DISPLAYS IT. IT IS NOT PUT INTO MEMORY. 

1025 . 103 

1028 tttmnt END OF MENU tTtttttTt 

1029 SELECT CURVE TO BE EDITED. 



1030 

1032 

1033 

1034 

1035 

1036 

1037 

1038 

1039 

1042 

1043 

1044 

1045 

1046 

1047 

1048 

1049 

1050 

1051 

1052 

1053 

1054 

1055 
1053 
1060 
1061 
1062 

1063 

1064 

1065 

1066 


. . SUBSCRIPT D IS HUriBER 
. . OF DRTR ITEMS IN FILE. 
. . D IS USED TO SET LOOP 
. . D IS USED TO SET LOOP 
. . LENGTH 


LINE 2 OF FILE HERDER. 
DISPLRV FILE HERDER. 


INPUT S VBSCRIPT FOR DRTR LIST- 

DISPLRV 10 DRTR ITEMS STRRTING RT RBOVE SUBSCRIPT. 
CHECK FOR END OF LIST. 

IF NO ENTRV/ DISPLRV NEXT 10. .ITEMS 

CHRNGE HRS BEEN MRDE - REDISPLRV ITEMS. 

IF RT END OF TRBLE. RETURN TO EDIT MENU. 


1067 . 

1063 CHRNGE COMMAND NUMBER TO FILE ID NUMBER. 
1069 EDIT SVSTEM OUTPUT CURVES 

1071 . 

1072 . 

1073 . 

1074 . 

1075 . 

1076 . 

1073 . 

1079 . 

1080 . 

1032 . 

1084 . 

1035 . 

1086 . 

1007 . 

1083 . 

1089 . 

1096 . 

1097 . 

1099 . 

1105 . 

1110 . 

1111 . 

1112 . 

1113 . 

1121 . 

1122 tm END OF EDIT ROUTINES ntt 
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1200 

1203 

1204 
1206 
1208 

1209 

1210 
1220 
1221 

1223 

1224 

1230 

1231 


INPUT CURVE FROM DISK 

■ 

SUBROUTINE fiT 2000 LISTS DISK D I RECTOR V. 


ESCRPE TO fIflIH MENU IF DESIRED. 



1232 OPEN COnnflND CHANNEL SO THAT ERRORS CRN BE READ. 

1234 OPEN FILE FOR READ. 

1235 READ HEADER. 

1236 DISPLRV HERDER AS READ. 

1237 . 

1238 SET D1 EQUAL TO FILE NUMBER. 

1239 SET D3 EQUAL TO NUMBER OF DATA ITEMS. 

1240 SUBROUTINE AT 1310 CHECKS FOR DISK ERROR & CLOSES FILES. 

1250 OPEN FILE FOR READ. 

1251 SSELECT CORRECT READ ROUTINE. 

1254 SSELECT CORRECT READ ROUTINE. 

1258 READ HEADER. 

1260 READ DATA. 

1262 . 

1263 DISPLRV DATA AS READ. 

1267 . 

1268 . 

1270 GOTO MAIN MENU 

1271 READ HEADER. 

1272 READ DATA. 

1274 . 

1279 . 

1281 . 

1284 GOTO MAIN MENU. 

1310 SUBROUTINE - READ DISK STATUS. 

1312 CLOSE CHANNELS. 

1314 IF FILE EXISTSi CONTINUE. 

1316 IF ERROR, PRINT TVPE OF ERROR AND START OVER. 

1320 . 

1321 ESCAPE IF DESIRED. 

1329 SET D1 EQUAL TO FILE NUMBER. 

1330 VERIFY D1 


. LIST OF FILES BV NUMBERS . 
. AND FUNCTION. 


Tttm END OF LIST ttTTtT 


1350 INITIATE INPUT OF FILE FROM DISK. 


S> 
• \ 


1332 . 

1333 . 

1334 . 

1335 . 

1336 . 

1337 . 

1338 . 

1339 . 

1340 . 
1344 . 
1346 . 
1348 . 


1400 . 

1402 IK****** PRINT DATA MENU ifi*#**#* 
1404 . 

I486 . 

1407 . 

1408 . 

1409 . 

1410 . 

1411 . 

1412 . 

1413 . 

1415 . 

1416 . 


1417 ttnnmt end of menu ttTnttTrt 

1420 SELECT FILE TO BE PRINTED. 

1422 CHECK FOR OUT OF RANGE. 

1424 CHECK FOR OUT OF RANGE. 

1426 . 

1428 GOTO PRINT ROUTINE FOR ALL FILES EXCEPT SVS. 
1430 GOTO PRINT ROUTINE FOR SYSTEM OUTPUT ONLY. 
1434 ENABLE PRINTER. 

1436 PRINT FILE NAME. 

1440 PRINT INFO. 

1441 PRINT FIRST 4 HEADER LINES. 

1442 SUBROUTINE AT 1700 SELECTS TITLE. 

1443 START DATA PRINTOUT. 


OUT. 
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1444 . 

1446 CnLCULflTE WAVELENGTH FOR DATA ENTRV. 

1443 TAB AND PRINT. 

1450 . 

1452 DISENABLE PRINTER 

1456 RETURN TO PRINT DATA MENU. 

1460 ENABLE PRINTER. 

1461 . 

1462 . 

1463 SET, TEMPERATURE DATA FROM FILE HEADER. 

1464 VERIFV LIMITS AND INCREMENT OF TEMPERATURE. 

1465 . 

1466 . 

1468 . 

1470 . 

1472 START LOOP TO PRINT DATA TABLE. 

1474 CALC. TEMPERATURE. 

1476 . 

1473 TAB AND PRINT 
1480 TAB AND PRINT 

1431 CHANGE 0 TO 1 PREVENT DIVISION BV ZERO ERROR. 

1432 TAB AND PRINT 

1483 CHECK FOR END OF PRINTOUT. 

1434 SECOND CHECK FOR END OF PRINTOUT. 

1435 . 

14?^8 

1490 DISENABLE PRINTER. 

1492 RETURN TO PRINT DRTFl MENU 


1700 SUBROUTINE SELECTS PROPER FILE NAMEISSS 

1711 FOR STRING VARIABLE. 

1712 . 

1713 . 

1714 . 

1715 . 

1716 . 

1718 . 

1719 . 

1730 tttt END OF SUBROUTINE tttt 


2000 SUBROUTINE READS DISK DIRECTORV m 
2002 . 

2004 . 

2006 THIS ROUTINE WAS FOUND IN A BOOK. 

2007 NO EXPLANATIONS WERE GIVEN. 

2008 . 

2009 . 

2010 . 

2015 . 

2020 . 

2030 . 

2040 . 

2050 . 

2060 . 

2070 . 

2080 . 

2031 . 

2090 . 

2100 . 

2110 . 

2111 . 

2112 . 

2113 . 

2114 . 

2120 . 

2130 . 

2140 . 

2150 . 

2151 . 

2160 . 

2164 . 

2166 . 

2168 . 

2170 . 
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2200 . 

2201 . 

2205 . 

2206 . 

2210 . 

2300 . 

2310 . 

2320 ttnt END OP DISK DIRECTORV SUBROUTINE ttm 


2400 m«m nuLTiPLV svs out curve bv cohstrnt w.mm 

2402 . 

2404 . 

2406 KT = TOP DET. CONSTANT i KB = EOT CONSTANT. 

2403 . 

2418 . 

2412 CHANGE HEADER LINE 7 TO INDICATE CONSTANT IN USE. 
2420 nULTIPLV ALL DATA BV NEW CONSTANTS. 

2422 . 

2423 . 

2424 . 

2425 . 

2426 TIME DELAV LOOP 
2460 GOTO MAIN MENU 

3000 OPEN COMMAND CHAN. 

3010 READ DISK STATUS. 

3020 CLOSE COMMAND CHAN. 

3030 DISPLAY DISK STATUS. 

3040 . 


3100 ***** EXAMINE FILE LOOKS AT ONE CHARACTER AT A TINE. 
3110 . 

3120 SEQUENTIAL, PROGRAM OR USER. 

3130 DISCARD ALL BUT FIRST LETTER. 

3140 TEST FOR VALID TYPE. 

3145 OPEN ERROR CHANNEL. 

3150 OPEN FILE 

3160 CHECK FOR EXISTANCE OF FILE. 

3170 READ A CHARACTER. 

3188 . 

3135 . 

3187 . 

3190 . 

3191 . 

3200 READ STATUS DATA. 

3210 IF ERROR, PRINT STATUS. 

3220 . 

3290 . , 

3291 . 


3300 t*r*T END OF PROGRAM REMARKS r*'r*1*1* T*1*t 

330 1 ttTttttttttTttttTfrtttTtTtTtttTttttttttttt 

3302 . 

3303 ******************************************************* 

3304 NEXT IS AUTO LIST WHICH RESPONDS TO A RUN COMMAND. 
4000 ******************************************************* 
8000 print: PRINT: print: PRINT 

8002 PRINT"********* REMARKS FOR PROGRAM ********" 

8003 PRINT''************** WHT-CA96F *************” 

8004 print: PRINT :PRINTTAB<6) "THIS LIST IS DATED 11/16/85" 

8006 PRINT :PRINTTAB<6) "THE PROGRAM IS DATED - REV. 12/09/85 
8008 PRINT : PRINT : PRINT : PRINT 

8020 INPUT"ENTER V FOR HARD COPV'jQf 
8022 IFQf="V"THEN0PEN4,4,0!CMD4: GOTO 8050 
8024 LIST 

8050 PRINT ! PR I NT : PR I NT : PR I NT 

8052 PRINT"********* REMARKS FOR PROGRAM ********" 

8053 PRINT"******:******** WAT-CA96F ***:**********" 

8054 print: PRINT :PRINTTABC6) "THIS LIST IS DATED 12/16/85" 

8056 PRINT :PRINTTAB(6) "THE PROGRAM IS DATED . r. . 12/09/85 
8058 PRINT! PRINT : PRINT : PR I NT 

8060 LIST 
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APPENDIX D3 


COMPUTER PROGRAM A/D READ-MIO 



1 REM PROGRAM fl/D READ-M10 

2 REM DATE 5/17/85 

3 POKE55.00!POKE56,145:REM LOWER TOP OF BASIC TO PROTECT MACODE 

4 REM CHANGED ROUNDOFF - AVOID STORAGE IN FLOAT PT OF FINAL ANSWER 

5 REM LINES1700-1740 ARE OK NOW 

7 REM ALT ROUNDOFF - LINES 2000-2018 

8 REM USING LINES 1700-1740 FOR ROUNDOFF 

9 DIM CF<20)/GF<20).DD<20)/B*<20)>E*a2),N#<12) 

10 DIM A4<255) 

11 P0KE57343,2 

12 DIM MA<256) 

14 DIM CC<3.90) 

15 DIM D<18)»Ga8>^Q:«18) 

16 SK»1 

17 G<0>»1 

18 G<1)=1 

20 0<2)-l.9531E-6 

21 G<3)»1.0255E-7 

22 G<4)»3. 1091E-9 

24 G<5>=9.3614E-2 

25 G(6)=4.8603E-7 

26 G<7>a2.1028E-8 

27 O<8)“9.821E-10 

28 G<9)-9.2525E-2 

29 REM ENTRV DATE 3-12-85 

30 G<12)»G<2)*G<5) 

31 Ga3)=Q<3)»0<5) 

32 0a4)=G<4)»G<5) 

33 0<16)«0<6)*0<9) 

34 0<17)-G<7)»0<9) 

35 G(18)-0<8}«0(9> 

64 B#<0)«"CHAN ZERO" :CF<0)»0 : GF<0)=8 : DD<0)=0 

65 B»<l)-"TOP LOW10X " :CFa)=3:GFa)«0:DD<l)«14 

66 B#<2)»"T0P LO " : CF<2)»1 ! GF(2)=>0 •• DD<2)® 4 

67 B*<3)-"T0P MID«10X":CF(:3)=3:GF<rj)ai:DD<3)=13 

68 B*<4)-"T0P MID " :CF(4)-1 :DF<4)»1 :DD<4)« 3 

69 B*<5)»”T0P HI*10X " :CF<5)=3 : GF<S)«2 • DD<5)=12 

70 Bf<6)“"T0P HI ":CF<6)»1 :GF<6)»2:DD<6)» 2 

71 B»<7)«"B0T LO»10X “ :CF<7)"4 : GF<7)»0 : DB<7)=»18 

72 BSC8)»"B0T LO " :CF<8)=2:GF<8)a0:DD<8)» 8 

73 B#<9)»"B0T MID»10X":CF<9)-4:GF<9)»1 :DD<9)=17 

74 B«(10)="BOT MID ” :CF<10)=2:GF<10)=1 :DDd0)= 7 

75 B$<in="BOT HI#10X " : CF< 1 1 )*4 : GF< 1 1 )»2 = DDCl 1 )=16 

76 B*U2>-"B0T HI " : CFa2)»2 : GF< 12>=2 ; DD< 12)=« 6 

80 IF(PEEK<38143)K>7THENIJU»*"LD" : GOTO350 

81 P0KE38142,! 

86 QOTO200 


100 PRINTCHR#< 147) = PRINT : PRINT : PRINT"***#*#*******#” ; 

101 PRINT" MAIN MENU ***************" 

102 print: PRINT: PRINT 

104 PRINTTAB<3)"1 - READ PROBE OUTPUT " 

106 PRINTTAB(3)"2 - LOAD TEMPERATURE RATIO CURVE" 

108 PRINTTAB<3)"3 - LORD MACH. PROGRAM" 

110 PRINTTRB<3)"4 - TRANSFER SEQ FILE TO ANOTHER" 

1 1 1 PRIHTTAB< 10) "DISK" : PRINT 

112 PRINTTAB<3)"5 - GET ZERO READINGS " 

114 PRINTTRB<;3)"6 - READ A/D - ANV SINGLE CHANNEL" 

116 PRINTTRB(3)"7 - PRINT ZERO READINGS" 

118 PRINTTRB<3)"0 - INIT ZERO RDG. TABLE" 

119 PRINTTRB(3)"9 - CONTINUOUS READ - NO HERDER" 

120 PRINT:pRINT:PRINTTRB<12)J :U»0: INPUT"COMMRND";U 
122 IFU»OTHEN10O 

124 IFU>9THEN100 

126-0NO-OOTO-180/2007-30O,-24e0.500;-600;160,800-/170 

160 OOSUB1O00 
162 GOTO 100 

170 WW»88:FL=»1 

171 OOTO709 

180 WW-0:FL»O 

181 OOTO700 
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200 PRINTCHR»< 147) SPRINT 

201 PR I NT "**«*» LOAD TEMPERATURE RATIO CURVE ♦»*!«*" 

202 Kf="SVSOUT-02" 

203 PRINTTAB<5)”FILE TO BE LORDED IS "K» 

204 PRINT :PRINTTRB<3) "*»»»#*****♦ MENU Jli#***#*#*#!**" = PRINT 

205 PRINTTRB<5)" 9 - MAIN MENU" 

206 PRINTTRB<5)" RETURN - CONTINUE " 

207 PRINTTRB<5)" C - CHANGE FILE NAME" 

208 PRINT :L4:»"" 

209 PR I NTTRB < 5 ) ; : I NPUT " COMMAND " ; L$ 

210 IFLSa"9"THEN100 

211 IFL»="C"THEN213 

212 G0T0220 

213 PRINT : INPUT"CHRNOE NAME TO "iKS 
220 INPUT'TRESS RETURN WHEN READVSR* 

222 PRINT !PRINT"L0AD1NG FILE "K$ " - PLEASE WRIT" 

226 OPEN2/8/2,"0:"+K#+’SS,R" 

228 F0RI»1T08 

230 IHPUT»2.CR» 

232 AWf(I)>MID«(CASi2> 

234 NEXT 

233 CL0SE2 

236 PRINT:IFAW*<2)OK»THENCLOSE2:PRINT"FILE "KJ" NOT FOUND" :GOTO204 
233 PRINT:PRINT"***#*r<intiiltiil(itiHl FILE HERDER 

239 PRINT 

240 FORUalTOS 

241 PRINT"RW»<"U")» "RW»CU) 

242 NEXT 

243 OPEN2/8,2."0:"+K#+"iS/R" 

246 F0RI=1T03 

247 INPUT#2.CR$ 

249 NEXT 

231 PR I NT "LORDING DATA" 

232 F0RI=1T084 

254 INPUT#2,CC<0/I) 

256 INPUT#2,CC<1,I) 

258 NEXT 

260 CL0SE2:pRINT"CALCULATING TABLE" 

261 TS=VRL<AW»<4)):REM LO TEMP 

262 TD»VHL<RW*<5)):REM TEMP INCREMENT 

263 TE=»VRL<RW$<6)):REM HIGH TEMP 

264 F0RI=»1T084 

265 IFCC<1,I)=0THENCC(1.I)»1 

266 CC<2. I)»TS+<<I-1)*TD) 

268 CC<3,I)»CC<0/l)/CCa,I) 

269 00T0272 

270 PRINTCC(2>I)” "CC(3,D 

272 I FCC < 2 . 1 ) «>>TETHEN275 
274 NEXT 

273 print: PRINT 

276 PRINT"RRTIO AT "CC<2,1)" DEG.= "CC<3,D 
273 PRINT"RATIO AT "CC<2a>" DEG.» "CC<3,D 
280 FORIalTO2500:NEXT 
299 GOTO 100 


300 PRINTCHR$< 147) : PRINT: print: PRINT 

302 PRINTTRB<3)"iHX<»*)|i»i**ill>t( LORD MRCODE 

303 print: print: PRINT 

328 PRINTTRB(5)"PRESS RETURN TO LORD MRCODE" 

330 V TO ESCAPE TO MENU";K$ 

- 332 IFK»»"V"THEN100 -- ■ --- " * 

333 PR I NT: PR I NT 

330 PRINT"CHECKING FOR PRESENCE OF MRCODE" 

331 P0KE38143,7 

332 0PEN2i 8> 2, "0: MRCODE/ S/R" 

353 F0RI-1T08 

354 INPUT#2/D$ 

355 R*<I)=MID*<D*/2) 

356 PRINTR$<I) 

337 NEXT 

333 IFR|:<2)O"MRC0DE"THENCL0SE2:pRINT"FILE NOT FOUND" : P0KE38 143/ 0 :GOTO390 

359 PRINT"RERD1NG DISK AND LORDING MEMORV" 

360 J=37888 



362 F0R1=IT0255 

364 INPUT#2.fl4 

365 P0KEJ.H4:J=J+i:PRINTJ" ”Fl4 

366 NEXT 
368 CL0SE2 

336 1 FuijJ = ” LD " THEHUU$= " ” : G0T08 1 
387 GOTO 100 

390 PRIHT-L$=»"**:iNPUT"RETRV IF V 
392 IFLl:='‘V'THEN350 
394 GOTO 100 


500 PRINTCHR$< 147) ‘PRINT SPRINT 

502 FRlHV'mMm ZERO READINGS ON ALL RANGES 'Mmm 
504 print: PRINT :PRINT”PLEASE WAIT" 

508 F0RI=1T012 
530 P0KE57343.GFd) 

532 FORU*1TO450:NEXT 

533 VB=»0 

534 F0RT«1T012 

535 CHcCFd) 

536 COSUB 1 145 

537 IFT<3THENR4=0 : 00T0542 

538 VB=»VB+R4 

540 PRINT"R4» ••R4“ VOLTS*' 

542 NEXT 

546 D(DDa))=<VB/10) 

547 PRINT'*DC‘DDa)")» **D<DDa))‘* “Bf 

550 NEXT 

551 FORU=lTO70O:NEXT:PRINTCHR$a47) :PRINT:PRINT:PRINT 

552 F0RI=1T012STEP2 

554 PRINTB$d)"« "DCDDa))*’ VOLTS 

555 PRINTB«(I + 1)”= “DCDDd + l ) ) '* VOLTS 

556 PRINT 
558 NEXT 

. 567 PRINT 

569 PRINT"READING TIME** : GOSUB1400 

570 P0KE33142.9:REn ZERO RDG. FLAG 

571 ND$=VD$ : NSJE=HSf 

572 K$a«»;lHPUT**DO OWER IF V ";K$ 

574 IFK$«'*V*'THENP0KE38142/ 0 : GOTO500 
576 A$a*'** : INPUT"PRINT LIST IF V ";fif 

573 IF A$-"V"THENGOSUB1800 
580 GOTO 100 


600 PR I NTCHR$ a 47 ): PR I NT 

602 PRINT‘*#4i*#5t( READ A/D # INDIVIDUAL RANGES mmr* 
604 MENU 

606 PRINT 

610 PRINT" 0 - CHANNEL ZERO" 

612 F0RI=1T012STEP2 

614 PRINTI" - "Bi<I)‘* : "I + l" - “BS^I + D 
616 NEXT _ 

613 PRINT ' ■' 

620 PRINT" 90 - MAIN MENU 

621 PRINT 

624 IFK-90THEN100 

625 IFK>12THEN600 

641 IFK>12THEN100 

642 GT=GF<K):CHkCF<K) 

660 PRINT 

662 GOSUB1120 

666 PRINT'*####3|t.4i "B^(K)" RANGE mmm” 

668 XRaR4 

669 *v'A=R4itiG<DD(K> > : VAl:=STR$<'^R) : VE$="AMP** 

670 IFK=0THEN'v*A3:=" " : VB$=" " 

671 PRINT"OUTPUT IS "XA" *v*OLTS" : PRINT" » "VA«" '*WBt 

672 PRINT 

680 K$=""' INPUT "MR IN MENU IF W'*;KS 
682 IFKJ="'t'"THEN100 
633 GOTO600 
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700 REM 

701 PRIHTCHRia47):PRINT:PRIHT:PRINT"##s(!**iti*i(i READ PROBE OUTPUT 

702 print: IFPEEKO 8142)=lTHENPRINT"ZER0 READINGS HAVE NOT BEEN TAKEN" :GOTO790 

703 INPUT" PRINTOUT WANTED IF V "JK$: PRINT 

704 FL=0: IFKf="V"THENFL=l 

705 IFFL=0THEN707 

706 G0SUB1565 

707 PR INTTABC3) "PRESS RETURN TO 

708 INPUT"START READING" ;K$ 

703 PRINT; PRINT" READING - PLEASE WAIT" 

710 GOSUB302 

712 PRINT"TOP DET OUTPUT » "XT" AMPS" 

722 PRINT"BOT DET OUTPUT * "XB" AMPS" 

724 I FXB=0THENXBa 1 : G0T0728 
726 GOTO730 

728 PR I NT: PR I NT "DATA m ERROR" 

730 RT=XTXXB 

732 BTb< INT< <XT/XB)*10000) )/'10000 

744 PRINT :PRINT"RATIO OF OUTPUTS IS "BT 

745 GOSUB1200 

747 GOSUB1240 

748 Wl$»STRf<TT) : IFTT$="<"THENWl»=MIDf CAWf <5)/2) 

749 I FTT* a " > " THENW 1 $=M I Df ( A Wf < 6 ) .■ 2 ) 

754 W2$=STR#<BB) : IFBBt="<"THENW2$«=MID*<AWJ:<5) . 2) 

755 IFBBi;*">"THENW2$»MID$< AW#<6) / 2) 

762 IFVRLCRWf<3)KlTHENWU:a"N0 TRBL" :W2$="N0 TRBL" :TT$=»"" :BB»="" 

764 W6*=3TR$<TM) 

766 IFVRL<RWf<3)KlTHENTMl:="":W6$="N0 TRBL":TM*="" 

767 IFTM$="<"THENW6$=MID*(flW*<4> , 2) 

768 IFTM$=">"THEHW6«=MID* (AW* C 6) / 2> 

780 PRINT"T-TOP«» "TT*+W1*" T-BOT= "BB$+W2* 

781 PRINT" INDICATED TEMP* "TM*+W6:f 

783 IFFL=0THEN786 

784 GOSUB1500 

785 IFWW*8aTHENG0T07S3 

786 PRINT;K»*"":INPUT"MENU IF V % OTHERWISE RESET" ;K* 

787 IFK*="V"THEN100 

788 GOTO703 

790 PRINT :PRINTTRB< 10); 

794 K»*"" : INPUT" MENU IF V ";K* 

795 IFKJ:="V"THEN100 

796 PRINT 

797 GOTO703 


800 PRINTCHR#< 147) "«******» INIT ZERO RDG. TABLE #***#*♦" 

810 F0RI=2T04 

012 Da)=0 

814 Da+10)*0 

816 NEXT 

320 F0R1*6T08 

822 D<I)*0 

824 D(I+10)*0 

826 NEXT 

829 P0KE33142.1 

330 GOTO 100 
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900 REM FIND BEST RANGE FOR TOP DET. 
902 T»0:TH»3 

904 OOSUB1100 

905 IFR4>9THEN910 
908 GOTO950 

910 TH«1 
912 OOSUBllOO 
914 IFR4>9THEN910 
916 GOTO950 
918 0T*1:TH*3 
920 COSUB 1100 
922 IFR4>9THEN926 
924 OOTO950 
926 TH»1 
928 GOSUB1100 
930 IFR4>9THEN934 
932 GOTO950 



934 0T=2:TH»3 
936 COSUB 11 00 
938 IFR4>9THEN942 
940 OOTO9S0 
942 TH=1 
944 OOSUB1100 

946 IF<R2+R3)>3630THENPRIMT"PUT OF RANGE" 

947 GOTO950 

948 FORV-1TO1000!NEXT:OOTO700 

950 BH-4 

951 OOSUB1110 

952 OOSUB1400 
954 IFR4>9THEN958 
956 OOTO970 

958 BH*2 

970 VT-0:VB"0 

971 REM READ 10 ALTERNATE READINGS OF EACH DET. 

972 FORL-1TO10 

973 CH»BH 

974 00SUB1145 

975 IFR4>9.90THENPRINT"BOTTOM OUT OF RANGE" : FORM- I TO 1000 : NEXT! GOTO700 
977 VB-VB+R4 

979 CH-TH 

980 GQSUBIMS 

982 IFR4>9.98THENPR1NT"T0P OUT OF RANGE" •' FORM* 1 TO 1000: NEXT: GOTO700 
984 VT-VT+R4 

990 NEXT 

991 REM NEXT SELECT GAIN SUBSCRIPTS 

993 PRINT 

994 IFTH«3THENQR«<4-OT)+10 

995 IFTHalTHENQR»(4-GT) 

996 IFBH«4THENLV»<8-OT)+10 

997 IFBH»2THENLV«(8-GT> 

1002 PRINT"TOP OUT= "VT/10" VOLTS RANGE"QR 

1003 PRINT"BOT OUT- "VB/10" VOLTS RANCE"LV 

1008 REM CALC. CURRENT OUT 

1010 XT-<(VT/10)-LCQRJ)*G(QR) 

1020 XB-<<VB/10)-D(LV))iKG(LV) 

1030 RETURN 


1100 REM SELECT INPUT CHANNEL 
1102 CH=TH 
1104 GOTO1120 
1110 CH-BH 

1119 REM SET AMPLIFIER RANGE 
- - 1120- P0KE57343,-QT 

1123 REM DELAV AFTER RELAV SWITCH 

1124 REM 300 COUNT LOOP IS 350 MSEC. 

1125 FORI- 1TO300: NEXT 
1140 REM READ A/D 
1145 P0KE57342,CH 
1150 A-57340 

1155 P0KEA,16 
1160 POKER, 23 
1165 POKER, 16 

1170 IF<PEEK<57342)AND2)=2THEN1 170 
1172 R2-0:R3-0 
1175 POKER, 24 
1130 R3=PEEK<R) 

1185 POKER, 25 

1186 R6-PEEK<R> 

1187 R7-1 

1188 1F<R6RND32)=32THENR7=-1 
1190 R2«<R6RND15)#256 

1192 R4*<R2+R3)*<10/4096)#R7 
1195 RETURN 
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1200 REM INTERPOLATION ROUTINE 

1201 TM»»'"' 

1202 IFRT<CC<3. 1)THENTM*="<" :W6i:=MID*<AUi:<5), 2) ’RETURN 

1 204 I FRT>CC <3 , 84 ) THENTM**" > " : M6J:=M I D$ < RW$ < 6 ) < 2 : RETURN 

1206 IFRT<CC<3. 15)THENI=»l :G0T01218 

1208 IFRT<CC<3,30)THENIal4:OOTO1218 

1210 IFRT<CC<3.43)THENI=29:G0T01218 

1212 IFRT<CC<3/60) THEN I »44 : GOTO 1 2 1 8 

1214 IFRT<CC<3,7S)THENI=59:00T01218 

1216 IFRT<CC<3.84)THENI=74:G0T01218 

1218 IFRT>CC<3/ I)THENIFRT<CC<3/ I+1)THEN1226 

1220 1*1+1 

1222 IFI>84THENTM#*">":W6$=MID$<HWI<6).2) :RETURN 
1224 G0T01218 

1226 RU"(RT-CC<3. D)/<CC<3, 1+U-CC<3, 1)) 

1228 TM*CC<2, I)+<<CC<2. I+1>-CC<2/ I))*RU) 

1230 TM*<INT<TM*10 ))/’10 
1232 RETURN 


1240 TTfa" " : IFXT<CC<0> 1 >THENTT*0 : TT$="<" : G0T01266 
1242 I FXT>CC<0 . 84 ) THENTT=0 : TTf * " > " : GOTO 1 266 

1244 IFXT<CC<0. 13)THENI=1 :GOTO1250 

1245 IFXT<CC<0, 30)THENI=14 : GOTO1250 

1246 I FXT<CC < 0 , 45 ) THEN I *29 : GOTO 1 250 

1247 IFXT<CC<0< 60)THENI=44 : G0T01250 

1248 I FXTCCC < 0 < 75 ) THEN I =59 : GOTO 1 250 

1249 IFXT<CC<0< 84)THENI=74 : GOTO1250 

1250 IFXT>CC(0< OTHENIFXT<CC<0< I + OTHEN1258 
1252 1=1+1 

1234 IF I>84THEHTT=0 = TT* * " > " GOTO 1 266 
1256 GOTO1250 

1238 RU=(XT-CC<0< O)/<CC<0.. I+O-CC<0< O ) 

1239 TT=CC<2< I )+< <CC<2< I+l )-CC<2< 1 ) )+:RU) 

1260 TT=(INT<TTitil0 ))X10 

1266 BB$=" •• : IFXB<CC< 1 , 1 >THENBB=0 • EB#*"<" : RETURN 


1267 IFXB>CC<1<84)THENBB=1 :BBf=">" : RETURN 

1268 IFXB<CC<1< 15)THENI=1 •G0T01274 

1269 IFXB<CC<1<30>THENI=14:GOTO1274 

1 270 I FXB<CC < 1 < 43 ) THEN I =29 : GOTO 1274 

1271 I FXB<CC< 1 < 60 ) THEN I =44 = GOTO 1 274 

1272 IFXB<CC< 1 < 73)THENI=59 : G0T01274 

1 273 I FXB<CC< 1 < 84 ) THEN I =74 : GOTO 1 274 

1274 IFXB>CC<1< OTHENIFXB<CC<1< I+1'<THEH1278 

1275 1=1+1 

1276 IFO84THENBB=0:BB$»">" : RETURN 

1277 GOTO 1274 

1278 RU=<XB-CC< 1 < O >/<CC< 1 < I+l )-CC< 1 < O ’> 
1280 BB=CC<2< I )+< <CC<2< I+l )-CC<2< I > )#RU) 
1282 BB=<INT<BB*10 ))/10 

1284 RETURN 


1400 REM ***** READ CLOCK ROUTINE ***** 

1420 SVS37913 
1422 B=37888 

1430 VR=PEEK<B+1O+<PEEK<B+12)*10> 

1432 MO=PEEK<B+9)+<PEEK<B+10)*10) 

1434 DV=PEEK<B+7>+< PEEK < B+8 ) * 1 0 ) 

1436 DA*PEEK<B+6) 

1438 HR=PEEK<B+4)+<<PEEK<B+3>-8)*10) 

1440 MN=PEEK < B+2 ) + < PEEK < B+3 ) * 1 0 ) 

1442 SC=PEEK<B)+<PEEK<B+O*10) 

1446 HR*=STR*<HR) 

1447 MNf*STR»<MN) : IFMN<10THENMN$=''0"+MN* 

1448 SCf=STRf <SC) : IFSC<10THENSCt="0"+SC* 

1449 HS$=HR*+'‘ : "+MID$<MN*.2>+'' = "+MID$<SC* < 2) 

1450 VR|:=STR»<VR) 

1451 MO*=STR*<MO> 

1432 DV|:=STR#<DV) 

1453 VD$=MID*<M0J<2)+"-“+riID*<DVf<2> + "-"+riID«<VR|:<2) 
1460 RETURN 



1300 REn PRINTOUT ROUTINE 

1504 UU=XT 

1505 GOSUB1700 

1506 GTJ:=G*:ET#=E$ 

1507 UU=0 

1508 UU=XB 

1509 GOSUB1700 

1510 GB*=G$:EB$=Ef 
1512 FORI »lTO40Gi: NEXT 
1525 OPEN4W/0 

1530 PRINT#4/GT#;ET$; 

1531 PRINT#4>CHR*<16>CHR$<49)CHRi;<50) QR; 

1532 PRINT#4,CHRT<16)CHR*<49)CHR$<54) GB*;EB$; 

1533 PRINT#4,CHR3:<16)CHR$<50)CHR$<56) LV; 

1534 PRINT#4,CHR$U6)CHR$<5nCHR#<52) TTJ+Wl#; 
1536 PRINT#4/CHR*<16)CHR*<52)CHRf<51) BB«+W2#; 
1344 PRINT#4/CHRf<16)CHRf<53)CHR»<49> BT; 

1346 PRINT#4,CHR$<16)CHR$<53)CHR*<57) TMF+W6*; 
1 548 PR I NT#4 , CHR# < 1 6) CHR$ < 54 > CHR* < 56 > HS$ 

1330 CL0SE4 
1560 RETURN 


1365 OPEN4.4,0 

1366 PRINT#4 

1 1367 PRINT#4.. "PROBE OUTPUT flS fl FUNCTION OF TIME" 

^ 1568 GOSUB1400 

1569 PRINT#4^ ■'VERR"VR"l10HTH"M0"DflV"DV 

1371 PRINT#4, "CRLIBRflTION FILE IS "RW$C2) 

1372 IFPEEK(38142)=1THEN1574 

1573 PRINT#4, "ZERO RERDIHGS TRKEN RT "NS$" HR ON "ND$ :G0T01576 
^ 1574 PRINT#4:PRINT#4,"2ER0 RERDIHGS HRVE NOT BEEN TRKEN" •PRINT#4 

1576 PRINT#4," TOP DIODE R BOT DIODE R <T TEMP) (B TEMP)"; 

1373 PRINT#4," RRTIO IND.TENP TIME" 

1383 PRINT#4:CL0SE4 

— 1590 RETURN 


1700 REM ROUNDOFF 

1701 B»="" :BR*="" :B8#="" :BL$="" :BL=0:BTJ="" :VV=0:G$3"" -E*®"" :P5=0 

1702 B*=STR»<UU) 

1712 F0RI=1T012:N$(I)="":NEXT 

1713 F0RI=1T012 

1714 N*<I)=MID$<Bf,I,l) 

1713 IFN$<I)=""THENI=I-1 ;G0T01726 

1716 IFNta)="E"THENI=I-l :G0T01724 

1718 IFN$<I)="."THENP3=I 

1720 NEXT: 1=1-1 

1722 GOTO 1726 

1724 E»=RIGHT|:<B#,4) 

1726 IFI>6THENIFN*:<7»"4"THENN6=VflL<NS<6)) :N6=N6+1 :NJ<6)=MID$<STR$<N6).2) 
1728 E=3:IFK3THENE»I 
1730 G«=LEFT$<B*/E) 

1732 GS»0*+N#<6) 

1734 IFP3>4THENGf=">E+2" 

1740 RETURN 


1800 OPEN4,4,0:CMD4 

1810 PRINT"LIST OF ZERO READINGS" : PRINT 

1811 IFND$»“"THENPRINT"2ER0S HAVE NOT BEEN TRKEN" :PRINT#4:CL0SE4: RETURN 

1812 PRINT"DHTE "VD»" TINE "NS*: PRINT 

1820 PRINT" RANGE DESCRIPTION VOLTS AMPS" 

1822 PRINT 
1830 F0RI=1T012 

1850 PRINTCHR*a6)CHR*<48)CHR*<50)DD<I>; 

1851 PRINTCHR*<16)CHR*<49)CHR*<48)BJ:<I); 

1 852 PR I NTCHR*< 1 6 )CHR*< 30 ) CHR* < 53 ) D < DD <1 ) ) ; 

1853 PRINTCHR*< 16)CHR*<52)CHR*<49)DCDD< I ) )»G<DD< I ) ) 

1860 NEXT 115 

1864 PRINT#4:CL0SE4 
1890 RETURN 
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2000 2: VV=0K7=Q*J$ 

2001 B$=STR$CUU) : C$=LEFT$(B$, 7) 

2002 F0RI=1T08 

2003 F$a)=MID$<C$.I>l) 

2004 IFF$C I )=“E“THENI=I-1 : GOTO2007 
2003 NEXT 

2007 E$=RIGHT$CBf.4) 

2008 B$=LEFT$<Cf. I) 

2010 IFI=9THENIFF$C8»"4*‘THENK7=VfiL<F$<7>) + l :F$<7)=riID$<STR$CK7>,2> 
2012 G$=D$+J$ 

2018 RETURN 


2125 PRINTCHR$< 147) ' PRINT : PRINT : PRINT 

2126 srve grin vrlues to disk mmmr* 

2127 print: PRINT 

2128 PRINT”WILL BE SRVED UNDER NRflE: "K$: PRINT 
2130 PRINT 

2132 INPUT'XHRNGE FILE NAME TO 

2133 PRINT:LI:="‘*: INPUT^TO CORRECT - ENTER V 

2134 IFLJf=="VTHEN2132 
2136 IFKf«M«THEN2160 
2140 GOTO2051 

2160 PRINT: PR I NT "FILE NEEDS R NRME" -'PRINT 
2162 G0T02132 


2400 PR I NTCHR$ <1 47 ): PR I NT : PR I NT : PR I NT : PR I NT 
2402 PRIHT"##i*f TRANSFER SEQ. FILE DISK TO DISK 
2410 K$«""-' INPUT "ENTER NRNE OF FILE";K$ 

2412 IFK:f=""THEN100 

2414 PRINT'PRINT"FILE REQUESTED IS"Kf 

2416 PRINT :R»="":PRINT"riflKE CHANGE IF Y";R$ 

2418 IFR$=»"Y"THEN2400 

2420 PRINT : PRINT" INSERT SOURCE DISK" 

2421 INPUT"PRESS RETURN WHEN RERDV";Rf 
2424 0PEN2/ 8; 2. "0 : "+K$+" . S. R" 

2426 F0RI=1T08 

2428 INPUT#2.C:$ 

2429 fl$a)=riID$CCi:.2) 

2430 NEXT 
2432 CL0SE2 
2440 F0RI=1T08 

2442 PR INTI" "R$<D 
2444 NEXT 

2450 PRINT:R:T=:"":INPUT»IF CORRECT FILE PRESS V";R$ 

2452 IFR:T="V"THEN 2460 
2454 GOTO100 

2460 PR I NT: PRINT "FINISH I NO FETCH OF DATA FROM SOURCE DISK" 

2470 IFR^ < 1 )="0"THENIFR|:<3)»"84"THENUU-168 : GOTO2480 
2472 UU=VRL(Ai:<3)) 

2480 0PEN2 .. 8 . 2 . " 0 : " +K|:+ " . S . R " 

2482 F0RI=1T08 
2484 INPUT#2.MM$<I) 

2436 NEXT 

2488 FORI=lTOUU 

2490 INPUT#2.MR<I) 

2492 NEXT 
2494 CL0SE2 

2500 PRINT :PRINT"PLRCE DESTINATION DISK IN DRIVE" 

2308 PRINT :PRINT"RNY FILE NAMED # "Kl" M WILL BE ERASED" : PRINT 
2310 R$="": print: INPUT "TO CONTINUE INPUT V ";R« 

2511 IFR$=""THEN2400 
2530 0PEN15>8,15>"S:"+K$ 

2532 CLOSE 15 

2540 0PEN2. 8. 2. "0 : "+K$+" . S. W" 

2342 F0RI=1T08 
2344 PRINT#2,MMi:a) 

2546 NEXT 

2348 F0RI=1T0UU 

2350 PRINT#2.MRa) 

2552 NEXT 
2538 CL0SE2 

2560 PR I NT "SRVE OF SEQUENTIAL FILE IS COMPLETE" 

2563 PRINT:R:f="":INPUT"COPV ONTO ANOTHER DISK IF V";RI 
2366 IFR$="Y"THEN2500 
2580 GOTO100 



2600 REfl CHECK ROUTINE FOR INTERPOLATION ROUTINE 

2610 INPUT "ENTER NUMBER” ;XB 

2620 G0SUB1267 

2630 PRINT" I=" I 

2632 PRINT"CC<2/n="CC<2.n 

2634 PRINT"flNSWRE IS "BB 

2636 GOTO2600 


6000 FORI»1TO100 
6010 POKE57343.0 
6020 FORT-1TO300 
6030 NEXT 
6033 P0KE57343,2 
6040 FORT=1TO300 
6050 NEXT 
6060 GOTO6000 
6080 GOTO6000 


7000 INPUT"ENTER TOP RDG";TR 

7005 MM=0: INPUT "MULTIPLY BY"; MM 

7006 flT=TR+HM*TR 

7010 INPUT"ENTER BOT RDG";BR 

7011 MN*0:INPUT"MULTIPLY BY"; MM 

7012 flB=BR+MM>tfER 

7010 PRINT#4, "FIRST TOP = "TR" SECOND TOP = "FIT 
7019 PRINT#4. "FIRST BOT = "BR" SECOND TOP »= "AB 
7020' PR I NT#4, "RATIO = "RT/flB:RT»HT/AB 
7030 GOSUB1200 

7040 PRINT#4. "INDICATED TEMP= "TM 
7050 PRINT#4:PRINT#4 


8000 REM mmm get DET temp 
8002 GT=0:CH=5 

8004 GOSUB1120 

8005 PRINT"CTEMP SENS OUT * "R4*-l" >" 

8006 PRINT"DET. TEMP = "<INT < a000»iR4:^-l )-2731 . 5) )/10" DEG C" 
8008 GT=0:CH=6 

8010 GOSUB1120 

8012 PRINT"HTR POWER = "< 16X5)»R4*-1" WATTS" 

3013 PRINT 

8014 FORI=1TO1500:NEXT 
8016 GOTO8000 


30000 INPUT"ENTER NUMBER" ;UU 
30010 GOSUB1700 
30040 PRINT"Gf;E$= "G«;E$ 
30050 GOTO30000 


30100 INPUT "ENTER NUMBER" ;UU 
30110 GOSUB2000 

30140 print"g*;e$="g*;e$ 

30150 GOTO30100 

40000 OPEN 15,8,15:1 NPUT# 1 5 , A* , B* , C$ , D|: 
40010 PRINTA#,B»,C»,Di: 

40020 CL0SE15 
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"vr. 


1 GOTO 40100: REM B/B REflD-M10 

2 THIS LISTING "fl/D RERD-REM" 

3 MRCODE SEEMS TO BE WIPED OUT BV BASIC VRRIRBLES IF NOT PROTECTED 

4 COMPUTER RETURNS fl LONG HUMBER OCCRSIONALLV 

5 . 

7 . 

8 . 

9 . 

10 . 

11 . 

12 . 

14 . 

15 . 

16 . 

17 

18 . 

20 . LISTING OF GRINS 

21 . 

22 . OF VARIABLE GRIN AMPLIFIERS 

24 . 

25 . USED WITH DETECTORS 

26 . 

27 . 

23 . 

29 

30 . 

31 . CALCULATION OF COMBINATION 

32 . GAIN VALUES 

3*^ • 

34 *. TWO AMPLIFIERS IN SERIES 

35 

64 . Bf INDICATES DET. AND GAIN 

63 . 

66 . 

67 . 

68 . GF<X) IS AMPLIFIER GAIN SETTING . 

69 . 

70 . CF<X) IS OUTPUT PORT TO USE 

71 . . 

72 . DD<X) IS NUMERIC IDENTIFICATION . 

73 . OF DET £ GRIN - USED IN PRINTOUT . 

74 . 

75 . 

76 

80 CHECKS FLAG TO SEE IF MRCODE HRS BEEN LORDED. 

31 PRINTPRINTPRINTPRINTPRINTPRINTPRINTPRINTPRINTPRIHT 
86 GO TO CALIBRATION CURVE LOAD ROUTINE. 


1 00 


101 « ■* 

102 * MAIN MENU * 
104 » » 
106 ■* * 

103 « « 
110 * ■* 
112 * * 
114 •* * 
116 •* ■* 
118 * •* 
119 * * 


120 

122 CHECK FOR PROPER INPUT 
124 CHECK FOR PROPER INPUT 
126 MAIN MENU ADDRESSES mW'**m* 

160 DISPLRV ZERO READINGS. 

162 GOTO MENU. 

170 SET HARD COPY FLAGS. 

171 PRINT DATA WITHOUT HEADER. 

180 IHIT FLAGS. 

181 PRINT DATA WITH HERDER. 
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122 


200 

201 

202 

203 

204 

205 

206 
207 
203 
210 
211 
212 
213 
220 
222 
226 
223 
230 

232 

234 

235 

236 

233 

239 

240 

241 

242 

245 

246 

247 

243 

251 

252 
254 
256 

253 
260 
261 
262 

263 

264 

265 

266 
263 

269 

270 

272 

274 

275 

276 

273 
280 
299 


CLEAR SCREEN. 

STRRT OF CURVE LOAD ROUTINE. 
DEFAULT FILE NAME. 

DISPLAV REMINDER. 


. OPTION TO CHANGE NAME OR 
. TO BVPASS LOAD OF FILE 


GOTO MENU. 


GIVES TIME TO INSERT DISK INTO DRIVE. 

OPEN FILE. 

! ' ’ read’ eight' item* HERDER* ‘ ! 

• • 

close' file! 

VERIFV THAT FILE WAS FOUND - IF NOT FOUND THEN STRRT OVER. 
! DISPLAV FILE HEADER ! 


REOPEN FILE.' 


. REREAD FILE HEADER 


! READ FILE DATA 


CLOSE FILE* 

VARIABLE "‘bottom TEMP IN TABLE ‘M FIND IN HERDER LINE 4 
VARIABLE =■ TABLE TEMP INCREMENT * FIND IN HERDER LINE 4 
VARIABLE » TOP TEMP IN TABLE » FIND IN HEADER LINE 4 


. CALCULATE RATIO VALUES 

m 

! < BVPASS DISPLAV) 


DISPLAV LOWEST RATIO IN TABLE. 
DISPLAV HIGHEST RATIO IN TABLE. 
TIME DELAV FOR VIEWING OF ABOVE. 
GOTO MAIN MENU. 


300 CLEAR SCREEN. 

302 "MACODE" IS A MACHINE LANGUAGE PROGRRM THAT 

303 ENABLES USE OF I/O PORT RND CLOCK. 

323 THIS PROGRAM WAS SUPPLIED BV THE SUPPLVER OF THE PORT. 
330 . 

332 . 

333 

350 ! MMESSAGE IS NOT CORRECT. 

351 SET FLAG INDICATING THAT MACODE HAS BEEN LOADED. 

354 ! " " LOAD* RLE* HEADER ! 

355 . 

356 . 

357 . 

358 . (CHECK TO VERIFV THAT 

359 

360 ! (STARTING ADDRESS) ! 





362 

364 ! LORD CODE STRRTING RT ! 

363 . ABOVE RDDRESS. 

366 . 

368 

386 CHECK FLAG. 

3B7 GO TO MAIN MENU IF FLAG HOT« TO LD. 
390 ESCAPE IF DESIRED. 

392 TRV ANOTHER DISK. 


500 CLEAR SCREEN 

502 GET ZERO READINGS ON ALL RANGES. 

504 . 

530 * SET AMPLIFIER GRIN <RT PORT) 

532 ■* TIME DELRV TO LET AMPLIFIER SETTLE. 

533 * INIT. OUTPUT REG. 

534 * MAKE 12 READINGS ON EACH RANGE. 

535 * CHANNEL TO BE READ. 

536 * SSUBROUTINE READS OUTPUT OF DETECTOR. 

537 * DISCARD FIRST TWO READINGS. 

533 ■* KEEP SUM OF READINGS. 

540 * DISPLRV READINGS AS MADE. 

542 m 

546 « CALC. AVERAGE OF 10 READINGS AND STORE RESULT. 

547 » DISPLRV RESULT. 

550 END OF READ ROUTINE. 

552 

554 

555 .' DISPLRV SUMMRRV TABLE ! 

556 . OF ZERO READINGS. 

553 

567 . 

569 READ TIME FROM CLOCK ON I/O PORT. 

570 SET FLAG THAT INDICATES THAT ZERO READINGS HAVE BEEN TAKEN. 

571 STORE TIME OF END OF READINGS. 

572 . 

574 IF DOING OVER RESET FLAG. 

576 HARD COPV PRINT 

578 IF VES - HARD COPV ROUTINE. 

580 GOTO MAIN MENU 


600 . 


602 ****»:» READ INDIVIDUAL RANGE mmm 


604 

606 

610 

612 

614 

616 

618 

620 

621 


. DISPLRV SELECTION MENU 

! B$<X) IS DEFINED IN LINES 64 
. THROUGH 76 


622 

624 

625 
641. 
642 
660 
662 
666 
663 

669 

670 

671 

672 
680 
682 
683 


SELECT RANGE TO READ. 

IF OUT OF RANGE START OVER 
NOT NEEDED 

SETUP GRIN AND CHANNEL 
MAKE ONE READING 


. DISPLRV 

. RANGE - AMPS - VOLTS 
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700 . 

701 . 

702 WHEN ZERO RDGS. TRKEN FLAG IS SET AT ADDRESS 33142 

703 SELECT DISPLAV OR HARD COPV. 

704 SET FLAG 

703 IF HARD COPV NOT WANTED THEN READV TO READ. 

705 PRINT HARD COPV HEADER. 

707 READV TO READ ON COMMAND. 

703 . 

709 . 

712 DISPLAV RESULT. 

717 SELECT PROPER RANGE AND READ OUTPUTS. 

722 DISPLAV RESULT. 

724 PREVENT DIVISION BV ZERO. 

726 . 

723 a 

730 CALC. RATIO. 

732 ROUND OFF RATIO. 

744 DISPLAV RATIO. 

745 LOOK UP TAEMPERATURE - INTERPOLATION ROUTINE. 

747 LOOK UP TEMPERATURE INDICATED BV INDIVIDUAL DETECTORS. 

748 m AW$CX) ARE FILE HEADER ITEMS LORDED IN LINE 232 

749 . RW$<5) IS BOTTOM TEMP OF CAL. * 

754 . * 

755 . RWf(6) IS TOP TEMP OF CAL. itiCHECK THESE OUT#*****:^: 

762 . « 

764 . it! 

766 . AW$<3) IS TABLE LENGTH 

767 . RWf<4) IS 

768 . 

780 DISPLAV DATA. 

781 . 

783 FLRG=0 - DON'T PRINT 

784 PRINT DATA 

783 PROVISION FOR CONTINUOUS READOUT. HAVE TO EXIT PROGRAM TO SET WW 

786 MUST REENTER PROGRAM WITH A GOTO 700 OR GOTO 100 ETC. 

787 . 

783 RESTARTS READ CVCLE. 

790 HOT USED PRINT 
794 

796 

797 


800 - 
310 . 


812 . 

SET ALL ZERO 

814 . 


816 . 

READINGS 

820 . 


822 . 

TO ZERO 

824 a 


826 .. 



829 RESET ZERO READING FLAG. 

830 GOTO MAIN MENU. 


900 FIND AMPLIFIER GRIN ROUTINE 

902 SET GAIN AND AMP. COMBINATION TO MAX. GRIN. 

904 MAKE READING. 

906 IF OUTPUT GREATER THAN 9 THEN GO TO NEXT LESS SENSITIVE COMBINATION. 

903 IF OUTPUT LESS THAN 9 THEN CHECK BOTTOM DET. OUTPUT 
910 . 

912 . 

914 . 

916 . 

913 . 

920 . 

922 . 

924 . 

926 . 

923 . 

930 . 

932 . 
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nA 


■ 


936 

933 

940 

942 

944 

946 

947 

949 

950 

951 

952 
954 
956 
958 

970 

971 

972 

973 

974 

975 
977 
979. 
980 
982 
984 

990 

991 

993 

994 
993 

996 

997 


TIME DELflV. 

START - SELECT GAIN FOR BOTTOM DETECTOR. 
READ CLOCK - TIME OF START OF DATA READ. 


READ ONE, THEN THE OTHER DETECTOR AND REPEAT TEN TIMES. 

SET CHANNEL FOR BOTTOM DET. 

READ BOT. DETECTOR. 

IF AMP NEAR SATURATION - GOTO MENU. 

VB IS SUM OF OUTPUTS. 

SET CHANNEL FOR TOP DET. 

READ TOP DET. 

IF AMP NEAR SATURATION - GOTO MENU. 

VT IS SUM OF TOP DET. READINGS. 

GAIN SUBSCRIPTS IDENTIFY AMPLIFIERS AND GAIN SETTINGS. 


1002 DISPLAY AVERAGE OUTPUT OF TOP DET. 

1003 DISPLAY AVERAGE OUTPUT OF BOT. DET. 
1008 . 

1010 CALC. TOP DET. CURRENT. 

1020 CALC. EOT. DET. CURRENT. 

1030 . 


1100 START OF DET. READ ROUTINE, mm 

1102 INPUT TOP DET. CHANNEL. 

1104 . 

1110 INPUT BOT. DET. CHANNEL. 

1119 . 

1120 SET AMPLIFIER GRIN. 

1123 . 

1124 . 

1125 . 

1140 . 

1145 SET CHANNEL 

1150 . I/O PORT READS DET. OUTPUT . 

1155 . DISABLE DDRD REGISTER.. 

1160 . START A/D.. 

1165 . RESET START FLIP FLOP.. 

1170 . WRIT FOR DATA.. 

1172 . INIT VARIABLES 

1175 . ENABLE LOW BYTE ONTO BUSS.. . 

1180 . READ LOW BYTE.. 

1185 . ENABLE A/D HI BYTE.. 

1186 . READ HI BYTE.. 

1137 . INIT SIGN FLAG.. 

1183 . TEST FOR SIGN OF ANSWER.. 

1190 . GET RID OF HIGH NIBBLE.. 

1192 . CALCULATE SIGNED DECIMAL ANS. . 

1195 
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1200 

1201 . INTERPOLATION ROUTINE FOR 

1202 . RATIO TABLE. 

1204 . 

1206 . HUNTS FOR APPROXIMATE PART 
1208 . OF TABLE AND THEN DOES A 
1210 . STRAIGHT LINE INTERPOLATION. 
1212 . 

1214 .. 

1216 .**TAKES RATIO OF DETECTOR 
1218 . OUTPUTS AND LOOKS UP TENP.iK* 
1220 . 

1222 . 

1224 . 

1226 . 

1228 . 

1230 ....ROUNDOFF OF TEMPERATURE. 
1232 


1240 

1242 . 

1244 . LOOK UP IND.TEMP OF TOP DET. 

1245 . 

1246 . STRAIGHT LINE INTERPOLATION 

1247 . 

1248 . 

1249 . 

1250 . 

1252 . 

1254 . 

1256 . 

1258 . 

1259 . TT IS INDICATED TEMP OF TOP DET. 

1260 . 

1266 


1267 

1269 . 

1270 . LOOK UP IND. TEMP OF BOT. DET. 

1271 . STRAIGHT LINE INTERPOLATION. 

1272 . 

1273 . 

1274 . 

1275 . 

1276 . 

1277 . 

1278 . 

1230 . 

1282 . 

1284 


1400 READ CLOCK ROUTINE *stf»!K***» 

1420 

1422 ... ; 

1430 . CLOCK IS ON I/O PORT BOARD. . 
1432 . 

1436 . DATA MUST BE READ OUT AS 
1438 . INDICATED ON LINES 1430 THRU . 
1440 . 1442. 

1442 . 

1446 . 

1448 . DATA IS THEN CONVERTED TO ! 

1449 . STRING VALUES FOR USE BV THE . 

1450 . PROGRAM. 

1451 . 

1452 . 

1453 . 

1460 
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1500 *»** HARD COPV PRINT ROUTINE 

1504 SET VALUE OF UU TO TOP DET OUT. 

1505 ROUNDOFF UU 

150S ANSWER RETURNS AS Of AND AS EXPONENT E^ 

1507 NOT NEEDED. 

1508 SET VALUE OF UU TO DOT. DET OUT. 

1509 ROUNDOFF UU 

1510 ANSWER RETURNS AS C$ AND AS EXPONENT Ef 
1512 DELAY LOOP. 

1525 OPEN PRINTER TO PRINT DATA. 

1530 . FOLLOWING CHRf COMMANDS FOR S01 PRINTER. 

1531 . TAB<12) FROM START OF LINE ! ! ! 

1532 . TAB06) 

1533 . TAB<23) CHR*<16) INDICATES THAT NEXT CHRJ VALUES 

1534 . TAB04) ARE TAB VALUES RELATIVE TO START OF LINE. 
1536 . TAB<43) 

1544 . TAB<51) 

1546 . TAB<59) 

1548 . TAB<68) 

1550 CLOSE PRINTER. 

1560 . 


1565 

1566 

1567 

1568 

1569 

1571 

1572 

1573 

1574 
1576 
1578 
1585 
1590 


OPEN PRINTER TO PRINT HEADER. 


READ THE CLOCK (TIME OF DAY) 


CHECK ADDRESS TO SEE IF FLAG INDICATES ZERO RDGS TAKEN. 


CLOSE PRINTER. 


1700 *»i*i»*»** ROUNDOFF ROUTINE 

1701 . INIT VARIABLES 

1702 . CONVERT UU TO A STRING 

1712 . INIT 12 LINE TABLE 

1713 

1714 . FIND IF EXPONENT USED 

1715 . LOCATE DECIMAL POINT 

1716 . FIND LENGTH OF NUMBER 
1718 . 

1720 

1722 . 

1724 . E$=EXPONENT IE E-19 

1726 . ROUND LAST NUMBER UP OR DOWN AS NEEDED. 

1728 . 

1730 . GET 5 CHARACTER MANTISSA. 

1732 . ADD 6TH CHARACTER. 

1734 . PRINTPRINTPRINTPRINT 

1740 tttmt END OF ROUNDOFF fTTtttt 


1800 Wm-M OPEN PRINTER TO PRINT ZERO READINGS ***» 

1810 . 

1811 IF ZERO READINGS HAVE NOT BEEN TAKEN DO NOT CONTINUE. CLOSE. 

1812 . 

1820 . 

1830 . 

1850 . TAB<2) 

1851 . TAB<10) 

1852 . TABC25) 

1853 . TAB<41) 

1860 . 

1864 CLOSE. 

1890 tttn END OF PRINT ROUTINE TTtTT 
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2000 

2001 . 

2002 . NUMBER ROUNDOFF ROUTINE - 

2003 . 

2004 . DOES NOT RLUflVS WORK DUE TO 
2003 . 

2007 . STORRGE METHOD USED BV 
2003 . 

2010 . COMPUTER. 

2012 . 

2018 


2123 . 

2126 

2127 . 

2123 . ROUTINE INCOMPLETE - 
2130 . 

2132 . 

2133 . 

2134 . 

2136 . 

2140 . 

2160 . 

2162 


2400 . 

2402 copy FILE ROUTINE mmmMf. UTILITV ROUTINE m 

2410 . 

2412 IF HO ENTRV 00 TO MRIN MENU. 

2414 VERIFIES FILE NRME OF FILE TO BE LORDED FROM DISK 
2416 CHRHCE TO CHRNGE ERROR 
2418 . 

2420 . 

2421 . 

2424 OPEN FILE TO RERD FROM DISK. 

2426 READ 3 STRING FILE HERSER. 

2423 . 

2429 ...REMOVE FIRST CHARACTER OF STRING. 

2430 . 

2432 CLOSE FILE. 

2440 DISPLRV FILE HERDER. 

2442 . 

2444 . 

2450 ESCRPE OPTION. 

2453 

2454 GOTO MRIN MENU. 

2460 . 

2470 CHECK FOR FILE FORMAT RND ID. 

2472 FILE LENGTH. 

2430 DPEN FILE TO RERD FROM DISK. 

2482 REREAD FILE HERDER. 

2434 . 

2486 . 

2438 RERD DRTR. 

2490 . 

2492 . 

2494 CLOSE FILE. 

2500 CHRNGE DISKS IF DESIRED. 

2308 CONFIRMS FILE NRME. 

2510 ESCAPE OPTION. 

2511 . 

2330 OPEN COMMAND CHANNEL RND ERASE FILE OF SAME NAME. 

2332 CLOSE COMMAND CHANNEL. 

2540 OPEN FILE FOR OUTPUT TO DISK. 

2542 WRRITE FILE HERDER. 

2344 . 

2546 . 

2548 WRITE DRTR. UU IS NUMBER OF DRTR ITEMS. 

2530 . 

2332 . 

2358 CLOSE FILE. 

2560 . 

2363 . ■■ 

2566 . 

2380 GOTO MRIN MENU. 
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2600 

** START 

OF SOME TROUBLESHOOTING ROUTINES. 

2610 

• 



m 

2620 

• 



m 

2630 

• 



m 

2632 

• 



m 

2634 

• 



m 

2636 





6000 

• 



m 

6020 

■ 



m 

6030 

■ 



W-* 

6035 

• 



m 

6040 

■ 



m 

6050 

■ 



ww. 

6060 

■ 



m. 

6080 





7000 

• 



m 

7010 

• 



m 

7011 

• 



m 

7012 

• 



m 

7018 

• 



m 

7018 

• 



m 

7020 

• 



m 

7030 

■ 



m 

7040 

• 



m 

7050 





8000 

wmmm 

READ 

. , . 

TEMPERATURE OF DETECTOR W.mmW.* 

6002 

SET GRIH 

1 AND 

CHANNEL. 

8004 

READ 

SENSOR OUTPUT. 

8005 

DISPLAY 

TEMP 

IN DEG. K 

8006 

DISPLAY 

TEMP 

IN DEG. C 

8008 

SET GAIN AND 

CHANNEL TO READ HEATER POHER. 

8010 

mm 

READ HEATER POWER, 

8012 

DISPLAY 

POWER- 

8013 

• 




8014 

TINE 

DELAY TO READ NUMBER. 

8016 

READ 

TEMP. AND POHER AGAIN. 


30000 »» CHECK IHTERPOLHTIOH ROUTINE, m 
30010 . 

30040 DISPLHV RUSHER. 


30100 »* CHECK INTERPOLATION ROUTINE. 
30110 . 

30140 DISPLflV ANSWER. 

30150 TRV AGAIN. 

40000 mmm. read disk error data. ««««» 

40010 DISPLHV DISK STATUS. 



PROGRRM 


'B/B RERD-m0' 


5/ 17/83 


REMARKS 


1200 REM INTERPOLRTIOH ROUTINE 
1210 IFRT<CC(6, nTHENTn»0:PRINT"UT" :OOTO1270 
1212 IFRT>CC<S/84)THENTri*0:PRINT"0T" :GOTO1270 
1220 IFRKCCCS/ 15)THEHI=0:GOTO1256 
1222 IFRKCCC6, 30)THENI=14 : C0T01256 
1 224 I FRT <CC ( 6^ 45 > THENI =29 : GOTO 1256 
1226 IFRTCCCCS/ 60>THENI=44 : G0T01256 
1228 IFRTCCCCS/ 75)THENI=59 : G0T01256 
1 230 I FRTCCC C 6/ 84 5 THENI=74 : GOTOl 256 

1252 1=1+1 

1253 IFI>84THEN1266 

1254 GOTO1250 

1256 IFRT>CCC6/ I>THENIFRT<CC(6, I+1)THEH1262 

1257 1=1+1 

1258 IFI>B4THEN1266 

1259 GOTO 1256 

1262 RU=CRT-CC<6/ I)>/CCC<6/ I+1>-CCC6/ D) 

1263 TM=CCC5/ D+CCCCC5/ I+D-CCC5.. I))ifRU) 

1264 TM=CINTCTM*10 ))/10 

1265 GOTO1270 ^ 

1266 TN=0 
1270 RETURN 




^eirse, Stsarc. t? 






RERDV. 
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PROGRRfl 


•'fl/D REflD-M10'‘ 


5 ^ 17/85 


REMARKS 


3 


f 


1400 

1420 

1422 

1430 

1432 

1434 
1436 

1435 
1440 
1442 
1460 


REHDV 


REM mmm refid clock routine wmm 

SVS37513 uses *' M#S4<?oe.‘' to U>e>ri A.T CUOCK 


B=37888 

VR=PEEK<B+in+<PEEK<B+12>*10) 
f1O=PEEK<B+8)+<PEEK<B+105*10) 
DV=PEEK < B+7 > + < PEEK < E+8 ) * 1 0 ) 
Dfl=PEEK<B+6> 

HR=PEEK<B+4>+<<PEEK<B+5)-8):+:l0) 
MN=PEEK < B+2 > + < PEEK < B+3 ) * 1 0 ) 
SC=PEEK<B)+<PEEK<B+1 >*10) 
RETURN 


■N 

>■ euT<»t»T rK-'ST fte" RCAO IKI 
EveM IP orgtv The Tlr^e 

J 



THIS ORpB^ 
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PPsOGRRM 'fl/D REHD-m0' - 5/17/85 


REMARKS 


1500 

1505 

1530 

1531 

1532 

1533 

1534 
1536 

1535 
1550 
1560 


REfl PRINTOUT ROUTINE 
OPEN4,4.0 

PRINT#4iCHR$<16>CHRf<4S)CHRf<49) XT; i 

PRINT#4.CHR$a6>CHRS:<49>CHR^<55> QR; Tab 17 

PRINT#4,CHRf<16)CHR$<:50>CHR$(5O> XB; tab 21 

PRINT#4.CHR$<16>CHRf<51)CHRJ<55) LV; Ta6 37 

PRINT#4^CHR$<16>CHR$<52)CHR$<52> RT; tab 44 

PRINT#4/CHR$<163CHRf<53)CHRf<52? TM; TAB TT 

PRINT#4,CHR$a6)CHRf<54>CHR$<53)HR;riN;SC TAgtr 


CL0SE4 

RETURN 


Tsb valutf prr*Ntci 


ilist next 
kvJo vsloej are 
tit v»[oei I 


Ch 1 Zi(^t> 



etc. 


REflDV. 


1565 nPEN4^4,0 

1566 PRINT#4 

1570 PRINT#4^" TOP DIODE RANGE BOTTOM DIODE RANGE RATIO"; 

1571 PRINT#4^" IND. TEMP TIME" 

15S0 PRINT#4 

1535 PR1NT#4:CL0SE4 
1530 RETURN 
1533 STOP 


READV. 






PROGRAM 


•^fl/D RERD-M10 


5/17/85 


REMARKS 


1800 OPEH4,4,0:CMD4 

1810 PRIHT”LIST OF ZERO REfiDIHGS" -PRINT -PRINT-PRINT^PRINT 
1820 PRINT"r<"2Z"> RANGE VOLTS" 

1830 F0RI=1T012 

1832 IFI = 1THEN22=14:K^^="T0P LON$10X" : GOTO 1849 

1833 IFI=2THENZZ=4:K^^="T0P LOW" :G0T01849 

1834 IFI=3THENZZ=13:K$="T0P MED^fr: 1 ON GOTO 1849 

1835 IFI=4THENZ2=3:K^="T0P MED" : G0T01849 

1836 IFI=5THENZZ=12:K$="T0P H I if 1 ON GOTO 1849 

1837 IFI=6THENZZ=2NC$="T0P HI":G0T01849 

1838 IFI=7THENZZ=18:K$="B0T LOW^:10X" -GOTO 1849 

1839 IFI=8THENZZ=8UC$=:"B0T LOW" :G0T01849 

1840 IFI=9THENZZ=17:K$="B0T MED:^10X" J G0T01849 

1841 IFI“10THENZ2s7:KS="BOT MED" : G0TO1849 

1842 IFI=11THENZZ=16:K4:="B0T HliflON" ^ G0T01849 

1843 2Z=6:K^="B0T HI" 

1849 PRINT 

1050 PR I NTCHR:f C 1 6 ) CHR-^ < 48 > CHRS C 50 > ZZ ; 

1851 PR I NTCHR^ <16) CHR$ < 49 ) CHR$ < 48 ) K4: ; 

1352 PR I NTCHR$ < 1 6 ) CHR$ < 50 > CHR$ < 53 ) D < ZZ) 

I860 NEXT " s> " 

1864 PRINT#4JCLOSE4 
1890 RETURN 


REHDV. 



LIST OF 

ZERO RERDIHOS 


D< 6 ) 

RRNGE 

VOLTS 

14 

TOP LOW:f!10N 

.712402344 

4 

TOP LOW 

. 0805664063 

13 

TOP MED?H10N 

.0366210938 

3 

TOP MED 

.0146484375 

12 

TOP HI»10X 

.0170893433 

2 

TOP HI 

.0122070313 

13 

EOT LOW*10X 

.0366210933 

6 

EOT LOU 

.0170393433 

17 

EOT HED*10X 

9.765625E-03 

7 

EOT MED 

.0122070313 

16 

EOT HI*10X 

9.765625E-03 

6 

EOT HI 

.0122070313 


S'^\ug\e 


133 




SEQ. FILE 'MflCOIlE'' 


WHICH IS USED BV THE I/O PORT. 


THE FOLLOWING IS fl LISTING OF THE 
'HRCODE' DECiriflL ADDRESS AND 


DATA PORTION OF SEQ. FILE NAMED 
T»cn TMC31 np rnnp RRE GIVEN. 


37889 

-■ 

1 

37945 


223 

38001 


37890 

— 

2 

37946 

- 

234 

38002 

- 

37891 

" 

1 

37947 

- 

234 

38003 

- 

37892 

- 

4 

37948 

- 

234 

38004 

- 

37893 

— 

9 

37949- 

- 

234 

38005 

- 

37894 

— 

1 

37950 

- 

234 

38006 

— 

37895 


3 

37951 

- 

234 

38007 • 

« 

37896 


1 

37952 

- 

173 , 

38008 

— 

37897 


5 i 

37953 

- 

252 

38009 

- 

37898 

““ 

0 

37954 


223 

38010 

- 

37899 

•* 

^ ! 

37955 

- 

41 

38011 

- 

37900 


8 1 

37956 

- 

15 

38012 

- 

37901 

— 

255 ! 

37957 

- 

157 

38013 

- 

37902 


0 

37958 

- 

0 

38014 

- 

37903 

— 

72 

37959 

- 

148 

TSOI'S - 

— - 

37904 

• 

173 

37960 

- 

169 1 

38016 

- 

37905 

- 

254 

37961 

- 

10 t 

38017 

- 

37906 

- 

223 

37962 

- 

141 

38018 

- 

37907 

- 

41 

37963 

- 

252 

38019 

- 

37908 

- 

1 

37964 

- 

223 

38020 

- 

37909 

- 

240 

37965 

- 

138 

38021 

- 

37910 

- 

249 

37966 

- 

208 

38022 

- 

37911 

- 

104 

37967 

- 

211 

38023 

- 

37912 

- 

96 

37968 

- 

104 

33024 

- 

37913 

- 

72 

37969 

- 

170 

38025 

— 

37914 

- 

138 

37970 

- 

104 

38026 

- 

37915 

- 

72 

37971 

- 

96 

38027 

- 

37916 

- 

169 

37972 

*- 

72 

38028 

- 

37917 

- 

0 

37973 

- 

138 ; 

35029 

- 

37918 


141 

37974 

- 

72 i 

38030 

- 

37919 

- 

252 , 

37975 

- 

169 1 

1 33031 

- 

37920 

- 

223 

37976 

- 

0 1 

! 38032 

- 

37921 

- 

162 

37977 

- 

141 1 

i 38033 

i. 

37922 

- 

13 

37978 

- 

252 

! 33034 

- 

37923 

- 

202 

37979 

- 

223 

' 35035 

— 

37924 

- 

142 

37980 

- 

162 

; 38036 

- 

37925 

— 

253 

37981 


13 

: 38037 

- 

37926 

- 

223 

37932 

- 

202 

! 38038 

- 

37927 

- 

32 

i 37983 

- 

142 

! 35039 

— 

37928 

- 

15 

1 37934 

- 

253 

38040 


37929 

- 

148 

1 37985 

- 

223 

38041 


37930 

- 

169 

37936 

— 

32 

35042 

- 

37931 

— 

1 

37987 

- 

15 

38043 

— 

37932 

- 

141 

' 37988 

- 

148 

38044 


37933 

- 

252 

37989 

- 

169 

38045 


37934 

- 

223 

37990 

- 

1 

38046 

* 

37935 

- 

234 

37991 

— 

141 

38047 

— 

37936 


169 

37992 

— 

252 

38048 


37937 

- 

16 

37993 

— 

223 

38049 

— 

37938 

- 

141 

37994 

— 

234 

38050 

_ 

37939 

- 

252 

37995 

— 

169 

38051 


37940 

- 

223 

37996 

— 

0 

38052 

— 

37941 


169 

37997 

— 

141 

38053 


37942 

- 

19 

37998 

— 

252 

35054 


37943 

- 

141 

37999 

— 

223 

38055 

•• 

37944 

- 

252 

38000 

— 

189 

35056 

- 


0 

‘38057 


141 

38114 

- 

5 

148 

38058 

- 

161 

38115 

- 

141 

141 

35059 

- 

148 

38116 

- 

252 

253 

38060 

- 

72 

35117 

- 

223 

223 

35061 

— 

173 

35118 

- 

169 

169 

38062 

- 

160 

38119 

- 

0 

2 

38063 

- 

148 

38120 

- 

141 

141 

38064 

- 

41 

38121 

- 

252 

252 

35065 

- 

15 

38122 

- 

223 

223 

38066 

- 

141 

38123 

- 

96 

234 

38067 

- 

253 

38124 

- 

72 

234 

38068 

- 

223 

38125 

- 

169 

234 

38069 

“ 

32 

38126 

- 

20 

169 

38070 

- 

225 

38127 

- 

141 

0- 

33071 

- 

148 

38128 

- 

252 

141 

38072 

- 

173 

38129 

- 

223 

252 

38073 

- 

160 

38130 

- 

169 

223 

38074 

- 

148 

38131 

- 

16 

138 

38075 

- 

41 

: 33132 

- 

141 

208 

38076 

- 

240 

38133 

- 

252 

216 

38077 

- 

24 

* 38134 

- 

223 

169 

38075 

- 

106 

' 38135 

- 

104 

16 

38079 

- 

106 

38136 

- 

96 

141 

33080 

- 

106 

33137 

- 

255 

252 

•35081 

- 

106 

38138 

- 

0 

223 

i 35082 

- 

g 

33139 

- 

0 

104 

^ 38053 

- 

16 

38140 

- 

85 

170 

38084 

- 

141 

■ 35141 

- 

191 

104 

35085 

— 

253 

i 38142 

- 

Q 

96 

38036 

- 

223 

38143 


196 

n 

35087 

- 

32 





38088 

— 

225 

i 




38089 

- 

148 

1 



n 

35090 

- 

173 




n 

38091 

- 

161 




u 

38092 


148 

! 




35093 

- 

41 

t 



n 

38094 

~ 

15 

! 



0 

ct 

38095 

- 

Q 




0 

35096 

- 

32 

1 



ocrr 

38097 

- 

141 




n 

38095 

- 


1 



0 

rj 

38099 

- 

223 




0 

38100 

- 

32 





38101 

— 

225 





38102 

- 

148 




0 

38103 

- 

169 




0 

38104 

- 

48 





38105 

- 

141 




87 

38106 

- 

253 




169 

38107 

- 

223 




255 

35108 

- 

32 

! 



141 

38109 

- 

225 




160 

; 38110 

- 

148 




148 

' 38111 

- 

104 




1 69 

38112 

- 

96 




255 

38113 

- 

169 





THIS PROCRRH SUPPLIED BV CGRS MICROTECH - MFG. OF I/O PORT 
THIS LIST IS FOR TROUBLESHOOTING PURPOSES ONLY 
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PROGRAM •'PM-TEMP4R-' 


1 1/07/85 


PAGE 1 


OF 


10001 REM PROGRAM - PM-TEMP4R - TEMPERATURE CORRECTION ROUTINE 

10002 REM date: 11/Q7/S5 

10003 D I MAR* < 93) 

10004 REM CALC. ROUTINE ENDS AT LINE 104r0 

10005 PRINTCHR*<147) : PRINT: PRINT 

10006 print: PRINTTRBO) "PROGRAM * PM-TEMP4R PRINT: PRINT 

10007 PRINTTRB<5) "THERMOCOUPLE CORRECTION ROUTINE" 

10008 REM ** PRINTOUT ACCESSED FROM LINE 11152 m 

10009 FOR I = 1 TO300 : NEXT : PR I NT : PR I NT 

10010 PRINTTRBO) "ENTERING TABLE DATA " 

10011 REM 

10013 DIMC9<42/3) 

10014 F0RI3=1T03 

10015 F0RI2=1T042 

10016 RERDC9<I2. 13) 

10017 NEXTI2 

10018 NEXTI3 

10020 GOSUB 12001 

10021 PRINT:PRINT:PRINTTAB<9)"L0RDING DEFAULT DATA" 

10022 G0SUB12230 

10024 print: print: PRINT 

10023 PRINTTRB<2)"MI!J0THERWISE CALCULATION MILL BE MADETT’ 

10029 PRINTTRB<7); 

10030 INPUT"ENTER V TO EDIT INPUT ";ZZ* 

10031 IFZZ*>""THEN12100 

10091 PRINT CHR*a47):PRINT:PRINT:PRINT:PRINT:PRINT 
10097 PRINTCHR*<147) :PRINT:PRINT:PRINT 

10099 PRINT" CORRECTION FOR UNCORTED TC": PRINT 

10100 PRINT " PROGRAM IS NOW RUNNING" : PRINT " PLEASE WAIT FOR RESULTS" 

10101 REM UTCCRl IS ORIGINAL PROGRAM 

10102 PRINTJPRINT 
10142 IFF2=0THEN10155 

10145 PRINT" INDICATED TEMPERATURE WAS GIVEN" 

10150 GOTO10153 

10155 PRINT"GR3 TEMPERATURE WAS GIVEN" 

10153 F9=l 

10163 C6=. 15 

10165 IFF1=1THEN10172 

10167 PRINT" C6= "C6:PRINT 

10170 GOTO 10202 

10172 IFF9*2THEN10182 

10175 C7=l.ll 

10177 PRINT" C7= "C7"C0NSTANT MULTIPLIER" : PRINT 
10180 GOTO1O202 
10182 C5=i.477 

10184 PRINT" C5= "C5:PRINT 
10200 REM 

10202 IFC1S0THEN10210 
*10203 'IFM1=0THEN10213 

10204 IFK1«*0THEN10220 

10205 GOTO10235 

10210 K9al 

10211 GOSUB10780 

10212 ClsC8 

10213 GOT010203 

10215 K9=2 

10216 GOSUB 10780 

10217 MlsCa 

10218 GOTO 10204 

10220 K9=3 

10221 G0SUB10780 

10222 KlsC3 

10235 REM READ G1.T4,T3>W 

10240 N2=300*O1»D1/M1 

10245 N3=C1»M17K1 

10260 IFF1»1THEN10275 

10265 N4sC6#<N2t<2/3) )*<N3t< 1/3) ) 

10267 GOTO 10320 

10275 IFF9=2THEN10305 

10280 IFN2>5500THEN10295 

10285 N4=. 43+C7*'. 48*SaR<N2)#<N3t. 31 ) 

10290 GOTO10320 

10295 N4=.43+C7*.174*CN2t.618)#<N3t.31) 

10300 OOTO10320 



PROGRAM •'PM-TEMP4A-' 


- 11 / 07/85 


PAGE 


OF 


10305 H4»C5*SQR<N2)»<N3t.3) 

10320 H=12*<N4)»K1/Dl 

10330 Bl=. 173!)fF4*El#fll*Dl/<N4#Kliffl2#12) 

10340 H5=144itsN4»Kl/<K2*Rlt2) 

10355 IFF2=-1THEH10505 
10370 U9=T1 
10380 GOSUB 10742. 

10390 T6=U6 
10400 B8=U3 


rl0410 B9=V9 
10420 D7=U7 


10430 PRIHT"GflS FLOW = ”G1 

10431 PRINT" INDICATED TEMP.. = "T6 

10432 PRINT"GHS TEMP = "T1 

10433 PRINT"DELTR T - RAD... = "D8 

10434 PRINT"DELTA T - COHD. . = "D9 

10435 PRIHT"TQTAL DEL. T. . . .' = "D7 

10436 PRINT"H BTU/HR S FR . . = "H 

10437 PRINT"REVNOLDS NUM. . . . = "N2 
10450 PRINT 

10470 PRINT!PRINT:GOTO11150 


10505 S1*.173*F4#E1*A1/A2 
10508 21=F5»E2*<<T7/100)t4) 

10510 S2=a/F4).*<F4)K<<Tl/'100)t4)-<l-R3)i*iF6+:<(T4/100) T4)-Z1) 
10515 Xl=CSli»!S2/H)+Tl 
10520 P1=S1,^H 


10525 X2=1+.04»P1*<<X1/'100)13) 

10530 X3=<Tl+Pl*C<<Tl/100)t4)+3#<Oa/100)t45)VX2 

10535 X4*X1-X3+T1 

10545 Q6=L1#(SQR<:N5))/12 

10550 Q7= . 5# < EXP < Q6) + IXEXP < Q6) J 

10555 Q8=<T1-T3)/'Q7 

10560 X5=X4+Q8 


10570 U9=X5 
10575 GOSUB 10742 
10580 V5=U6 
10535 X6=X5+T1-V5 
10590 U9=X6 
10595 GOSUB10742 
10600 V6=U6 
“10605' 11=1 
10610 FORJ1=1TO10 

1 0620 X7= < X5-X6 ) * < T 1 -V6 ) X < V5-V6 ) +X6 
10625 U9=X7 


10630 GOSUB 10742 
10635 D7=U7 
10640 D8=U3 
10645 D9=V9 
10650 V7»U6 . 

10652 X8=T1+D7 
10655 S3=<V7-T1) 

10660 IFABS<S3K.02THEN10700 

10665 11=11+1 

10670 X5=X6 

10675 V5=V6 

10680 X6=X7 

10685 V6=V7 

10690 HEXTJl 

10700 PR I NT "GAS FLOW 

10701 PRINT" INDICATED TEMP.. 

10702 PRINT"GA3 TEMP 

10703 PRINT"DELTA T - RAD... 

10704 PRINT"TOTAL DEL. T 

10705 PRINT" ";H;N2;D9JD7 

10706 PRINT"H BTUXHR S FR . . 

10707 PR I NT "REYNOLDS NUM 

10720 PRINT 

10730 GOTO11150 


"G1 

"T1 

"X3 

"D3 

"D7 

"H 

"N2 


10742 B6=(K2*<Rl't2))/(Dl*12) 

1 0744 B7= < < N4i»iK 1 * 1 2 ) XD 1 + . 00692»F4* < < U9X 1 00 ) t3 > #E 1 ? 

10746 .Z1=F5*E2*<<T7X100) T4.'> 

10743 B8=. 173»El*<F6.+'<l-R3>*<<T4/100)'r4)-F4*<<U9X100M4>+21) 
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PROGRRM '■PM-TEMP4fl^ 


11/07/85 


PAGE 


OF 5 


10750 U5=U9+<B3/B7) 

10752 H6=<L1*SQR(B7/B5))/12 
10754 H7=. 5*<EXP<H6>+1/EXP<H6) ) 

10756 H8=<T3-U5)7N7 

10753 U6=U5+H3 
10760 U7SU9-IJ6 
10762 U8=U9-U5 
10764 V9=U5-U6 
10766 RETURN 

10730 IFTK400THEN10786 

10781 IFT1>4499THEN10736 

10782 T2=»INT<T1/100) 

10783 T9=<T1/100)-T2 

10784 C8=C3<< T2-3 > / K9 ) +T9>^ < C9 < < T2-2 ) ; K9 ) -C9 < < T2-3 > K9 ) ) 

10785 GOTO 10787 

10786 PRINT"T1MU3T BE BETNEEN 400 AND 4499 R" 


10787 RETURN 

10300 DATA. 2395. . 2397. . 2405. . 2417. . 2435. . 2458. . 2485. .2512 

10301 DATA .2542.. 2574.. 2607. .2638.. 2670.. 2700.. 2727.. 2752 

10802 DATA .2278..2S00..2823.. 2846.. 2868.. 2892.. 2915. .2933 

10803 DATA ,2962..2984..3010..3034..3062..3090..3120..3153 

10304 DATA ,3192.. 3237.. 3287.. 3342. .3414.. 3480.. 3588.. 3698 

10305 DATA .3328.. 3977 

10806 DATA .0356. .0418. .0478.. 0536. .0591.. 0643. .0692. .0733 
10307 DATA .0781.. 0820. .0861.. 0898.. 0934.. 0968.. 1002.. 1036 
10803 DATA .1072. .1100.. 1130.. 1158.. 1186.. 1214.. 1242.. 1270 
10309 DATA .1296..1325..1343..1373..1398..1422..1446..1470 

10810 DATA .1425.. 1520.. 1542.. 1568.. 1592. .1616.. 1640.. 1665 

10811 DATA .1690.. 1713 

10312 DATA .0118. .0143. .0168.. 0191.. 0213.. 0237.. 0258. .0280 

10313 DATA .0303. .0325.. 0347. .0368.. 0337.. 0407.. 0426.. 0446 

10314 DATA .0464.. 0481.. 0497.. 0515.. 0532.. 0545.. 0560.. 0573 

10315 DATA. 0537. .0601.. 0615.. 0628.. 0642.. 0655.. 0668.. 0632 

10316 DATA .0694.. 0707.. 0720.. 0732.. 0745.. 0758.. 0772.. 0784 

10317 "DFITA '.0798;. 0311 " 

11100 RETURN 


11150 

11151 

11152 

11153 

11154 
11232 
12001 

12004 

12005 

12006 
12007 


C1=VAL< AAi:<64) ) : ril=VAL<AA$<63) ) : K1=VALC AAf <72) ) 

REM RESET THESE SO THAT NEW VALUE WILL REFLECT INEXT INPUT. 
GOTO 12400 

REM PRINT ROUTINE AT 12400 
GOTO 10023 


D1 

FI 


RAf<l)=' 
RR«(5)=' 
HRf<:6)='' 
RA$<7)='' 
RR$<9)=”F4 


OD OF TC WIRE 

- FLOW DIRECTION 
<PERPENDICULRR =1)" 
(PARALLEL FLOW = 0)'' 

- TOTAL ANGLE FACTOR 


" : ART < 4 ) = " . 04 " : RR T ( 2 ) = '■ I NCHES " 
" :RRJCS)=". 1" 


’ :nRJ<12)=".35" :RRf a0>=" 


12003 fiR$<ll)=" " 

12010 HRf(13)="El - EMISSIVITV OF TC 

12011 RRf(14>='' ":RR$(15)=" " 

12013 RR*<17)="R1 - RADIATION AREA 

12014 RRfa8)=" ''••flA«a9)=" " 

12016 flR«<21)='"R2 - CONVECTION AREA 

12017 RRi:(22)a" "••flRf<23>=" " 

12019 flR$<25)a*'F2 - TVPE OF TEMP. INPUT 

12020 flRf(26)=" (GAS TEMP GIVEN = 0>" 

12021 RRf(27)=" (IND TEMP GIVEN =1)" 

12022 RR$(29)="K2 - THERM. COHD.- WIRE 

12023 RR$(30)=" (BTU/HR FT DEG. R )" 

12024 flR«(31)=" " 

12025 flRf(33)=''Rl - RADIUS OF WIRE 

12026 raj: ( 34 >*■' INCHES ":RR^<35>='' ” 

12023 RRf(37)="Ll - TC LENGTH TO SUPPORT 
12029 flR$(38)«''INCHES "!RA$(39)=’’ " 

12031 RR.r<41)*"R3 - RBSORBTIVITV OF GAS 

12032 RRfC42)=" " ••RR«<43>=” " 

12034 flR$(45)="F5 - FLAME ANGLE FACTOR 

12035 RRJ:(46)»" (RADIATION ANGLE FACTOR 


" :RA|:(16)=".3" 
"!RRf(20)="r' 
■•:RR$(24)="1" 
":RAf(23>="0" 

" :RA$(32)="14" 

":RA4:(36)=".020 

'•:AR$<40)=".70" 

"■•RR$(44)="0" 

":RR$(43)=".05" 


12036 flR|:(47)='' " 

12037 RRf(49)="E2 - FLAME EMISSIVITV " :flflt<52)="1.0" 

12033 Rflf(50)=" ":Rflf(51)=" " 

12040 RRf(53>="T7 - FLAME RAD. TEMP " :Rfl$(56)=“3960" 

12041 Rflf(54)=" (DEGREES R )" 

12042 RR'f(53)=" " 


:RR|:(3>=" 
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12043 

12044 

12045 

12046 

12047 
12043 

12049 

12050 

12052 

12053 

12055 

12056 

12057 
12053 

12059 

12060 
12061 
12062 

12064 

12065 
12067 
12063 

12069 

12070 
12100 
12101 
12102 

12103 

12104 

12105 

12106 
12107 
12103 
12110 
12111 
12112 

12113 

12114 

12115 

12116 
12117 
12113 

12119 

12120 
12130 
12140 

12230 

12231 

12232 

12233 

12234 

12235 

12236 

12237 
12233 

12239 

12240 

12241 

12242 

12243 

12244 

12245 

12246 

12247 
12243 

12249 

12250 

12251 

12252 
12260 
12300 

12302 
12304 
12306 

12303 
12310 


RR$<57>="F6 - NfiLL ANGLE FACTOR 
AA$<53)=" CRADIRTION FROM WALL) " 
RR$C59)=” " 

AR|:(61)«”C1 - SPECIFIC HEAT 
RRf(62)=*’ <BTU/# DEG, R > 

RR$<63)=” ” 

RR:rC65)="f11 - VI3C0SITV (#/FT HR) 
RA^(66)=” ":flA$<67) = ’* " 

RA$C69)="K1 - THERMAL COHD. 

RA^(70)=" ••:AR$<71)==" ” 

RR|:(73)-”Tl - GIVEN TEMP. 
flR$<74)=’* <TEMP INPUT DEG R)“ 
RR|:<75)=" '• 

RAi:C77)=''Gl - MASS FLOW 
RA:f<78)='' <#/SEC SQ FT )“ 

RR$<79)-" » 

RA^<81)=”T4 - WALL TEMP <DEG R) 
RR|:<32)=" *»:RA$<33)=" " 

RA^C35)=“T3 - SUPPORT TEMP 
AR|:<:86)=” •'••RRf(37) = " ” 
flR:f<S9)='’W - CONTINUATION CODE 

RR|:C90)=" <NOT USED IN THIS PROGRAM 
AR$C91)=” " 

RETURN 

PRINTCHR$< 147) J PRINT : PRim**WmmW3E 

PRINT 

GOTO 12300 

PRINT- PRINT jPRINTnO USE DESCRIPTIVE 

22$=‘"*:IHPUT"ENTER ID OF CHANGE"; 22$ 

IF22$="£"THEN10091 

IF22$>" "THEN12349 

IF22$=""THENI«r 

PRINT -PRINT 


" :AA$<60)=".S8" 

":RA$C64)="0" 

":RR$<68)="0“ 

":RA$<72)="0" 

":HR$<76)="2760" 

":AR$<80)=:"167'‘ 

" JRR$<:34)="2100" 

"*RR$C38)="21O0" 

" :RR$<92) = *’1" 

VERSION)" 

£ TO ESCAPE EDIT 

EDIT ROUTINE PRESS RETURN" 


PRINTCHR$a47):PRINT:PRINT"r*5^mUSE £ TO ESCAPE EDIT R0UTINE.^;m:.4::f:" 

print: PRINT 

F0RI=IT039STEP4 

PRINT" "RR$CI);RR$(I+3) 

PRINTAfl$a + l) 

PRINTRR$a+2) 

22$="0": INPUT "CHANGE TO ";22$ 

IF22$=:"£"THEN12130 
I F22$> " 0 " THENRflf < QQ+3 ) =22$ 
print: PRINT 
NEXT 

GOSUB12230 
GOTO 10023 
D1=VAL<RA$<4)) 

F1=VRL<RA$(8)) 

F4=VAL<AA$<12)) 

El=VRLCAR$a6)) 

Rl=VAL<Afl$<20)) 

R2=VAL<RA$<24)) 

F2=VRL<RA$<23)) 


K2=VRL<RR$<32)) 

Rl=VAL<Afi$<36)) 

L1=VRL<AR$(40)) 

R3=VRLCRR$C44)) 

F5=VRL<RR$(43)) 

E2=VRL<RRf(52)) 

T7=VRL(RR$<56)) 

F6=VRL<RR$<60)) 

Cl=VflL<RA$<64)) 

Ml=VRL<Rfl$(63)) 

K1=VALCAR$C72)) 

T1=VRL<RR$<76)) 

Gl=VRL<Rfl$(80)) 

T4=VRL(RA$<84)) 

T3=VRLCRR$C83)) 

W =VAL<RR$C92)) 

RETURN 

PRINT" 

PRINT :PRINT"D1= 
PRINT"E1= "El 
PRINT :PRINT"R1= 
PRINT"K2= "K2 
PRINT :PRINT"R1= 


\ 


DATA SUMMRRV" 

"Dl; :PRINTTAB<10)"F1= 

"Rl; :PRINTTRBa0)"R2= 

"Ri; :PRINTTRB<10)"L1= 


I 


"Fl; :PRINTTAB<20)"F4= 
"A2; :PRINTTRBC20)"F2= 
"Ll; :PRINTTAE(20)"A3= 


"F4 

••pc 

"A3; 
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PAGE 


- 11 /O'; 


•85 


12311 PRINTTflBC30)’*F5= »F5; :PRIHTTRBC70) "E2= "E2 
12315 PRINT‘*T7= "T7 
12318 PRINT :PRINT'T6= "F6 

12320 PRINT”C1= ”Ci; :PRINTTRB<10)”M1= "Ml ; :PRINTTRB(20)"K1= "fa; 

12321 PRINTTRB<30)"T1= "T1 

12324 PRIHT"G1= "GU :PRINTTRBa0>"T4= "T4; :PRINTTRB<20>"T3= "T3; 

12325 PRINTTRB<30)"W= ”W 

12326 GOTO12103 

12349 IF22$=""THEN 

12350 IF22.r="Dl"THEHI=?l. 

12351 IF22$="F1”THENI==5 

12352 IF22f-"F4*'THENI=9 

12353 IF22$="E1"THENI-13 

12354 IFZ25^="Rr*THENI=17 

12355 IFZ2T="R2"THENI=21 

12356 IFZ2^="F2"THENI=25 

12357 IFZZi:="K2”THEMI=29* 

12358 IFZ2I^=''R1" THEN 1=33 

12359 IF2Z$="L1”THENI=37 

12360 IFZZ5^="R3"THENI=41 

12361 IFZZT="F5"THENI=45 

12362 IFZZ$="E2"THENI=49 

12363 IFZZ$="T7"THENI=53 - 

12364 IFZZ$="F6"THENI=57 

12365 IFZZ$="C1"THENI=61 

12366 IFZ2^="rir’THENI=65 

1 2367 I FZ2J = " K 1 " THEN I =69 

12368 IF2Z$="T1"THENI=73 

12369 IFZZf="Gr*THENI=77 

12370 IFZZ:T="T4"THENI=S1 

12371 IF2Z$="T3"THENI=85 

12372 IFZZ:5="W"THENI=89 

12373 IFZZT="£"THEN12100 
12379, PR I NT -PRINT 

12380 PRINTRR$a)iRRfa+3) 

12382 PRIMTRfl$a+n jpRINTRR$<I+2) 

12383 QQ$=""-INPUT"ENTER NEW VRUJE ";QQ$ 

1 2384 I FQQ$> " " THENRR$ < I +3 ) =QQ^ 

12385 GOSUB12230 
12390 GOTO12100 
12400 OPEN4.4.0:CMD4 

12402 PRINT"DRTR OUT - PROGRRM m PM-TEriP4R m VERSION DRTED 11/07/85 



APPENDIX E - SENSOR STRUCTURAL AND OPTICAL ANALYSIS: 
CYLINDRICAL VERSUS CONICAL 


STRUCTURAL INTEGRITY OF THE SAPPHIRE CRYSTAL THERMOMETER 

A first-order analysis was made of the sensor in the typical gas turbine 
combustor exit environment. The sensing element must withstand gas pressure 
loading of the flowing hot gas stream, engine vibration, and occasional 
particle impact. For this first-order analysis, however, only the gas 
pressure loading will be considered. The best representation of the sensing 
element would be a cylinder in crossflow. A uniform velocity is assumed along 
the immersed length of the cylinder. A local Mach Number of 0.2 and a total 
pressure of 400 psi were taken as a worst-case steady load. The maximum 
steady-state stress was found by considering the sensor as a cantilever beam 
with uniform load. See references 18 and 19. 


Maximum Shear Stress , T max: 

4 V 

t max = r T 
3 A 

where A is the cross sectional area 

V is the shear force 
for a uniformly loaded beam, 

V = ql 

where q is the load per unit length 
1 is the length 

but q = CD X d X (p^ - P^) x (unit length) 

CD is the drag coefficient 

d is the diameter 

P^ is the total pressure 

P^ is the static pressure 

For a Mach Number of 0.2, (P. - P ) = 1.029 

’ t s 

and (P^ - P ) = 0.029 (400) 
t s 

q = 1.16 X d X (0.029 X 400) x (unit length) 
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For this loading, the maximum shear stress : 


t max = 4/3 V/A 


T max = 22.85 1/d (psi) 


The maximum tensile and compressive stresses are found from a circular beam 
equation: 


(a x) Max = 


32 M max 

n d^ 


^2 

and where M max, bending moment = 


(a x) Max = 68.55 (psi) 
d"^ 


Thus, both shear and tensile stress are dependent on the length-to-diameter 
ratio. Example 1/d values were chosen and stresses were calculated. This is 
shown in Table E-1. 


Table E-1. Stresses in a Cylindrical Sensing 
Element Due to Gas Flow Loading. 


Sensor 

Shear Stress 

Tensile Stress 

l/d 

psi 

psi 

5 

114 

1,714 

9 

206 

5,553 

17 

388 

19,811 


PHYSICAL PROPERTIES OF SAPPHIRE PRODUCTS 


Information channels were established with two sapphire products vendors, 
one university professor of materials science, and two GE materials engineers 
concerning the tensile strength and the modulus of rupture of single crystal 
aluminum oxide over the range of temperature from 0 to 1700® C. 
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What has been learned at this point is that sapphire exhibits a rupture 
strength that is dependent on both temperature and crystal orientation. In 
addition, sapphire exhibits a property of "slip" in various planes of the 
crystal structure. The slip results in deformation and is similar to "creep" 
in metal alloys. The stress required to cause slip ("flow stress") is 
dependent on both temperature and crystal orientation. At temperatures over 
about 800® C, the flow stress and orientation appear at this time to be con- 
trolling factors in sensor design. 

Referring to Table E-1, the tensile stress is dominant for the useful 
range of sensor 1/d (length to diameter ratio). The heat transfer studies 
showed that radiation loss is of greater importance than conduction loss for 
1/d of 4 or greater. 

The physical properties of single crystal sapphire were further investi- 
gated, and it was learned from Professor Heuer*s work at Case-Western Reserve 
University that slip along the basal plane at elevated temperatures would 
occur at a flow stress of approximately 700 to 800 psi at 1700® C. This 
orientation should be avoided. The prismatic slip plane is the second weakest 
plane, and would require approximately 14,000 psi to cause slip at 1700® C. 
This is encouraging. 


Cylindrical Element Stress Problems 

The length-to-diameter ratio for the sensor was considered. From the 
standpoint of keeping conduction loss errors small, a 1/d of 4 or greater 
is sufficient, as was mentioned in last month* s narrative. However, there is 
a real need to achieve an immersion distance sufficient to bring the sensor 
out of the boundary layer and into the **pitch line** region where a profile 
temperature peak is predicted. However, this is expected to increase the 
bending stress. 

As an example, for an immersion of 1 inch and a cylindrical sensor diame- 
ter of 0.040 inch, a 1/d of 25 would be needed (see Figure E-1). At the Mach 
0.2 and 400 psia flow condition, the tensile stress at the base of the cylin- 
der would reach 42,840 psi. This is considered too high for long life of the 
sensor. 


Cylindrical Element Light Pipe Problems 

It is still considered essential to maintain the light pipe properties 
of the cylindrical sensing element. The combustor discharge region is still 
expected to create solid deposits on a sensing element, just as on turbine 
parts. To prevent these deposits from causing attenuation of the internally 
reflected rays, an optical coating was proposed. This coating would maintain 
light pipe properties and prevent loss of calibration. 

Two major optical coatings manufacturers were contacted. One company 
has tested coatings to 900® C, but recommends a 600® C limit for good life. 
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Figure E-1, Cylindrical Sensor. 



6.25° 



Figure E-2. Conical Sensor. 



The second company did not expect any of its products to be applicable and 
suggested we contact a fiber optics manufacturer. 


Cone-Shaped Element as Solution to the Problem 

The idea occurred that a cone-shaped sensing element could solve two 
problems at once. A solid cone of 6.25 degrees included angle, with a tip 
truncated to 0.040 inch diameter, was selected for analysis purposes. The 
cone was presumed to be immersed 1 inch and held at the base (see Figure E-2) . 
The beam analysis was used to calculate the maximum tensile stress using the 
same gas loading conditions as were used for the cylindrical sensor. 

A tensile stress of only 1,016 psi was calculated. The beam analysis 
was checked by a different approach and appears to be correct. A shear stress 
of 76 psi was also calculated. This solves the first problem - stress. 

The second problem to be solved is that of maintaining light transmission 
from the emitting cavity at the sensor tip to the detector, even in the pres- 
ence of solid deposits. The cone shaped sensing element solves this problem 
by alleviating the need to depend on the light pipe mode for light transmis- 
sion. This is because the cone has enough of a solid angle to permit an 
imaging device to be used; that is, a lens. Figure E-3 shows how this could 
be accomplished. The lens, aperture, and spacing are such that only rays that 
emit from the tip and do not intersect the cone surface cannot reach the fiber 
optic cable face. Thus, foreign matter deposits will have no effect on the 
calibration. 
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APPENDIX G 


- ELECTRONIC DIAGRAMS. 
AND PANEL LAYOUT 
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Figure G-1. Front Panel Layout. 



piaERort'c 

IHPOT 

0^.0 

lO) 

0 0 


MANoAL| 
Auto « (\ ^) I 


o 


Fl/S £ 
/'/a AS ’’5 


cuTVrCT TePticc-i 50Trwi;c. KTf^.r'.A'r, 

“OPO;?CE. BOT.D'OD= 

C.AL!3l^/'’L^■ M x. lO x\0 A^'r^EK?, SPARL 

O 0 O 0 O 0 0 0 

]To\/ 

(_IW£ c^CP 


Figure G-2. Rear Panel Layout. 
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APPENDIX H - PHOTOCURRENT CALCULATIONS 


Calculations were made of photocurrent versus temperature for each of the 
elements of the two-color detector. The spectral sensitivity characteristics 
of the detectors were taken from the spec sheets furnished by the vendor, 
EG&G. Planck's law was used in a computer program that divides the silicon 
devices' spectrum into 43 wavelength intervals from 0.300 to 1.14 micrometers. 
The detected energy in each of the 0.02 micrometer bands was summed to fom 
the total response. Losses for each of the optic elements in the sensor 
system were estimated and taken into account. Table H-1, attached, is a 
condensed tabulation of the results . Column 1 is the temperature in degrees 
C. Column 2 is the photocurrent for the bottom diode element. The third 
column is the photocurrent for the bottom diode element. And the fourth 
column lists the ratio of top diode element to bottom diode element photo 
currents . 


Table H-1. Photocurrent and Ratio Versus Temperature. 


Temperature 
Degrees C 

Photocurrent, Top 
Diode Element Amps 

Photoccurrent , Bot . 
Diode Element Amps 

Current Ratio 
I-Top/I-Bottom 

500 

1.6609E-10 

6.0021E-11 

2.767 

700 

8.1257E-09 

2.2862E-09 

3.554 

900 

1.1207E-07 

2.5405E-08 

4.411 

1100 

7.4998E-07 

1.4063E-07 

5.333 

1300 

3.1908E-06 

5.0515E-07 

6.317 

1500 

1.0026E-05 

1.2621E-06 

7.360 

1700 

2.5454E-05 

3.0077E-06 

8.463 

1900 

5.5296E-05 

5.7463E-06 

9.622 



APPENDIX I - OPERATING INSTRUCTIONS 


1. Assemble probe and fiber optic cable. See Figure 12, "Optical 
Gas Temperature Sensor and Holder Assembly" and Appendix F, 
"Mechanical Detail of Sensor Assembly." 

2. Assemble electronics system. See Figure 14, "1700° C Gas 

Temperature System - Electronics Block Diagram." 

3. Apply line power: 115 volts alternating current, 60 hertz to 

the fully connected system. Give the thermostated detector 
module about 15 minutes for temperature stabilization before 
taking readings. 

4. To take data, load the program labeled "A/D READ-MIO" and run 
it. When run, the program will request a calibration file. 
Select the file corresponding to the probe and fiber optic 
cable being used. 

5. Shield the probe from stray energy and from elevated temper- 
ature and take a "zero-reading set" according to the program 
menu. 

6. Take sensor temperature data by selecting from the menu. The 
sensor temperature data is displayed on the screen; the data 
can be printed if the print feature is selected. * 

7. The gas temperature can be calculated from the sensor temper- 
ature if the convection and radiation heat transfer environment 
is entered. Load program PM-TEMP4A and enter the data from the 
keyboard as requested. 
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APPENDIX J - SCANNING ELECTRON MICROGRAPHS 


A set of scanning electron micrographs are presented for Cone Elements 
1 and 3, taken after the oven and propane-air flames testing given in Table 
11, Test Summary Chart. Figures J-1 and J-2 are similar views of the Cone 
Element 1 tip. The solid sapphire cone apparently had a very small tip frac- 
ture during manufacture of about 80 to 85 microns diameter. This small 
truncation was coated, as was the forward third of the cone. The sputtered 
70% platinum and 30% rhodium coating conformed to the striated fracture 
surface as well as the fine-ground cone side. The thickness was estimated by 
the vendor at the time of the coating to be 2000 Angstroms. 

This was followed by an aluminum oxide topcoat of 1.5 microns (estimated 
by vendor). Figures J-1 and J-2 show excellent adherence with little or no 
chipping or spalling of the coatings. 

Cone Element 3 is shown in Figures J-3, J-4, J-5, and J-6. This cone 
element had a 3500 Angstrom coating of 94% platinum and 6% rhodium with no 
overcoat. These figures, just as the previous figures, were taken after high 
temperature exposure. However, as shown on Table 11, a peak temperature of 
1750® C was reached, and 107 cycles from room temperature to 1700® C were also 
run. Figures J-3 and J-4 each show some debris attached. This is believed to 
be dust that was statically attracted between test time and SEM time. The 
photos show good adhesion of the metallic coating. Figures J-5 and J-6 were 
taken at the temperature transition region about one-half inch from the tip. 
Figure J-5 shows a darkened appearance at the top of the figure where the 
temperature did not reach as high a value as the lower half of the figure. 
The peak temperature at the transition region was at least 100® C cooler than 
the tip. Comparing Figures J-4 and J-6 leads one to conclude that the coating 
was flame polished near its melting point at the tip. Figure J-4 shows a much 
smoother surface texture than Figure J-6. 
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Figure J-4. 


Cone Element of No. 3 Tip. 
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Figure J-5 


Cone Element of No. 3 Side. 



Figure J-6. Cone Element of No. 3 Side 
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